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space  between  two  surfaces  that  are  (a)  separated  or  (b)  contact  directly. 

This  objective  was  achieved  by  deriving  microstructural  contact  laws  for  the  combination  of  two  cemented  elastic 
spherical  grains  subject  to  normal,  shear,  or  torsional  loading.  The  cement  was  treated  as  an  elastic,  elastic-plastic, 
purely  plastic,  and  viscous  material.  The  macroscopic  constitutive  laws  were  derived  from  the  above  microstructural 
contact  laws  by  using  the  formulas  which  relate  the  effective  moduli  of  the  random  packing  of  identical  spheres  to 
the  normal  and  tangential  contact  stiffnesses.  The  theoretical  solutions  were  positively  supported  by  experiments 
performed  on  artificial  cemented  aggregates,  as  well  as  by  othCT  experiments  performed  of  real  cemented  rocks. 

The  main  relevance  to  the  Air  Force  mission  is  through  a  quantitative  description  of  cemented  geomaterials  such  as 
asphalt  cement.  The  results  have  been  used  at  the  Wright  Laboratory  to  model  asphalt  concrete*s  behavior. 
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1.  SUMMARY 


1.1.  OBJECTIVE 

The  overall  objective  of  the  proposed  research  was  to  provide  a  quantitative  description  of 
the  microstructural,  as  well  as  macroscopic  mechanical  behavior  of  particulate  materials  with 
intergranular  cementation  in  a  wide  range  of  strain  amplitudes  and  for  various  types  of  the 
cementation  material.  By  cementation,  a  material  is  meant  which  fills  the  space  between  two 
surfaces  that  are  (a)  separated  or  (b)  contact  directly  (Figure  1.1). 

This  objective  was  achieved  by  deriving  microstructural  contact  laws  for  the  combination  of 
two  cemented  elastic  spherical  grains  subject  to  normal,  shear,  or  torsional  loading.  The 
cement  was  treated  as  an  elastic,  elastic -plastic,  purely  plastic,  and  viscous  material.  The 
macroscopic  constitutive  laws  were  derived  from  the  above  microstructural  contact  laws  by 
using  the  formulas  which  relate  the  effective  moduli  of  the  random  packing  of  identical  spheres 
to  the  normal  and  tangential  contact  stiffnesses. 


Figure  1.1.  Left;  Separated  grains,  grains  with  a  point  contact,  and  precompacted  grains.  Right;  Normal, 
shear,  and  torsional  deformation  of  two  cemented  grains. 

The  major  thrusts  of  the  research  were  in  the  theoretical  solutions  to  the  problems  of 
deformation  of  two  cemented  elastic  spherical  particles  subject  to  loadings  of  various  types  and 
amplitudes,  as  well  as  for  various  mechanical  models  of  the  cement  material. 

1.2.  Accomplishments 

Micromechanics.  The  micromechanical  solutions  deal  with  the  interaction  of  two  elastic 
particles  which  are  cemented  at  their  contact.  Typically,  the  contacting  surfaces  are  assumed  to 
be  spherical.  However,  the  solutions  obtained  are  general,  and  can  be  used  for  any  non- 
conforming  surfaces.  The  main  assumption  is  that  the  contact  area  is  small,  as  compared  to  the 
particle  size,  so  that  the  particle  can  be  treated  as  an  elastic  half-space  (in  the  three-dimensional 
case)  or  as  a  half-plane  (in  the  two-dimensional  case).  As  a  result  of  the  solutions  were 
obtained:  stress  distributions  at  the  contact,  and  normal,  shear,  and  torsional  stiffnesses  of  a 
two-grain  combination.  In  all  cases  the  solution  was  reduced  to  an  ordinary  integral  equation 
which  had  to  be  solved  numerically.  The  following  problems  were  solved: 

•  Normal,  shear,  and  torsional  deformation  of  two  initially  undeformed  elastic  cemented 
grains  (3-D  and  2-D)  with  elastic  cement. 

•  Easy-to-use  formulas  for  the  normal  and  shear  stiffnesses  of  two  initially  undeformed  elastic 
cemented  grains  (3-D)  with  elastic  cement.  These  formulas  have  a  closed  analytical  form. 

•  Normal,  shear,  and  torsional  deformation  of  two  precompacted  elastic  cemented  grains  (3-D 
and  2-D)  with  elastic  cement.  In  this  case,  the  grains  were  first  normally  pressed  together, 
and  the  cement  was  deposited  afterwards. 
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•  Normal  and  shear  deformation  of  two  initially  undeformed  elastic  grains  with  elastic-plastic 
cement. 

•  Normal  deformation  of  two  initially  undeformed  elastic  grains  with  perfectly  plastic  cement 
(cement  squeezing). 

•  Normal  deformation  of  two  initially  undeformed  elastic  grains  with  viscous  cement  (a 
hydroelasticity  problem). 

•  Normal  deformation  of  two  initially  undeformed  elastic  grains  where  cement  is  elastic  and 
evenly  covers  the  surface  of  spherical  grains.  When  a  load  is  applied,  a  new  contact  surface 
develops  between  the  shells  of  cement.  This  solution  allows  for  introducing  pressure- 
dependence  in  the  elastic  properties  of  cemented  aggregates. 

Macroscopic  Laws.  The  micromechanical  solutions  were  used  to  obtain  constitutive 
macroscopic  laws  for  the  deformation  of  cemented  particulate  materials.  It  was  assumed  that 
the  particulate  aggregate  is  a  random  dense  pack  of  identical  elastic  spheres  with  porosity  0.36 
and  coordination  number  9.  The  following  constitutive  laws  were  derived: 

•  Small-strain  effective  elastic  moduli  for  the  initially  undeformed  aggregate  with  elastic 
cement. 

•  Small-strain  effective  elastic  moduli  for  the  precompacted  aggregate  with  elastic  cement. 

•  Small-strain  effective  elastic  moduli  for  the  initially  undeformed  aggregate  with  elastic 
cement  where  pressure  dependence  is  taken  into  account. 

•  Small-strain  effective  elastic  (complex)  moduli  for  the  initially  undeformed  aggregate  with 
viscous  cement.  The  solution  provides  for  calculating  velocity-frequency  dispersion  and 
attenuation  of  acoustic  waves. 

•  Large-strain  stress-strain  laws  for  the  aggregate  with  elastic-plastic  cement  for  various  types 
of  loading  (e.g.,  hydrostatic,  uniaxial). 

•  Large-strain  stress-strain  law  for  the  hydrostatic  compaction  of  the  aggregate  with  perfectly 
plastic  cement. 

Experiments.  Elastic  wave  velocities  were  measured  on: 

•  epoxy-cemented  glass-bead  samples; 

•  ice-cemented  glass-bead  and  sand  samples; 

•  glass-bead  samples  with  viscous  cement  (uncured  epoxy). 

Large-strain  and  failure  experiments  were  performed  on  epoxy-cemented  glass  samples. 
The  strength  was  monitored,  and  the  large-strain  bulk  modulus  was  calculated. 

Verification.  The  theoretical  solutions  were  positively  supported  by  the  above- 
mentioned  experiments  performed  on  artificial  cemented  aggregates,  as  well  as  by  other 
experiments  performed  of  real  cemented  rocks.  Experiments  on  visualizing  stress  distribution 
in  contacting  cemented  elastic  cylinders  were  performed  at  the  University  of  Rhode  Island. 
These  experiment  supported  the  theoretical  predictions  as  well. 

Relevance  and  Application.  The  main  relevance  to  the  Air  Force  mission  is  through  a 
quantitative  description  of  cemented  geomaterials  such  as  asphalt  cement.  Our  colleagues  at 
Wright  Laboratory  used  our  solutions  to  model  asphalt  concrete's  behavior. 

1.3.  Personnel  Supported 

Senior  Personnel:  Dvorkin,  Jack,  PI;  Nur,  Amos,  Co-PI;  Mavko,  Gary  (all  US). 
Postdocs:  Yin,  Hezhu  (non-US),  Ph.D.  in  Geophysics,  Stanford,  1992. 

Students:  Frank  Liu  (non-US),  Ph.D.  in  Geophysics,  Stanford,  1994;  Maureen  Jacoby 

(US),  M.S.  in  Geophysics,  Stanford,  1994;  Packwood,  James  (US),  Ph.D.  candidate. 

1.4.  PUBLICATIONS 

1 .  Dvorkin,  J.,  and  Nur,  A.,  1993,  Rock  Physics  for  Characterization  of  Gas  Hydrates,  in 
The  Future  of  Energy  Gases,  USGS  Publications,  293-298. 

2.  Dvorkin,  J.,  Nur,  A.,  and  Yin,  H.,  1994,  Effective  Properties  of  Cemented  Granular 
Materials,  Mechanics  of  Materials,  18,  351-366. 
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3.  Yin,  H.,  and  Dvorkin,  J.,  1994,  Strength  of  Cemented  Grains,  Geophysical  Research 
Letters,  21,  903-906. 

4.  Dvorkin,  J.,  Yin,  H.,  and  Knight,  R.,  1994,  Seismic  Detection  of  Residual 
Contaminants,  64  Annual  Meeting  of  the  Society  of  Explorational  Geophysicists, 
Expanded  Abstracts,  Los  Angeles,  584-586. 

5.  Dvorkin,  J.,  and  Yin,  H.,  1995,  Contact  Laws  for  Cemented  Grains:  Implications  for 
Grain  and  cement  Failure,  Int.  J.  Sol.  Struct.,  32,  2497-25 10. 

6.  Jacoby,  M.,  Dvorkin,  J.,  and  Liu,  F.,  1996,  Elasticity  of  Partially  Saturated  Frozen  Sand, 
Geophysics,  61,  288-293. 

7.  Dvorkin,  J.,  and  Nur,  A.,  1996,  Elasticity  of  High-Porosity  Sandstones:  Theory  for  Two 
North  Sea  Datasets,  Geophysics,  in  press. 

8.  Elata,  D.,  and  Dvorkin,  J.,  1996,  Pressure  Sensitivity  of  Cemented  Granular  Materials, 
Mechanics  of  Materials,  in  press. 

9.  Dvorkin,  J.,  1996,  Large  Strains  in  Cemented  Granular  Aggregates:  Elastic-Plastic 
Cement,  Mechanics  of  Materials,  in  press. 

10.  Sienkiewicz,  F.,  Shukla,  A.,  Sadd,  M.,  Zhang,  Z.,  and  Dvorkin,  J.,  1996,  A  Combined 
Experimental  and  Numerical  Scheme  for  the  Determination  of  Contact  Loads  Between 
Cemented  Particles,  Mechanics  of  Materials,  22, 43-50. 

11.  Dvorkin,  J.,  and  Yale,  D.,  1996,  Plastic  Compaction  of  Cemented  Granular  Materials, 
Computers  and  Geotechnics,  in  press. 

1.5.  Interactions  and  Transitions 

Workshops  and  Conferences: 

•  Society  of  Explorational  Geophysicists  Development  and  Production  Forum,  Missouri,  June 
1993:  Role  of  Cement  in  Rocks 

•  American  Geophysical  Union  1993  Annual  Meeting,  San  Francisco,  December  1993: 
Strength  of  Cemented  Grains 

•  Caltech,  Symposium  on  the  Dynamic  Failure  Mechanics  of  Modern  Materials,  February 
1994:  Failure  and  Strength  of  Cemented  Grains 

•  UC  San  Diego,  Mechanics  and  Statistical  Physics  of  Particulate  Materials,  June  1994: 
Cement  among  Grains 

•  SEG  1994  Annual  Conference,  Los  Angeles,  October  1994:  Seismic  Detection  of  Residual 
Contaminants 

•  Phillips  Oil  Company,  and  Conoeo  Oil  Company,  Oklahoma,  January  1995:  Effect  of 
Cementation  and  Diagenesis  on  Physical  Properties  of  Rocks 

•  Symposium  on  Mechanics  of  Effective  Media,  University  of  Karlsruhe,  January  1995: 
Cemented  Geomaterials 

•  The  Bay  Area  Geophysical  Society,  Invited  Lecture,  San  Ramon,  CA,  March  1996: 
Elasticity  of  High-Porosity  Sands 

Transitions: 

•  USAF  Laboratories:  Wright  Laboratory,  Tyndall  AFB,  Dr.  Jeff  Rish  and  Dr.  Han  Zhu  — 
used  the  cementation  theory  (viscoelastic  cement)  to  describe  the  deformation  of  asphalt 
concrete. 

•  Universities:  The  University  of  Rhode  Island,  Dr.  Arun  Shukla:  Theoretical  and 
experimental  (photoelasticity)  investigation  of  the  dynamic  interaction  of  cemented  disks. 

1.6.  PATENTS  AND  INVENTIONS 
No  patents  or  inventions  resulted  form  this  effort. 
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2.  TECHNICAL  PART 


In  this  part  the  technical  results  are  presented  in  the  form  of  journal 
ARTICLES.  The  numeration  of  formulas  and  figures  is  internal  for  every 

SECTION. 


2.1.  Effective  properties  Of  Cemented  Granular 

MATERIALS 

ABSTRACT 

An  analytical  model  is  developed  to  describe  the  effective  elastic  properties  of  a 
cemented  granular  material  that  is  modeled  as  a  random  packing  of  identical  spheres.  The 
elastic  moduli  of  grains  may  differ  from  those  of  cement.  The  effective  bulk  and  shear 
moduli  of  the  packing  are  calculated  from  geometrical  parameters  (the  average  number  of 
contacts  per  sphere  and  porosity),  and  from  the  normal  and  tangential  stiffnesses  of  a  two- 
grain  combination.  The  latter  are  found  by  solving  the  problems  of  normal  and  tangential 
deformation  of  two  elastic  spherical  grains  cemented  at  their  contact.  A  thin  cement  layer  is 
approximated  by  an  elastic  foundation,  and  the  grain-cement  interaction  problems  are 
reduced  to  linear  integral  equations.  The  solution  reveals  a  peculiar  distribution  pattern  of 
normal  and  shear  stresses  at  the  cemented  grain  contacts:  the  stresses  are  maximum  at  the 
center  of  the  contact  region  when  the  cement  is  soft  relative  to  the  grain,  and  are  maximum 
at  the  periphery  of  the  contact  region  when  the  cement  is  stiff.  Stress  distribution  shape 
gradually  varies  between  these  two  extremes  as  the  cement’s  stiffness  increases.  The 
solution  shows  that  it  is  mainly  the  amount  of  cement  that  influences  the  effective  elastic 
properties  of  cemented  granular  materials.  The  radius  of  the  cement  layer  affects  the 
stiffness  of  a  granular  assembly  much  more  strongly  than  the  stiffness  of  the  cement  does. 
This  theoretical  model  is  supported  by  experimental  results. 

INTRODUCTION 

The  mechanical  characteristics  of  geologically-derived  granular  materials  such  as  rock, 
soil,  and  concrete  can  be  strongly  affected  by  the  properties  and  position  of  intergranular 
bond  material.  Microstructural  contact  laws  of  grain-to-grain  interaction  provide  an 
indispensable  basis  to  theoretical,  as  well  as  to  numerical  models  of  deformation  of 
granular  media. 

In  spite  of  significant  progress  achieved  in  modeling  macroscopic  deformation  of 
granular  media  (e.g.,  Cundall  and  Strack,  1979;  Nemat-Nasser,  1983;  Nelson  et  al.,  1988; 
Cundall  et  al.,  1989),  the  models  of  microscopic  intergranular  interaction  require  further 
development.  The  main  challenge  is  in  predicting  contact  stiffness  and  strength  from  the 
geometry  of  a  contact  region  and  from  the  mechanical  properties  of  grains  and  intergranular 
cement.  All  rigorous  theoretical  investigations  of  deformation  of  aggregates  of  elastic 
grains  are  based  on  the  classical  solutions  to  the  problems  of  normal  (Hertz,  1882)  or 
oblique  (Mindlin,  1949;  Walton,  1978)  interaction  of  elastic  spheres.  Digby  (1981)  solved 
the  problem  of  interaction  of  spherical  particles  that  are  initially  bonded  together  across 
small  areas.  These  solutions  resulted  in  a  number  of  theoretical  and  numerical  works  that 
examine  low-strain,  as  well  as  high-strain  deformation  of  granular  media  (the  latter  models 
include  relative  sliding  of  grains).  Many  of  these  works  are  discussed  in  Stoll  (1989)  in 
relation  to  the  acoustic  properties  of  sediments. 

Only  recently,  the  effect  of  intergranular  cementation  on  the  elastic  as  well  as  inelastic 
behavior  of  particulate  materials  started  to  attract  the  attention  of  researchers.  Bruno  and 
Nelson  (1991)  used  a  two-dimensional  discrete  element  procedure  to  investigate  rock 
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failure  under  tensile,  uniaxial  compression,  and  biaxial  loading.  In  the  elastic  domain, 
contact  stiffnesses  were  assumed  to  be  linear  functions  of  the  Young’s  and  shear  moduli  of 
the  cement,  and  the  thickness  and  width  of  cementation  bonds.  An  experimental  study  by 
Bernabe  et  al.  (1992),  conducted  on  synthetic  cemented  granular  materials  (precompacted 
Ottawa  sand  with  halite  and  silica  glass  as  the  cement),  included  triaxial  compression  tests. 
It  has  been  found  that  a  small  amount  of  cement  can  significantly  increase  the  strength  of  a 
granular  material  if  it  is  precisely  deposited  at  previously  formed  grain-to-grain  contacts. 

Trent  (1989),  and  Trent  and  Margolin  (1992)  numerically  investigated  the  behavior  of 
cemented  granular  materials  under  both  low-strain  and  high-strain  loads.  The  materials  are 
composed  of  circular  particles  which  are  glued  together  with  elastic  bonds.  They  show  that 
the  effective  properties  of  samples  are  governed  by  the  properties  and  distribution  of 
individual  intergranular  bonds. 

The  impressive  progress  of  numerical  methods  combined  with  the  increasing  power  of 
computers  and  computer  visualization  (Williams  and  Mustoe,  1993)  challenges  us  to 
fundamentally  understand  and  theoretically  describe  interparticle  contact  laws.  Indeed,  it  is 
important  to  know  how  stresses  are  transmitted  among  deformable  particles  through 
deformable  interparticle  bonds.  Two  cemented  grains  have  to  be  necessarily  considered 
deformable  when,  for  example,  they  have  a  point  contact  with  cement  deposited  around  it. 
If  these  grains  are  absolutely  rigid,  the  two-particle  combination  will  not  deform  at  all.  In 
this  case  there  is  no  way  of  finding  the  contact  stiffnesses  by  examining  the  deformation  of 
the  cement  alone. 

Dvorkin  et  al.  (1991),  examined  the  normal  interaction  of  two  spherical  elastic  grains 
and  an  elastic  cementation  layer  between  them.  Their  results  indicate  that  a  thin  cement 
layer  subject  to  normal  and  shear  load  can  be  approximately  treated  as  an  elastic  foundation. 
By  using  this  approximation  and  if  the  width  of  the  cemented  zone  is  small  compared  to  the 
grain  radius,  Dvorkin  et  al.  (1991)  were  able  to  reduce  the  problem  of  grain-cement 
deformation  (where  the  grains  are  deformable)  to  an  ordinary  integral  equation  for  the 
normal  stresses  at  the  cemented  interface.  The  stiffness  of  a  cemented  system  increases 
with  the  radius  of  the  cement  layer  and  with  its  stiffness.  The  former  factor  is  the  most 
important:  the  small  increase  of  the  cementation  content  results  in  significant  growth  of  a 
contact  zone  between  two  contacting  grains  and,  therefore,  dramatically  increases  the 
stiffness  of  a  two-grain  system  as  well  as  the  macroscopic  stiffness  of  a  cemented 
particulate  material  (even  if  the  cement  is  relatively  soft). 

These  theoretical  predictions  were  supported  by  the  experiments  of  Liu  et  al.  (1991) 
where  a  granular  material  was  represented  by  identical  glass  beads,  and  the  process  of 
cementation  was  simulated  by  freezing  the  capillary  water  accumulated  at  the  bead  contacts. 
The  measured  compressional  velocity  in  the  material  had  small  constant  values  before 
freezing  and  increased  sharply  at  the  freezing  point.  This  effect  is  apparently  due  to  the 
solidification  of  the  pendular  rings  which  form  ice  intergranular  cementation.  Qualitatively 
similar  to  these  are  our  measurements  of  velocities  in  epoxy-cemented  glass  beads 
(discussed  below).  The  important  message  of  the  experiments  is  that  the  velocity  increase 
with  cementation  is  very  large  between  0%  and  10%  saturation  and  small  between  10%  and 
50%  saturation  (Figure  1).  Therefore,  cementation  can  dramatically  increase  the  stiffness 
of  a  granular  material  only  if  it  is  placed  around  grain-to-grain  contacts. 

In  this  paper  we  concentrate  on  the  theoretical  modeling  of  effective  elastic  properties  of 
random  packing  of  identical  elastic  spheres  with  elastic  cement  at  their  contacts.  The 
effective  bulk  and  shear  moduli  of  a  packing  are  calculated  from  the  average  number  of 
contacts  per  sphere  and  porosity  of  the  aggregate,  and  from  the  normal  and  tangential 
stiffnesses  of  a  two-grain  combination  (Digby,  1981).  These  stiffnesses  are  found  by 
solving  the  problems  of  normal  and  tangential  deformation  of  two  elastie  spherical  grains 
cemented  at  their  contact.  The  approach  used  in  obtaining  this  solution  is  similar  to  that  of 
Dvorkin  et  al.  (1991):  a  thin  cement  layer  is  approximated  by  an  elastic  foundation,  and  the 
grain-cement  interaction  problems  are  reduced  to  linear  integral  equations. 
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Our  mathematical  approach  to  solving  these  elastic  contact  problems  is  similar  to  those 
used  in  many  contact-dominated  theories.  Among  them  are  solutions  that  describe  the 
effective  viscosity  of  concentrated  suspensions  (e.g.,  Frankel  and  Acrivos,  1967; 
Goddard,  1977;  Nunan  and  Keller,  1984);  adhesive-contact  interaction  (e.g.,  Johnson,  et 
al.,  1971);  and  thermal  and  electrical  conduction  through  a  granular  material  (e.g., 
Batchelor  and  O'Brien,  1976). 


Figure  1.  P-wave  and  S-wave  velocities  measured  in  epoxy-cemented  glass  beads  at  varying  cement 
saturation  at  hydrostatic  confining  pressure  30  MPa.  The  circles  indicate  the  measurement  points.  The 
experiments  have  been  conducted  using  the  ultrasonic  pulse  transmission  method.  The  experimental 
data  is  compared  to  our  upper  and  lower  bound  theoretical  predictions  —  solid  lines  obtained  under  two 
different  assumptions  of  cement  placement  around  the  grains. 

The  solution  reveals  a  peculiar  pattern  of  normal  and  shear  stresses  distribution  at  the 
cemented  grain  contacts:  the  stresses  are  maximum  at  the  center  of  the  contact  region  when 
the  cement  is  soft  relative  to  the  grain,  and  are  maximum  at  the  periphery  of  the  contact 
region  when  the  cement  is  stiff.  Stress  distribution  shape  gradually  varies  between  these 
two  extremes.  This  result  provides  an  important  insight  into  possible  modes  of 
intergranular  cement  failure:  relatively  stiff  cement  will  yield  at  the  periphery  of  a  cement 
layer.  Our  theoretical  solution  shows  that  it  is  mainly  the  amount  of  cement  that  influences 
the  effective  elastic  properties  of  cemented  granular  materials.  The  radius  of  the  cement 
layer  affects  the  stiffness  of  a  granular  assembly  much  more  strongly  than  the  stiffness  of 
the  cement  does. 

Our  theoretical  prediction  of  the  effective  properties  of  cemented  granular  materials  is 
applicable  to  a  random  packing  of  identical  spheres.  The  solutions  to  the  problems  of 
normal  and  shear  deformation  of  a  two-grain  combination  are  more  general  and  thus  can  be 
used  as  constitutive  contact  laws  in  rigorous  numerical  simulations.  Still,  caution  should 
be  exercised  when  substituting  these  contact  laws  into  mean-field  approximations,  such  as 
the  above-mentioned  Digby's  solution.  Such  approximations  may  become  invalid  in 
highly-disordered  or  contact-depleted  systems  (e.g.,  Goddard,  1990). 

Finally,  the  problems  of  normal  and  tangential  deformation  of  two  cemented  elastic 
cylinders  (the  plane  case)  are  considered  in  the  Appendix. 

Effective  Properties  of  Identical  Sphere  Packings 

The  effective  elastic  properties  of  a  random  packing  of  identical  spherical  particles  can 
be  expressed  through  its  porosity  (j),  coordination  number  C  (the  average  number  of 
contacts  per  sphere),  the  radius  of  a  particle  R,  and  the  normal  (S„)  and  tangential  (5^) 
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stiffnesses  of  a  two-sphere  combination  (Digby,  1982;  and  Winkler,  1983).  The  normal 
and  shear  stiffnesses  are  defined  as  the  ratios  of  a  corresponding  force  increment  to  the 
displacement  of  the  sphere  center  relative  to  the  contact  region  (Figure  2): 


35'  '  dt' 


Effective  bulk  ( ^^^)and  shear  ( moduli  are: 


_C^ 

l27tR  ” 


207tR  "  2  " 


(1) 


(2) 


Figure  2.  Normal  and  shear  deformation  of  a  two-grain  combination. 

Therefore,  the  effective  compressional  and  shear  velocities  can  be  found  as 


”  207tRp 


(S„+-S,),  V^  = 


C  3 

-(S„+-S,), 
20kRp  2  " 


where  p  is  the  density  of  the  grain  material. 

NORMAL  DEFORMATION  OF  TWO  CEMENTED  SPHERES 


Normal  Stiffness  and  Effective  Bulk  Modulus 

The  solution  to  the  problem  of  the  normal  deformation  of  two  cemented  grains 
(Dvorkin  et  al.,  1991)  is  based  on  the  assumption  that  a  contact  region  is  small  as  compared 
to  the  grain.  Therefore,  normal  displacements  v  of  the  surface  of  an  elastic  grain  due  to  a 
concentrated  force  P  can  be  expressed  in  the  (x,y,z)  coordinate  system  (Figure  3)  as  those 
of  the  surface  of  an  elastic  half-space  (Johnson,  1985); 


v(x,y)= 


(l-v)f^ 

27tG  r ' 


(3) 


where  G  and  v  are  the  shear  modulus  and  Poisson's  ratio  of  the  grain,  respectively. 

We  assume  that  the  cemented  region  is  symmetrical  with  respect  to  the  line  that 
connects  the  centers  of  the  two  contacting  spheres  (Figure  4a).  Therefore,  the  contact 
region  on  the  surface  of  the  sphere  is  a  circle  of  radius  a  (Figure  4b),  and  contact  stresses 
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and  displacements  are  axisymmetrical,  and  depend  on  the  radius  r  only.  Normal  stresses 
p{r)  are  related  to  normal  displacements  v(r)  as  (Timoshenko  and  Goodier,  1970) 


v(r)=-^^— ^  {d(p  J p(-\jr^+s^ -2rscos(p  )ds, 

7tG  i  i 


rcos<p+^Ja^-r^sin^(p 


(4) 


where  integration  is  conducted  inside  the  circle  Irk  a  (Figure  5). 


Figure  3.  Concentrated  forces  acting  upon  the  surface  of  an  elastic  half-space.  A  cemented  contact 
region  between  two  grains  is  assumed  to  be  small  as  compared  to  the  grain.  On  the  basis  of  this 
assumption,  displacements  of  the  surface  of  the  grain  are  related  to  stresses  through  the  formulas  for  an 
elastic  half-space. 

The  displacement  of  the  center  of  the  spherical  grain  relative  to  the  median  plane  of  the 
cement  layer  S  can  be  related  to  the  displacements  V’(r)  of  the  surface  of  the  cement  layer 
and  v(r)  of  the  surface  of  the  grain  (Figures  4b,  c,  and  d)  as 


5=v(r)-y(r).  (5) 

The  results  of  Dvorkin  et  al.  (1991)  show  that  a  thin  cement  layer  subject  to  normal  and 
shear  load  can  be  approximately  treated  as  an  elastic  foundation.  Therefore,  normal 
stresses  p  acting  upon  the  surface  of  the  grain  are  related  to  the  displacements  V  of  the 
surface  of  the  cement  layer  as 


p(r)=- 


2G,(l-V,)V(r) 
1-2  h(r)’ 


(6) 


where  and  are  the  shear  modulus  and  Poisson's  ratio  of  the  cement,  and  h(r)  is 
half-thickness  of  the  cement  layer.  For  the  case  under  consideration  we  have 


h(r)=R[£+U^n 

L  K 


(7) 


By  assuming  in  addition  that  shear  stresses  at  the  grain  surface  do  not  significantly 
influence  its  normal  deformation  (Johnson,  1985),  we  can  combine  equations  (4)  -  (7)  into 
the  following  integral  equation: 
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n  rcos(p+^ja^-r^sin^(p 

S+V{r)=-A]d(p  J 

0  0 


+s^-2rscos(p) 


1  rs 


ds. 


2G,  (l-v)(l-v,) 
TtG  1-2  V, 


The  constant  S  here  can  be  determined  from  the  resulting  compressional  force  F: 


u 

F= J  p{r)2Krdr=- 


4;rG,(l-vJ  fV(r)r 
1-2  V  1  h{r) 


dr. 


B 

4"^ 

Cement 

f 

B 

A 

a 

A 

J 

Figure  4.  A  cemented  contact  region  between  two  spherical  grains:  a)  The  cement  forms  an 
axisymmetrical  thin  layer,  b)  Normal  deformation  of  the  grain  and  the  cement  --  5  is  the  displacement 
of  the  center  of  the  sphere  relative  to  the  median  plane  of  the  cement  layer  AA  (this  plane  moves  into 
position  BB).  c)  Normal  axisymmetrical  displacement  of  the  cement  layer,  d)  Normal  axisymmetrical 
displacement  of  the  surface  of  the  grain. 


Figure  5.  The  region  of  contact  on  the  grain  surface  lrl<  a  in  the  plane  z  =  0. 

After  normalizing  the  above  equations,  and  using  (1),  we  arrive  at  the  following 
scheme  of  determining  normal  stiffness 
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1 .  Find  an  as  yet  unknown  function  H{t)  from  the  integral  equation 


tcos(p+^ja^- 


A^+H(t)=-Ajd(p  J 


'  ’’  H{^Jt^+s^-2tscos(p) 

£+^(r^ ^  -2t5cos  (p) 


ds. 


where  a  =  al  R,  and  A,,  is  an  arbitrarily  chosen  non-zero  constant. 

2.  Calculate  integral 


I  e+t^  /  2  “  ■ 

3.  Find  normal  stiffness  as 
47tRG,(l-v,)  k 

1-2  V,  Ao’ 

Formula  (2)  now  yields  the  following  final  expression  for  the  effective  bulk  modulus  of 
a  random  sphere  packing  of  porosity  (j)  and  coordination  number  C : 

1-2 V,  3(1+£)^  a/’ 

where  e  is  proportional  to  the  minimal  thickness  of  the  cement  layer  (at  r  =  0),  and  e  =  0 
for  the  case  where  grains  have  direct  point  contacts. 

The  above  integral  equations  allow  for  a  straightforward  numerical  solution  using  the 
quadrature  method  (e.g.,  Delves  and  Mohamed,  1985). 

Normal  Stress  Distribution  at  the  Contact 

We  used  the  above  equations  to  calculate  normal  stress  distribution  at  the  contact  of  two 
cemented  elastic  spheres.  In  the  following  example,  the  Poisson's  ratios  of  both  the  grain 
and  the  cement  materials  were  constant:  v  =  =  0.28,  the  densities  were  also  constant 

and  equal  2480  kg/m^.  The  minimal  separation  between  grains  was  zero  ( e  =  0),  the  ratio 
of  the  radius  of  the  cement  layer  to  the  grain  radius  was  0.33.  We  considered  six  cases 
where  the  rigidity  of  the  cement  material  relative  to  that  of  the  grains  varied  from  very  soft 
cement  ( '  ^pgrain  =  0-2)  to  relatively  stiff  cement  ( =  1.4),  where 
is  the  compressional  velocity  in  a  material. 

The  results  reveal  a  peculiar  picture  of  normal  stress  distribution  at  the  cemented  grain 
contacts:  the  stresses  are  maximum  at  the  center  of  the  contact  region  when  the  cement  is 
soft  relative  to  the  grain,  and  are  maximum  at  the  periphery  of  the  contact  region  when  the 
cement  is  stiff.  Stress  distribution  shape  gradually  varies  between  these  two  extremes  as 
the  cement's  stiffness  increases.  These  results  allow  for  a  simple  physical  interpretation: 
When  a  soft  cement  layer  is  confined  between  two  rigid  grains,  its  maximum  compression 
occurs  at  its  thinnest  part  -  near  the  point  of  direct  grain-to-grain  contact.  If  the  cement  is 
stiff,  the  problem  approaches  that  of  a  rigid  punch  (the  cement)  penetrating  an  elastic  half¬ 
space  (the  grain).  In  the  latter  case,  normal  stress  becomes  infinite  at  the  edge  of  the  punch 
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(Johnson,  1985)  which  explains  the  observed  maximum  of  the  normal  stresses  at  the 
periphery  of  the  contact  region.  One  practical  implication  of  the  observed  normal  stress 
distribution  is  that  a  relatively  stiff  intergranular  cement  layer  will  tend  to  yield  at  its  edges. 


0  0.5  1 

Coordinate  along  cement  layer 
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Figure  6.  Normal  stress  distribution  along  the  radii  of  cement  layers  for  varying  relative  rigidity  of  the 
cement  and  the  grain  materials  --  from  very  soft  to  very  stiff  cement.  Stresses  are  normalized  by  the 
average  stress.  The  horizontal  axis  is  the  normalized  distance  along  the  radius  of  the  cement  layer. 
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TANGENTIAL  DEFORMATION  OF  TWO  CEMENTED  SPHERES 


Tanjgential  Stiffness  of  a  Two-Grain  Combination 

Our  approach  to  solving  the  problem  of  tangential  deformation  of  two  cemented  grains 
with  the  relative  displacement  2t  in  the  x  direction  (the  coordinate  system  on  the  surface 
of  a  grain  is  shown  in  Figure  3)  is  similar  to  the  one  used  for  the  problem  of  iiormal 
deformation.  Again,  we  assume  that  the  contact  region  is  small  compared  to  the  grain  and, 
therefore,  the  formulas  for  an  elastic  half-space  can  be  used  to  relate  stresses  on  the  grain 
surface  to  its  displacements. 

A  concentrated  tangential  force  when  applied  to  an  elastic  half-space  (Figure  3), 
produces  the  following  displacements  of  the  surface  (Johnson,  1985): 


M=-%;-[(l+^)+a-2v)(l-^)], 

4;rG  r  r  r 


w 


2kG  r 


AkG 


(8) 


at  z  =  0 . 

It  follows  from  the  above  equations  that  the  problem  under  consideration  is  not 
axisymmetrical  as  in  the  case  of  normal  compression.  Indeed,  a  tangential  deformation  of  a 
two-grain  combination  in  the  x  direction  will  produce  displacements  in  the  cement  layer  not 
only  in  the  x  but  also  in  the  y  direction.  The  latter  are  symmetrical  with  respect  to  the  x 

axis  and  thus  are  zero  at  y  =  0  and,  generally,  non-zero  elsewhere. 

Similar  to  (4),  the  following  integral  equations  can  be  obtained  to  relate  tangential 
stresses  and  qy  distributed  in  the  circle  of  radius  a,  to  the  corresponding  displacements 
u  and  w: 


^  2%  L{r,6,(p) 

u{r,  0)=:—  jdcp  j  [q,(r^,  0,)(1-  vsin>)-^/r„  vsmq)cos(p]ds, 

2nG  Q  Q 

.  Ill  Lir,Q,(p) 

w{r,6)= - jd(p  j[qy(r^,e^)(l-vcos^(p)-q^(r^,d^)vsm(pcos(p]ds, 

2^^  0  0 

where 


(9) 


L(r,e,(p)=a[Jl-(-f  sin^(6+(p)  +-cos(e+(p)], 

V  a  a 

n — X - — — -  ^  ,  r sin 0-1- 5 sine? , 

f  =Jr  +s  - 2r5COs( 9+(p) ,  d. = arctan( - - - ), 

^  ^  rcosO-scos(p 

r=Jx^+y^,  0=arctan(— ). 

X 

The  domain  of  integration  in  the  z  =  0  plane  is  shown  in  Figure  7. 

Similar  to  the  above  case  of  normal  deformation,  the  following  coinpatibility  equations 
can  be  written  for  the  tangential  displacements  of  the  cement  and  the  grain: 
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t=u{r,d)-U(r,d),  w(r,6)=Wir,6), 

where  U  and  W  are  the  tangential  displacements  of  the  surface  of  the  cement  layer  in  the  x 
and  y  directions,  respectively,  and  T  is  the  tangential  displacement  of  the  center  of  the 
grain  along  the  x  direction. 


Figure  7.  The  domain  of  integration  used  in  equations  (9). 

Employing  in  addition  the  assumption  that  a  cement  layer,  when  deforming  in  shear, 
can  be  treated  as  an  elastic  foundation: 


we  can  transform  (9)  into  the  following  system  of  two  integral  equations: 


T+Uir,d)= 

— ^\d(p  f[  ^^-^^^^^(l-vsinV)-^.^^'’^"~^sin^cos^]tf.y, 

^  J  h(r)  h{r^) 


0  0 


W{r,e)= 


J  ^  J  h{r)  h(r?i 


0  0 


HO 


vsm(pcos(p]ds. 


For  the  case  of  two  spherical  grains  in  contact,  h  is  given  by  (7). 

In  these  equations,  t  can  be  determined  from  the  resulting  tangential  force  T: 


a  2k 

T=j  jq^ir,d)rdedr 

0  0 


-Gjj 


Uir,9)r 

hir) 


dddr. 


The  resulting  tangential  stiffness  can  be  found  from  (1). 


The  above  system  of  two  integral  equations  for  determining  the  tangential  stiffness  of  a 
two-grain  combination  can  be  simplified  and  reduced  to  one  integral  equation.  It  follows 
from  (8)  that  a  force  acting  in  the  x  (0  =  0)  direction  produces  tangential  displacements  u 
and  w  with  the  ratio 
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w(r,d)_  vsin0cos0 
u(r,6)  l-vsin^0 

This  ratio  monotonously  increases  with  the  increasing  v.  Its  maximum  value  for  v  =  0.4 
is  about  0.25.  Therefore,  a  tangential  force  acting  on  an  elastic  half-space  produces  lateral 
displacements  that  do  not  exceed  25%  of  the  displacements  in  the  parallel  direction.  By 
rewriting  the  first  integral  equation  in  (10)  as 


u(r,  6)=- 


G. 


iTt  L{r,d,(p) 

Jl 


27cG 


“ ( ''c ’ ^ (1  -  V sin^ ^) - vsin (pco^(p]ds. 


HO 


h{0 


we  conclude  that  the  contribution  of  these  lateral  displacements  is  no  more  that  25%  of  that 
of  the  displacements  in  the  parallel  direction.  Therefore,  the  total  error  of  neglecting 
displacement  w  does  not  exceed  6.25%. 

Based  on  this  approximation,  we  can  reduce  (10)  to  a  single  integral  equation 


'C+U{r)=-h,\dq>  j 


’’  U  ^ -Irscos  ^ )  (1-  V  sin^  (p) 


ds. 


The  constant  T  can  be  determined  from  the  resulting  tangential  force  T : 


T=  f  q^{r)2Krdr=-27tGA - — dr. 

0/?[£-hi(-^)^] 

L  K 

After  normalizing  the  above  equations  for  tangential  deformation,  we  arrive  at  the 
following  scheme  of  determining  tangential  stiffness  of  a  two-grain  combination: 


1.  Find  an  as  yet  unknown  function  H(t)  from  the  integral  equation 


tcos(p+^[^-t^sm^  (p 


\  +  H{t)=-0\d(p  J 


H(Vr^  +5^-2t^cos^)(l-  vsin^^) 
e-i-^(r^  4-s^  -2tscos(p) 


ds. 


where  a  =  a  /  R,  and  Aq  is  an  arbitrarily  chosen  non-zero  constant. 

2.  Calculate  integral 


I  e+t^  /  2  “  ■ 

3.  Find  tangential  stiffness  as 

S=-2kRG^—. 

A 
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Formula  (2)  now  yields  the  following  final  expression  for  the  effective  shear  modulus 
of  a  random  sphere  packing  of  porosity  ^  and  coordination  number  C : 


3C(l-<^)  k 
20(l+e)^ 


where  can  be  found  from  the  above  solution  for  the  normal  deformation  of  two 
cemented  grains. 


Shear  Stress  Distribution  at  the  Contact 

We  used  the  above  equations  to  calculate  shear  stress  distribution  at  the  contact  of  two 
cemented  elastic  spheres.  All  parameters  in  the  following  example  are  identical  to  those  in 
the  above  example  for  normal  stresses. 

The  results  (Figure  8)  are  similar  to  those  for  normal  stresses:  shear  stresses  are 
maximum  at  the  center  of  the  contact  region  when  the  cement  is  soft  relative  to  the  grain, 
and  are  mayimnm  at  the  periphery  of  the  contact  region  when  the  cement  is  stiff.  The  shape 
of  the  stress  distribution  curve  gradually  varies  between  these  two  extremes  with  the 
increasing  stiffness  of  the  cement. 

EXAMPLES  -  EFFECTIVE  PROPERTIES  OF  CEMENTED  SPHERICAL 
GRAINS 


Effect  of  the  Amount  of  Cement 

In  the  following  example  we  explore  the  influence  of  the  amount  of  ceinentation  on  the 
effective  properties  (compressional  and  shear  velocities)  of  a  random  packing  of  identical 
spheres.  The  porosity  of  the  packing  (without  cementation)  was  0.36,  the  coordination 
number  was  assumed  9.  The  properties  of  the  grain  material  were:  density  2480  kg/m^, 
bulk  modulus  49.9  GPa,  and  shear  modulus  26.2  GPa  -  parameters  typical  for  glass.  The 
density  of  the  cement  material  was  1 160  kg/m^.  We  explored  four  cases  of  the  increasing 
stiffness  of  the  cement:  a.  =  3.4  GPa,  =  1.0  GPa;  b.  =  6.8  GPa,  =  2.0 
GPa;  c.  =  13.6  GPa,  G^  =  4.0  GPa;  and  d.  =  27.2  GPa,  =  8.0  GPa.  The 
second  case  here  corresponds  to  typical  epoxy  parameters.  The  ratio  of  the  cement  layer 
radius  a  to  the  grain  radius  varied  from  0.01  to  0.5. 

Plots  in  Figure  9  show  that  both  compressional  and  shear  velocities  dramatically 
increase  with  the  increasing  radius  of  the  cement  layer.  The  influence  of  this  parameter  is 
clearly  much  stronger  than  that  of  the  cement  stiffness:  even  a  small  increase  in  the  amount 
of  relatively  soft  cement  gives  a  significant  rise  to  both  and  V^.  These  theoretical 
predictions  are  well  supported  by  the  earlier  experiments  of  Liu  et  al.  (1991)  which  are 
discussed  in  the  Introduction,  and  by  our  experiments  on  epoxy-cemented  glass  beads. 

The  effective  Poisson’s  ratio  of  the  cemented  sphere  packings  increases  with  the 
increasing  amount  of  cementation  (Figure  10).  The  effect  of  cement  stiffness  is  opposite  to 
that  in  the  case  of  compressional  and  shear  velocities:  the  Poisson's  ratio  decreases  with 
the  increasing  stiffness. 


Effect  of  Cement  Stiffness 

We  explored  the  effect  of  cement  stiffness  on  compressional  and  shear  velocities 
(Figure  11)  in  the  above-described  random  packings  of  identical  spheres  for  constant 
amount  of  cement  {a/  R  =  0.1).  The  parameters  of  the  grain  material  were  kept  constant 
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with  bulk  and  shear  moduli  49.9  and  26.2  GPa,  respectively,  whereas  the  stiffness  of  the 
cement  varied  from  =3.4  GPa  and  =  1.0  GPa  to  =  68  GPa  and  =  20  GPa. 
This  dramatic  increase  of  cement  stiffness  results  in  only  about  a  15%  increase  in 
compressional  and  a  20%  increase  in  shear  velocity  -  an  effect  much  smaller  than  that  of 
increasing  cement  content. 


Coordinate  along  cement  layer  Coordinate  along  cement  layer 

(normalized)  (normalized) 


0  0.5  1  0  0.5  1 

Coordinate  along  cement  layer  Coordinate  along  cement  layer 

(normalized)  (normalized) 


0  0.5  1 

Coordinate  along  cement  layer 
(normalized) 


0  0.5  1 

Coordinate  along  cement  layer 
(normalized) 


Figure  8.  Shear  stress  distribution  along  the  radii  of  cement  layers  for  varying  relative  rigidity  of 
the  cement  and  the  grain  materials  -  from  very  soft  to  very  stiff  cement.  Stresses  are  normalized 
by  the  average  stress.  The  horizontal  axis  is  the  normalized  distance  along  the  radius  of  the  cement 
layer. 


17 


Figure  9.  Compressional  and  shear  velocities  in  random  packings  of  identical  cemented  spherical 
particles  as  functions  of  the  amount  of  cementation.  The  parameters  of  the  grain  material  are  constant; 
the  parameters  of  the  cement  vary  from  soft  ("a")  to  stiff  ("d")  cement. 

EXPERIMENTAL  RESULTS 

In  order  to  check  the  validity  of  our  theoretical  model,  we  have  conducted  ultrasonic 
(0.42  MHz)  pulse  transmission  experiments  on  epoxy-cemented  glass  beads.  We  used 
spherical  glass  beads  of  0.42  -  0.5  mm  diameter.  The  material  contains  no  more  than  two 
percent  of  irregularly-shaped  grains,  and  is  reasonably  free  of  sharp  angular  particles,  or 
particles  showing  milkiness,  surface  scoring,  and  foreign  matter.  The  porosity  of  the 
uncemented  glass  bead  sample  was  about  36  percent  —  close  to  the  porosity  of  a  dense 
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random  pack  of  identical  spheres.  The  cemented  sample  was  prepared  by  mixing  glass 
beads  with  a  given  volume  of  epoxy.  The  mixture  was  placed  into  a  holder  where  the 
epoxy  solidified. 

The  measured  compressional-  and  shear-wave  velocity  is  compared  to  the  theoretical 
predictions  in  Figure  1.  The  upper  theoretical  curve  (upper  bound  estimate)  corresponds  to 
the  case  where  all  cement  is  deposited  precisely  at  grain  contacts,  the  coordination  number 
here  was  chosen  as  9.  The  lower  theoretical  curve  (lower  bound  estimate)  was  calculated 
for  the  case  where  a  contact  cement  region  is  formed  as  the  intersection  of  cement  shells 
that  are  uniformly  deposited  on  grain  surfaces.  For  a  random  pack  of  identical  spheres,  the 
average  number  of  contacts  per  grain  is  between  8  and  9.  Therefore,  we  choose  8  as  the 
coordination  number  for  this  lower  bound  estimate.  The  measurements  are  within  the 
predicted  bounds. 


Figure  10.  Poisson’s  ratio  of  random  packings  of  identical  cemented  spherical  particles  as  a  function  of 
the  amount  of  cementation.  The  parameters  of  the  grain  material  are  constant;  the  parameters  of  the 
cement  vary  from  soft  ("a")  to  stiff  ("d")  cement. 

CEMENT  DEPOSITION,  POROSITY  AND  EFFECTIVE  PROPERTIES 

The  following  example  illustrates  how  two  possible  schemes  of  cement  deposition  may 
affect  the  porosity  and  the  elastic  moduli  of  a  granular  material.  Both  grain  and  cement 
materials  were  identical  (glass),  which  is  relevant  to  the  process  of  quartz  cementation  of 
unconsolidated  well-sorted  sands.  The  process  of  cementation  starts  at  the  critical  porosity 
point  (about  36%  for  a  random  identical  sphere  packing)  at  which  both  bulk  and  shear 
moduli  are  zero  in  an  unconsolidated  granular  material.  The  first  model  (Arrangement  1)  is 
where  aU  added  cement  is  accumulated  at  the  intergranular  contacts.  In  the  second  model, 
(Arrangement  2)  the  cement  evenly  covers  the  grains  and  the  contact  cement  layer  is  formed 
as  an  intersection  of  these  cement  shells.  It  is  clear  that  in  the  first  arrangement,  the  size  of 
a  contact  cement  layer  is  larger  (for  given  porosity)  than  in  the  second  arrangement.  This 
effect  results  in  a  steeper  increase  of  the  effective  moduli  with  decreasing  porosity  in  the 
first  model  as  compared  to  the  second  one  (Figure  12).  Both  bulk  and  shear  moduli  steeply 
increase  at  the  very  beginning  of  cement  deposition  because  even  a  small  addition  of  cement 
at  intergranular  contact  points  leads  to  the  fast  growth  of  the  radius  of  the  cement  layer.  It 
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is  important  to  note  that  our  theoretical  model  is  applicable  only  for  small  amounts  of 
cement  at  intergranular  contacts. 
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Figure  1 1.  Compressional  and  shear  velocities  in  a  random  packing  of  identical  cemented  spherical 
particles  as  a  function  of  cement  stiffness.  The  ratio  of  the  radius  of  the  cement  layer  to  grain  radius  is 
constant  and  equals  0. 1 . 

CONCLUSIONS 

A  theoretical  solution  is  given  to  the  problems  of  normal  and  tangential  deformation  of 
two  cemented  elastic  grains.  The  elastic  moduli  of  grains  may  differ  from  those  of  cement. 
This  solution  is  used  to  describe  the  effective  elastic  properties  of  a  cemented  granular 
material  that  is  modeled  as  a  random  packing  of  identical  spheres.  The  effective  bulk  and 
shear  moduli  of  the  packing  are  calculated  from  porosity  and  coordination  number,  and 
from  the  normal  and  tangential  stiffnesses  of  a  two-grain  combination.  A  thin  elastic 
cement  layer  is  approximated  by  an  elastic  foundation,  and  the  grain-cement  interaction 
problems  are  reduced  to  linear  integral  equations.  The  theoretical  solution  shows  a  peculiar 
pattern  of  normal  and  shear  stresses  distribution  at  the  cemented  grain  contacts:  the  stresses 
are  maximum  at  the  center  of  the  contact  region  when  the  cement  is  soft  relative  to  the 
grain,  and  are  maximum  at  the  periphery  of  the  contact  region  when  the  cement  is  stiff. 
Stress  distribution  shape  gradually  varies  between  these  two  extremes  as  the  cement's 
stiffness  increases.  The  solution  reveals  that  the  radius  of  the  cement  layer  has  a  much 
stronger  effect  on  the  stiffness  of  a  granular  assembly  than  does  the  stiffness  of  the  cement. 
This  result  provides  an  important  insight  into  possible  modes  of  intergranular  cement 
failure:  relatively  stiff  cement  will  fail  at  the  periphery  of  a  cement  layer.  The  results  of  our 
theoretical  solution  to  the  problem  of  deformation  of  a  two-grain  combination  provide 
constitutive  contact  laws  that  can  be  used  in  rigorous  numerical  simulations.  It  is  important 
to  keep  in  mind  that  in  this  paper  we  present  the  effective  elastic  constants  of  a  composite 
material  that  have  been  obtained  from  the  static  analysis.  These  constants  cannot  be  used  in 
predicting  wave  velocities  if  frequencies  are  so  high  that  wavelengths  are  comparable  to  the 
individual  grain  size.  However,  static  elastic  constants  can  be  used  with  confidence  if  the 
wavelength  is  significantly  larger  that  the  individual  grain  size.  In  our  experiments,  the 
wavelengths  were  about  twelve  times  the  average  grain  diameter. 
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Porosity 


Figure  12.  Bulk  and  shear  moduli  of  a  random  packing  of  identical  glass-cemented  glass  beads  for  two 
models  of  cement  deposition.  The  moduli  are  normalized  by  those  of  glass.  The  porosity  of  the 
aggregate  reduces  with  the  increasing  amount  of  the  cement. 
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APPENDIX:  DEFORMATION  OF  TWO  CEMENTED  CYLINDERS  -- 
PLANE  CASE 

The  problem  of  normal  deformation  of  two  cemented  elastic  cylinders  was  considered 
by  Dvorkin  et  al.  (1991).  The  assumptions  in  this  solution  are  similar  to  those  in  the  above 
case  of  two  spherical  grains:  a  thin  elastic  cement  layer  between  two  elastic  grains  is 
approximated  by  an  elastic  foundation,  and  a  contact  region  is  assumed  to  be  small  as 
compared  to  a  grain,  so  that  stress  distribution  in  the  grain  is  close  to  one  in  a  half-plane 
with  identical  surface  deformation. 

In  the  plane  case,  normal  displacements  of  the  grain  surface  v(x)  are  related  to  normal 
stress  p{x)  as  (Johnson,  1985) 


v(x)=- 


1-v 

kG 


a 

J  p(5)ln|x— 5|<is+const, 

~a 


(Al) 


where  x  is  a  coordinate  along  the  contact  interface  (consider  Figure  4  with  x  instead  of  r). 
By  using  equations  (5)  and  (6)  with  x  instead  of  r,  we  arrive  at  the  following  integral 

equation  for  V(x) : 
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V(x)=A  J^^^ln|x-5|J5+const,  (A2) 

-a  ^ 

where  V  is  the  normal  displacement  of  the  surface  of  the  cement  layer  (Figure  4),  and  h  is 
given  by  (7).  The  constant  in  the  right-hand  side  of  (A2)  can  be  found  using  the  condition 
of  a  given  integral  compressional  force  F  per  unit  length  of  an  elastic  cylinder: 


2GAI-V.)  }V{s) 

(1-2 V,)  lh(s) 


ds. 


(A3) 


Equation  (A2)  with  condition  (A3)  can  be  easily  solved  by  the  quadrature  method.  The 
solution  results  in  normal  stress  distributions  similar  to  those  given  in  Figure  6  (Dvorkin  et 
al.,  1991). 

The  problem  of  tangential  deformation  of  two  cemented  cylinders  is  identical  to  that  of 
normal  deformation.  Indeed,  equation  (A2)  will  hold  if  we  substitute  tangential 
displacement  u  instead  of  v  and  shear  traction  q  instead  of  p  (Johnson,  1985).  The  final 
set  of  equations  for  the  tangential  displacement  of  the  surface  of  the  cement  layer  U  is: 


U f  ^ -—In [x-5|<i.y+ const, 
nG  iHs)  '  ' 


where  T  is  an  integral  tangential  force  per  unit  length  of  an  elastic  cylinder. 
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2.2.  CONTACT  LAWS  FOR  CEMENTED  GRAINS: 
IMPLICATIONS  FOR  GRAIN  AND  CEMENT  FAILURE 


ABSTRACT 

Analytical  solutions  are  presented  to  predict  the  intergranular  contact  load  transfer  in 
cemented  granular  media  where  both  grain  material  and  cement  are  elastic.  The  grains  can 
be  separated,  have  a  direct  point  contact,  or  be  compacted  prior  to  cement  deposition.  For 
all  these  cases  contact  stress  distributions  are  obtained  for  normal,  tangential,  and  torsional 
deformation  of  two  cemented  deformable  grains.  An  important  result  is  that  intergranular 
cement,  even  if  very  soft,  is  load-bearing.  Thus  cementation  reduces  contact  stress 
concentration  (as  compared  to  direct  Hertzian  interaction).  Contact  stresses  are  maximum 
near  the  center  of  the  contact  region  when  the  cement  is  soft  relative  to  the  grains,  and  are 
maximum  at  the  periphery  of  the  contact  region  when  the  cement  is  stiff.  These  results 
allow  us  to  predict  the  following  modes  of  static  and  dynamic  failure  of  the  grains  and 
intergranular  bonds  in  a  particulate  material:  (1)  Uncemented  grains  will  tend  to  shatter 
whereas  cemented  grains  will  stay  intact,  and  the  cement  will  fail.  (This  conclusion  is 
supported  by  hydrostatic  loading  experiments  where  intensive  crushing  of  uncemented 
glass  beads  was  observed  at  about  50  MPa,  whereas  grains  cemented  at  their  contacts  with 
small  amounts  of  epoxy  stayed  intact.)  (2)  Where  intergranular  cementation  is  present, 
grain  failure  may  still  be  expected  if  the  cement  is  strong  and  stiff.  In  this  case,  grain 
damage  will  be  initiated  at  the  periphery  of  the  cement  layer.  (3)  Yielding  of  a  cement 
material  that  is  soft  (as  compared  to  the  grain  material)  will  initiate  at  the  center  of  the 
contact  region,  whereas  stiff  cement  will  yield  at  the  periphery. 

INTRODUCTION 

The  study  of  dynamic  and  static  load  transmission  in  granular  media  is  important  in 
different  branches  of  engineering  and  geophysics.  Artificial  and  natural  granular  materials 
behave  as  good  shock  attenuators  and  as  such  are  used  to  isolate  shock-sensitive 
instruments.  Such  materials  can  be  used  to  protect  underground  facilities  from  explosion- 
associated  damage.  Many  rocks  can  be  treated  as  granular  media.  Understanding  and 
quantifying  shock  wave  propagation  in  rocks  helps  in  detecting  and  locating  nuclear 
explosions. 

The  mechanical  static  and  dynamic  characteristics  of  granular  materials  can  be 
significantly  affected  by  intergranular  cementation.  One  way  to  understand  and  quantify 
this  dependence  is  by  studying  grain-to-grain  microscale  interaction  because  in  a  dry 
granular  material  load  is  transferred  mainly  through  contact  mechanisms  between 
neighboring  particles  (e.g.,  Oda  et  al.,  1982;  Shukla  et  al.,  1988;  Sadd  et  al.,  1989). 

Recent  advancements  in  studying  the  effect  of  intergranular  cementation  on  the  elastic 
and  inelastic  behavior  of  particulate  materials  include  a  numerical  and  experimental 
investigation  by  Bruno  and  Nelson  (1991),  an  experimental  study  by  Bemabe  et  al.  (1992), 
and  numerical  (distinct  element  method)  works  by  Trent  (1989),  and  Trent  and  Margolin 
(1992).  One  important  detail  that  has  been  missing  in  all  these  cementation-related  works  is 
a  theoretical  description  of  interparticle  contact  laws.  Dvorkin  et  al.  (1991  and  1994)  have 
examined  normal  and  shear  stress  transmission  between  deformable  particles  through 
deformable  interparticle  cement  bonds.  Two  cemented  particles  could  be  separated  or  have 
a  point  direct  contact  (Figure  1).  Three  main  assumptions  were:  (1)  a  thin  cement  layer 
subject  to  normal  and  shear  load  can  be  approximately  treated  as  an  elastic  foundation,  (2) 
the  width  of  the  cemented  zone  is  small  as  compared  to  the  grain  radius  so  that  a  grain  can 
be  treated  as  an  elastic  half-space,  and  (3)  there  is  no  slip  between  grains  and  at 
grain/cement  interfaces.  The  first  assumption  has  been  justified  by  an  approximate 
analytical  solution  for  the  normal  and  shear  deformation  of  a  thin  smooth  cement  layer 
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(Dvorkin  et  al.,  1991).  Based  on  these  assumptions  the  problem  of  grain-cement 
deformation  has  been  reduced  to  an  ordinary  integral  equation  for  the  normal  and  shear 
stresses  at  the  cemented  interface. 


Figure  1.  Left:  Separated  grains,  grains  with  a  point  contact,  and  precompacted  grains.  Right: 

Normal,  shear,  and  torsional  deformation  of  two  cemented  grains. 

The  solution  has  revealed  a  peculiar  pattern  of  normal  and  shear  stresses  distribution  at 
the  cemented  grain  contacts:  the  stresses  are  maximum  at  the  center  of  the  contact  region 
when  the  cement  is  soft  relative  to  the  grain,  and  are  maximum  at  the  periphery  of  the 
contact  region  when  the  cement  is  stiff  (Figure  2).  The  physical  explanation  of  this  pattern 
is  as  follows:  When  a  soft  cement  layer  is  confined  between  two  rigid  grains,  strain  and 
stress  in  it  are  maximum  at  its  thinnest  part  —  near  the  center  of  contact.  When  the  cement 
is  stiff  relative  to  the  grain  material,  its  action  on  the  grain  is  close  to  that  of  a  rigid  punch 
penetrating  an  elastic  half-space  (e.g.,  Johnson,  1992).  In  this  case  stress  concentration  is 
expected  at  the  periphery  of  the  punch. 

The  macroscopic  stiffness  of  a  cemented  system  increases  with  the  increasing  radius  of 
the  cement  layer  and  with  its  increasing  stiffness.  The  former  factor  is  the  most  important: 
the  small  increase  of  the  cementation  content  results  in  significant  growth  of  a  contact  zone 
between  two  contacting  grains  and,  because  the  cement  is  load-bearing,  dramatically 
increases  the  stiffness  of  a  two-grain  system  and  thus  the  macroscopic  stiffness  of  a 
cemented  particulate  material.  This  conclusion  holds  even  if  the  cement  is  relatively  soft. 

In  this  paper  we  concentrate  on  the  description  of  normal,  shear,  and  torsional  stress 
transmission  between  two  elastic  grains  that  were  initially  precompacted  by  a  normal  force 
to  develop  a  finite  (Hertzian)  direct  contact  area.  Afterwards  cement  was  deposited  around 
the  initial  contact  zone  (Figure  1).  Our  analytical  solutions  have  been  obtained  under  the 
same  assumptions  as  in  the  uncompacted  case.  The  results  show  that  in  soft  cement, 
normal  and  shear  stresses  developed  due  to  the  continuing  deformation  of  two  grains  are 
maximum  at  the  periphery  of  the  initial  contact  zone  between  two  grains.  These  stresses 
reduce  towards  the  periphery  of  the  cement  layer.  However,  in  stiff  cement,  contact 
stresses  are  maximum  at  the  periphery  of  the  cement  layer. 

Below,  we  give  solutions  to  the  problems  of  normal  and  shear  deformation  of  two 
precompacted  cemented  spherical  grains.  Similar  solutions  for  two-dimensional  grains 
(cylinders)  are  given  in  Appendix  1.  Torsional  deformation  of  two  spherical  cemented 
grains  (both  uncompacted  and  precompacted  cases)  is  explored  in  Appendix  2.  The 
theoretical  results  are  compared  to  our  experimental  data. 
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We  also  discuss  the  implications  of  the  obtained  solutions  for  the  failure  of  cemented 
granular  materials. 
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Figure  2.  Normal  (A  and  B)  and  shear  (C  and  D)  stresses  in  the  cement  along  the  radii  of  the  cement 
layers  between  two  spherical  grains  with  a  point  contact.  The  Poisson’s  ratios  of  both  grain  and 
cement  materials  are  constant  and  equal  0.28.  The  stiffness  of  the  cement  varies  relative  to  that  of  the 
grains:  from  very  soft  cement  (A  and  C)  to  very  stiff  cement  (B  and  D).  is  compressional  wave 

velocity  in  the  material.  The  ratio  of  the  radius  of  the  cement  layer  to  the  grain  radius  is  0.33. 
Stresses  are  normalized  by  the  average  stress.  The  radial  coordinate  r  (Figure  1)  is  normalized  by  the 
radius  of  the  cement  layer. 

NORMAL  AND  SHEAR  DEFORMATION  OF  TWO  PRECOMPACTED 
GRAINS 


Precompacted  Grains 

Consider  two  identical  elastic  spherical  grains  of  radius  R  that  are  normally  pressed 
together  by  force  P^to  form  a  circular  contact  area  of  radius  a .  Then 


a=> 


3Ph^(1-v) 

'  8G  ’ 


(1) 


where  G  and  v  are  the  grain  shear  modulus  and  Poisson's  ratio,  respectively  (e.g., 
Johnson,  1992).  This  intergranular  force  P^^  can  be  easily  calculated  from  the  hydrostatic 
pressure  acting  at  a  random  pack  of  identical  spherical  grains: 
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(2) 


AkR^P^ 


where  C  is  the  average  number  of  contacts  per  grain  (about  9)  and  (j)  is  the  porosity  of  the 
pack  (about  36%). 

From  Hertz's  solution,  the  initial  distribution  of  normal  stresses  on  the  surface  of 
the  precompacted  grains  is: 


AG 


7tR{l-v) 
p^=0,r>  a. 


-r^ ,  r<a'. 


(3) 


The  corresponding  initial  displacements  of  the  grain  surface  are: 


a 


Vn(r)  =  —  -—^  0<r<a; 


Vh('')= 


2R 


tcR 


(2a^  - ) arcsin(— )  +  arJl-—. 


,  a<  r<b. 


(4) 


After  the  initial  compaction,  elastic  cement  is  uniformly  added  around  the  direct  grain 
contact  zone  to  increase  the  radius  of  the  contact  area  to  b  (Figure  3).  It  follows  from 
equations  (4)  that  the  half-thickness  of  the  undeformed  cement  layer  between  two 
precompacted  spherical  grains  is: 


^(0= 


a 


kR 


a 


J 


arctan 


a<r<b. 


(5) 


Figure  3.  Left:  Contact  area  between  two  precompacted  cemented  grains  with  cement  added  after 
compaction.  Right:  Cement  layer,  its  initial  thickness  and  its  displacements. 

Normal  Deformation 


Consider  now  additional  normal  loading  of  the  grains  with  force  AP  that  produces 
additional  displacement  8  of  each  grain's  center  towards  their  contact  zone.  The  resulting 
total  normal  displacements  v(r)  of  the  grain  surface  within  the  area  of  the  initial  direct 
contact  are: 
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(6) 


v(r)  =  6-  0  <  r  <  a. 

^  2R 

Within  the  cemented  zone,  a<r<b,  the  total  normal  displacements  of  the  grain 
surface  v(r)  are  related  to  the  normal  displacements  of  the  cement  as: 

V(r)  =  v(r)-Vu(r)-5  +  ^,  (?) 

K 

where  V{r)  is  the  normal  displacement  of  the  cement  surface  (at  the  cement-grain  interface) 
relative  to  the  median  plane  of  the  cement  layer  (Figure  3). 

By  treating  the  cement  layer  as  an  elastic  foundation,  we  have  the  following  expression 

for  normal  stresses  p{r)  at  the  cement-grain  interface: 


p{r)  =  - 


2G,(l-vJ  V(r) 
1-2V,  h„ir) 


,a<r<b, 


(8) 


where  and  are  the  cement's  shear  modulus  and  Poisson's  ratio,  respectively. 

Finally,  treating  the  grain  as  an  elastic  half-space,  we  have  the  following  expression 
that  relates  normal  stresses  p(r)  to  normal  displacements  v(r)  of  the  grain's  surface 
(Timoshenko  and  Goodier,  1970): 


rcos^p+Vfe^-r^sin^f _ 

v{r)-- — —\d(p  f  p(-\jr^+s^-2rscos(p)ds,  (9) 

i  i 

where  the  integration  is  inside  the  circle  |r|  <  b  (Figure  4). 

Now  we  combine  equations  (6),  (7),  (8),  and  (9)  to  obtain  the  following  integral 
equation  for  normal  stress  p(r): 


v(r)=^^  jd^  J p{^jr^+s^-2rscos(p)ds  = 

r 


rcosip+-\Jb^-r^  sin^  <p 


kG 


S - ,  0<r<a; 

2R 


VhG)  +  d- 


a 


l-2v 


R  2G(1-vJ 


p(r)h„{r),  a<r<b. 


(10) 


Let  us  present  the  normal  stress  p{r)  as  the  sum  of  the  initial  Hertzian  stress  py^{r)  and 
the  additional  normal  stress  /(r) : 


P(0  =  Ph(0 +  /('')• 

It  follows  then  from  equations  (4),  (9),  (10),  and  (11)  that 


(11) 
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2R 
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(12) 


This  last  equation  can  be  easily  solved  using  the  numerical  quadrature  method  (e.g.,  Delves 
and  Mohamed,  1985). 

The  constant  5  can  be  found  from  the  resulting  force  AP : 


h 

AP  =  J  f{r)2Krdr. 

0 


(13) 


Figure  4.  Domain  of  integration  in  equation  (9). 

Shear  Deformation 

Our  treatment  of  the  shear  deformation  of  two  cemented  grains  is  analogous  to  that  of 
the  normal  deformation;  we  assume  that  after  cement  deposition  a  tangential  force  Q  is 
exerted  on  the  system  to  produce  the  tangential  displacement  T  of  each  grain's  center 
relative  to  the  contact  zone.  The  resulting  tangential  displacements  u{r)  of  the  grain  surface 
within  the  area  of  the  initial  direct  contact  are: 

u{r)  =  t,0<r<a.  (14) 

Within  the  cemented  zone,  a  <  r  <  Z? ,  the  tangential  displacements  of  the  grain  surface 
u{r)  are  related  to  the  tangential  displacements  of  the  cement  as: 

U(r)  =  u(r)-T,  (15) 

where  U(r)  is  the  displacement  of  the  cement  surface  (at  the  cement-grain  interface)  relative 
to  the  median  plane  of  the  cement  layer. 

By  treating  the  cement  layer  as  an  elastic  foundation  in  shear,  we  have  the  following 
expression  for  tangential  stresses  qir)  at  the  cement-grain  interface: 

g(r)  =  -G,-^,a<rSb.  (16) 

hH(r) 
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Now  treating  the  grain  as  an  elastic  half-space,  we  have  the  following  approximate 
expression  that  relates  shear  stresses  q{r)  to  tangential  displacements  u{r)  of  the  grain's 
surface  (Dvorkin  et  al.,  1994): 


-I- 

j q(^jr^+s^ -Irscoscp  )(1  -  vsin^  (p)ds, 

0 


(17) 


where  the  integration  is  inside  the  circle  |r|  <  b  (Figure  4). 

Finally,  we  combine  equations  (14),  (15),  (16),  and  (17)  to  obtain  the  following 

integral  equation  for  shear  stress  q{r) : 


rcos<D+-\/fc^  -r^  sin^ _ 

jqi-^r^  +■5'^ -2rscos(p  )(1  -  vsin^  (p)ds  = 
0 


T,  0<r<a; 

t--^q{r)hf^(r),  a<r<b. 


(18) 


Again,  this  integral  equation  can  be  easily  solved  using  the  numerical  quadrature  method. 
The  constant  T  can  be  found  from  the  resulting  force  Q  : 


b 

Q  =  j  q{r)27trdr. 
0 


Contact  Stresses 


(19) 


Normal  stresses  in  the  contact  zone  of  two  precompacted  grains  are  given  in  Figure  5. 

The  ratio  of  the  cement  radius  b  to  the  grain  radius  R  was  0.7,  the  radius  a  of  the 
initial  zone  of  direct  grain  contact  was  0.105  of  R  (0.15  of  b).  The  elastic  moduli  of  the 
grains  were  constant  and  equal  to  those  of  glass  (bulk  modulus  49.9  GPa  and  shear 
modulus  26.2  GPa).  The  bulk  and  shear  moduli  of  the  cement  varied  from  1.7  GPa  and 
0.5  GPa,  respectively,  to  170.0  GPa  and  50.0  GPa.  In  the  second  example  (6.8  GPa  and 
2.0  GPa)  the  elastic  moduli  of  the  cement  are  identical  to  those  of  epoxy. 

When  the  cement  is  relatively  soft,  we  observe  the  concentration  of  the  additional 
normal  contact  stresses  at  the  periphery  of  the  direct  grain  contact  zone  (r  =  a).  The 
physical  meaning  of  this  effect  is:  In  the  soft  cement  layer  stresses  are  maximum  at  its 
thinnest  part  -  near  its  inner  boundary  (r  ^  a,  r>  a).  In  the  direct  grain  contact  region  we 

have  the  condition  of  a  constant  additional  displacement  d  —  analogous  to  that  on  the  face 
of  a  rigid  punch  penetrating  an  elastic  half-space.  This  is  the  reason  for  the  additional 
contact  stress  concentration  on  the  direct  contact  side  ( r  a,  r  <  a)  of  the  contact  region 

0<r<b. 

As  the  cement  becomes  stiffer,  its  action  on  the  grain  approaches  that  of  a  rigid  punch 
of  radius  b.  Therefore,  we  observe  stress  concentration  on  the  periphery  of  the  cement 
layer  at  r  =  ^? . 

Shear  stresses  for  the  six  above  cases  are  given  in  Figure  6.  The  patterns  of  stress 
distribution  and  their  physical  meaning  are  similar  to  those  observed  for  normal  stresses.  It 
is  important  to  notice  that  in  this  case  we  have  used  the  no-slip  condition  in  the  direct  grain- 
to-grain  contact  zone.  Clearly,  this  condition  is  not  valid  if  shear  stresses  are  high  and  slip 
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does  occur.  However,  we  justify  this  approximation  by  the  fact  that  for  small  tangential 
forces  exerted  on  contacting  bodies,  stress  distributions  that  are  produced  from  the  no-slip 
solutions  are  fairly  close  to  those  produced  from  the  partial-slip  solutions  (see  examples  in 
Johnson,  1992). 


Figure  5.  Additional  normal  stresses  /(r)  along  the  radii  of  cement  layers  for  varying  relative  rigidity 
of  the  cement  and  grain  materials.  Stresses  are  normalized  by  the  average  additional  stress.  The 
horizontal  axis  is  the  normalized  distance  along  the  radius  of  the  cement  layer. 
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Figure  6.  Shear  stresses  qir)  along  the  radii  of  cement  layers  for  varying  relative  rigidity  of  the 
cement  and  grain  materials.  Stresses  are  normalized  by  the  average  shear  stress.  The  horizontal  axis  is 
the  normalized  distance  along  the  radius  of  the  cement  layer. 

EXPERIMENT 

Our  experiment  has  been  conducted  on  hydrostatically  precompacted  glass  beads 
cemented  with  epoxy.  The  cemented  samples  were  prepared  by  mixing  glass  beads  (0.2  - 
0.3  mm  in  diameter)  with  a  given  volume  of  epoxy  (10,  25,  50,  and  100%  of  the  pore 
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space  volume).  They  were  compacted  at  20  MPa  prior  to  epoxy  solidification. 
Compressional  and  shear  (VJ  wave  velocities  have  been  measured  at  1  MHz 

frequency  using  the  pulse  transmission  technique. 

The  above-derived  contact  laws  for  the  normal  and  shear  deformation  of  two 
precompacted  cemented  grains  have  been  used  to  calculate  the  effective  bulk  and  shear 
moduli  of  a  random  packing  of  identical  elastic  spheres.  These  theoretical  predictions  are 
compared  to  the  experimental  results  in  Figure  7.  The  agreement  between  the  theory  and 
the  experiments  is  within  15%  accuracy. 


Figure  7.  Experimental  and  theoretical  results  on  precompacted  epoxy-cemented  glass  beads.  Circles 
show  experimental  points,  solid  lines  are  theoretical  predictions. 

IMPLICATIONS  FOR  GRAIN  AND  CEMENT  FAILURE 

One  important  result  following  from  the  above-derived  grain-cement-grain  contact  laws 
is  that  intergranular  cement,  even  if  very  soft,  is  load-bearing.  Thus  cementation  reduces 
contact  stress  concentration  (as  compared  to  direct  Hertzian  interaction).  As  a  result,  even 
small  amounts  of  relatively  soft  cement  may  prevent  grain  failure:  uncemented  grains  will 
tend  to  shatter  whereas  cemented  grains  will  stay  intact,  and  the  cement  will  fail.  This 
conclusion  has  been  supported  by  hydrostatic  loading  experiments  where  intensive 
crushing  of  uncemented  glass  beads  was  observed  at  about  50  MPa,  whereas  grains 
cemented  at  their  contacts  with  small  amounts  of  epoxy  (10%  of  the  pore  space  volume) 
stayed  intact  (Figure  8).  The  observed  effect  follows  from  our  theory  of  cemented  granular 
materials:  stress  concentration  is  high  at  the  contacts  of  uncemented  grains,  whereas  even 
small  amounts  of  relatively  soft  cement  result  in  a  more  uniform  stress  distribution  over  a 
larger  contact  area.  The  total  force  transmitted  between  two  grains  must,  of  course,  remain 
the  same  in  the  cemented  and  uncemented  cases;  however,  the  stresses  at  the  cemented 
interface  deerease  dramatically  as  compared  to  the  Hertzian  contact  case. 

The  contact  stress  distributions  as  obtained,  imply  the  following  modes  of  static  and 
dynamic  failure  of  the  grains  and  intergranular  bonds  in  a  particulate  material:  (1) 
Uncemented  grains  will  tend  to  fail  whereas  cemented  grains  will  stay  intact,  and  the 
cement  will  fail.  (2)  Where  intergranular  cementation  is  present,  grain  failure  may  still  be 
expected  if  the  cement  is  strong  and  stiff.  In  this  case,  grain  damage  will  be  initiated  at  the 
periphery  of  the  cement  layer.  (3)  Yielding  of  a  cement  material  that  is  soft  (as  compared  to 
the  grain  material)  will  initiate  at  the  center  of  the  contact  region,  whereas  stiff  cement  will 
yield  at  the  periphery. 
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Figure  8.  Porosity  versus  hydrostatic  confining  pressure  in  water-saturated  randomly  packed  glass 
beads  that  were  a)  uncemented  and  b)  cemented  by  epoxy  at  their  contacts  (Yin  and  Dvorkin,  1994). 
Porosity  was  measured  by  the  volume  of  expelled  fluid.  In  the  uncemented  case,  a  sharp  porosity 
decrease  is  observed  at  about  50  MPa.  The  decrease  is  associated  with  the  crushing  of  grains.  The 
cemented  grains  (the  volume  of  the  epoxy  accounted  for  only  10  percent  of  the  pore  space)  did  not 
crush.  The  photos  showed  that  in  the  latter  case,  the  grains  stayed  intact  with  the  failure  being 
localized  witWn  the  epoxy. 

CONCLUSIONS 

All  the  solutions  presented  in  this  paper  have  been  applied  to  the  interaetion  of  two 
circular  grains.  Note  that  the  obtained  contact  law  solutions  are  more  general  and  can  be 
applied  to  cemented  elastic  grains  with  arbitrarily-shaped  contact  surfaces  as  long  as  the 
underlying  assumptions  are  valid.  Specifically,  the  cement  layer  has  to  be  thin  and  small  as 
compared  to  a  grain,  and  in  addition  contacting  grain  surfaces  have  to  be  smooth,  so  that 
the  elastic  foundation  approximation  is  valid.  The  shape  of  grain  surfaces  will  affect  the 
thickness  of  the  cement  layer,  and  the  initial  stress  distribution  at  the  direct  contact  area  of 
precompacted  grains. 

In  the  example  presented  in  Figure  7  we  have  used  formulas  for  the  effective  bulk  and 
shear  moduli  of  a  random  packing  of  identical  spherical  grains  as  given  in  Digby  (1981) 
and  Winkler  (1983).  Again,  our  contact  laws  can  be  used  to  describe  the  effective 
properties  of  more  complicated  arrangements  of  cemented  grains.  Specifically,  these 
contact  laws  can  be  plugged  into  sophisticated  numerical  schemes  (e.g..  Distinct  Element 
Method). 

One  important  assumption  that  has  been  used  in  treating  the  problems  of  shear  and 
torsional  deformation  of  precompacted  cemented  grains  is  the  no-slip  condition  in  the  direct 
grain-to-grain  contact  zone.  This  assumption  is  valid  only  if  stresses  are  small,  so  that 
stress  distributions  obtained  under  the  no-slip  assumption  are  close  to  those  obtained  by 
considering  the  relative  sliding  of  contacting  grains  (see  examples  in  Johnson,  1992).  At 
the  same  time,  the  computed  shear  stress  distributions  shown  in  Figure  6  indicate  that  the 
sliding  will  initiate  at  the  periphery  of  the  direct  grain-to  grain  contact  area  only  if  the 
cement  is  very  soft  as  compared  to  the  grain  material.  In  most  cases,  we  do  not  expect  to 
encounter  high  stress  concentration  and  sliding  in  the  mentioned  area. 

The  presented  solutions  can  describe  the  dynamic  interaction  of  cemented  grains  only  if 
the  quasi-static  approximation  is  valid,  i.e.  the  wavelength  is  significantly  larger  than  the 
characteristic  microscopic  dimension  (the  grain  radius).  This  was  the  case  in  the  ultrasonic 
pulse  transmission  experiment  described  above:  the  wavelength  was  approximately  ten 
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times  the  average  grain  diameter.  These  solutions  can  be  also  applied  to  describing  shock 
wave  propagation  in  a  particulate  medium  if  the  length  of  the  pulse  is  much  larger  than  the 
grain  radius. 
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APPENDIX  1:  TWO  PRECOMPACTED  CYLINDERS 

Precompaction 

Consider  two  identical  elastic  circular  cylinders  (plane  strain)  of  radius  R  that  are 
normally  pressed  together  by  force  P^^  to  form  a  contact  area  of  half- width  a .  Then 


V  ;rG 


(Al.l) 


(Johnson,  1992).  The  initial  distribution  of  normal  stresses  p^^  on  the  surface  of  the 
precompacted  cylinders  is; 


Ph  =■ 


=  0,  |x|  >  a. 


^[a^  -  ,  |x|  <  a; 


(A1.2) 
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Here  x  is  the  coordinate  along  the  median  line  of  the  cement  layer  in  the  plane 
perpendicular  to  the  cylinder's  axis,  analogous  to  the  r  coordinate  in  the  three-dimensional 
case. 

By  treating  the  cylinder  as  an  elastic  half-plane,  we  have  the  following  expression  for 
the  initial  normal  displacements  v„(r)  of  the  cylinder's  surface  from  normal  stresses  p^ir) 
(Johnson,  1992): 


1-v 

nG 


b 


-b 


const. 


(A1.3) 


After  the  initial  compaction,  cement  is  uniformly  added  around  the  direct  contact  zone  to 
increase  the  half-width  of  the  contact  area  io  b.  It  follows  from  equation  (A  1.3)  that  the 
half-thickness  of  the  undeformed  cement  layer  between  two  precompacted  circular 
cylinders  is: 


h^{x)= 


-s^  In 

nR  •’ 


x-s 

a-s 


\ds-\- const,  a  <\x\<b. 


Normal  Deformation 


(A1.4) 


Consider  now  additional  normal  loading  of  the  cylinders  with  force  AP.  The  resulting 
total  normal  displacements  v(x)  of  the  cylinder  surface  within  the  area  of  the  initial  direct 
contact  are: 


v(x)-  v(iz)  = 


,  |x|  <  a. 


(A1.5) 


Within  the  cemented  zone,  a<\x\<b,  the  total  normal  displacements  of  the  grain 
surface  v(jc)  are  related  to  the  normal  displacements  of  the  cement  as: 

V(x)  =  v(x)  -  via)  -  [v^(x)  -  v^(a)],  a<\x\<  b,  (A1.6) 

where  V(r)  is  the  normal  displacement  of  the  cement  surfaee  (at  the  cement-cylinder 
interface)  relative  to  the  median  line  of  the  cement  layer. 

By  treating  the  cement  layer  as  an  elastic  foundation,  we  have  the  following  expression 

for  normal  stresses  pix)  at  the  cement-grain  interface: 


pix)  =  - 


2G,(1-VJ  Vjx) 

l-2v^  hfjix) 


a  <  |x|  <  b. 


(A1.7) 


Equation  (A  1.3)  now  takes  the  following  form: 

1- 

v(x)= - -  [p(5)ln|x-5l(i^+ const.  (A1.8) 

Now  we  combine  equations  (A1.5),  (A1.6),  (A1.7),  and  (A1.8)  to  obtain  the  following 
integral  equation  for  normal  stress  pix) : 
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v(x)-  v(a)  =-  — -[  j -  s\ds-  J /7(^)ln|x  -  s\ds]  = 


nG 


2  2 
a  —  X 


2R 


X  <a\ 


(A1.9) 


l-2v 


2Gn-vJ 


pix)h^  (x),  a  <  l^:]  <  b. 


Let  us  now  present  the  normal  stress  pix)  as  the  sum  of  the  initial  Hertzian  stress 
Ph(x)  and  the  additional  normal  stress  f(x): 

Pix)  =  Pii(x)  +  fix).  (ALIO) 

It  follows  then  from  equations  (AL2),  (A1.3),  (A1.9),  and  (ALIO)  that 


b 

J  /(5)ln|x  -  s\ds+  const  = 

-b 


[O,  |x|<<2; 


-  kG  1-2v, 
l-v2G,(l-v,) 


f(x)hfj(x). 


a  <  <  b. 


(Al.ll) 


This  last  equation  can  be  solved  using  the  numerical  quadrature  method. 

The  constant  in  the  left-hand  part  of  equation  (Al.l  1)  can  be  found  from  the  resulting 

force  AP  : 


b 


AP  =  J  f(x)dx. 

-b 

(A1.12) 

Shear  Deformation 

In  the  two-dimensional  case  under  consideration,  the  treatment  of  the  shear  deformation 
of  two  cemented  grains  is  identical  to  that  of  the  normal  deformation:  Equation  (AL8)  is 
replaced  by: 

1  V  ^ 

u{x)= - f^(5')ln|jr-5|ci5-i- const; 

(A1.13) 

equation  (A  1.5)  is  replaced  by 

u{x)  -  u(a)  =0,  |x|  <  a; 

(A1.14) 

equation  (A  1.6)  is  replaced  by 

U(x)  =  u{x)  -  u{a),  a  <  |x|  <  b; 

(A1.15) 
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and  equation  (A  1.7)  is  replaced  by 


q{x)  =  -G,^,a<\x\<b.  (A1.16) 

h„(x) 

The  definitions  of  the  tangential  displacements  u(x)  of  the  cylinder's  surface,  tangential 
displacements  C/(x)  of  the  cement,  and  shear  stresses  q{x)  in  contact  zone  are  the  same  as 
in  the  spherical  case.  Again,  as  in  the  spherical  case,  we  assume  the  no-slip  condition  in 
the  direct  contact  zone.  This  assumption  is  valid  if  shear  stresses  are  small. 

Similar  to  the  case  of  normal  deformation,  we  arrive  at  the  following  integral  equation 

for  stresses  q(x) : 

b 

J  ^(5)ln|x  -  const  = 

-b 

0,  (A1.17) 

-^^^;^q{x)hft{x),  a<\x\<b, 

where  the  constant  in  the  left-hand  part  of  equation  can  be  found  from  the  resulting  shear 
force  Q: 


b 

Q=  jq(x)dx.  (A1.18) 

-b 

Normal  and  shear  stress  distributions  in  this  two-dimensional  case  are  qualitatively 
similar  to  those  obtained  for  two  precompacted  spheres. 


APPENDIX  2;  TORSIONAL  DEFORMATION  OF  TWO  SPHERICAL 
GRAINS 

Uncompacted  Grains 

Consider  a  twisting  moment  M  applied  to  the  system  of  two  cemented  spherical  grains 
(Figure  1).  Shear  stress  components  q{r)  acting  on  the  grain-cement  interface  are 
axisymmetrical,  and  the  only  non-zero  displacement  of  the  grain  surface  u(r)  is 
perpendicular  to  the  polar  radius  r  (Figure  A2.1). 

The  displacements  u{r)  are  related  to  stresses  q{r)  as  (Johnson,  1992): 


M(r)= 


rcos 


J 


^jr^+s^-2rsco5(p 


(A2.1) 


where  integration  is  inside  the  circle  |r|  <  a  in  the  plane  of  contact  (see  Figure  4). 
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Figure  A2.1.  Torsional  displacements  and  shear  stresses  in  the  plane  of  grain-cement-grain  contact. 
The  problem  under  consideration  is  treated  in  a  linear  approximation,  so  we  assume  the  contact  zone  is 
planar. 

The  tangential  displacement  of  the  grain  surface  u(r),  and  the  tangential  displacement 
of  the  cement  layer  surface  relative  to  its  median  plane  U(r)  are  related  to  the  respective 
angles  of  rotation  j8(r)  and  B(r)  as: 

u(r)  =  rpir),  U(r)  =  rB(r).  (A2.2) 

The  condition  of  displacement  compatibility  between  the  grain  and  the  cement  surface 
is: 


{3(r)  =  B(r)  +  r,  (A2.3) 

where  7  represents  the  grain’s  rigid  rotation. 

Treating  the  cement  layer  as  an  elastic  foundation  in  torsion,  we  have  the  following 
relation  between  stress  q{r)  and  angle  B(r): 


^  h{r)  h(r) 


(A2.4) 


where  h(r)  is  the  half-width  of  the  cement  layer.  For  a  cement  layer  between  two  spheres: 


h(r)  =  R[e  +  -{-fls  =  ^ 
2  R  R 


(A2.5) 


where  is  the  minimum  half-separation  between  the  spheres. 

Finally,  equations  (A2.1),  (A2.2),  (A2.3),  (A2.4)  ,  and  (A2.5)  yield  the  following 
integral  equation  for  B(r) : 


r[7-bB(r)]= 


rcos<p+^ 


’’  B{-\jr^+s^-2rscos(p)(r- scos(p) 


kG 


\d(p  J 


1  /  S  M 


ds. 


(A2.6) 


The  constant  7  in  equation  (A2.6)  can  be  found  from  the  resulting  twisting  moment 


M: 
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(A2.7) 


The  resulting  distributions  of  the  twisting  shear  stresses  along  the  radius  of  a  cement 
layer  are  qualitatively  similar  to  those  shown  in  Figure  2. 

Precompacted  Grains 

The  treatment  of  this  problem  is  analogous  to  that  used  in  the  problems  of  normal  and 
shear  deformation  of  two  precompacted  grains.  The  resulting  integral  equation  for  the 
twisting  shear  stresses  q(r)  in  the  contact  zone  is: 


,  rcosW^^-.^sin^^^(^^2^^2_2^^COS.p)(r-^COS<p) 


—jd(p  j 

yr,  r<a; 


^jr^+s^—2rscos(p 


ds  = 


(A2.8) 


The  constant  7  here  can  be  found  from  the  resulting  twisting  moment: 


b 

M=27t^q{r)r^  dr.  (A2.9) 

0 

The  distributions  of  twisting  shear  stresses  are  similar  to  those  shown  in  Figure  6. 

This  solution  uses  the  no-slip  condition  in  the  direct  grain-grain  contact  zone.  The 
assumption  is  valid  if  shear  stresses  are  small. 
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2.3.  STRENGTH  OF  CEMENTED  GRAINS 


ABSTRACT 

We  conducted  compaction  tests  (isotropic  drained  loading)  on  randomly  packed  glass 
beads  that  were  a)  uncemented  and  b)  cemented  by  epoxy  at  their  contacts.  In  the  latter 
case,  the  volume  of  the  epoxy  accounted  for  10  percent  of  the  pore  space.  Intensive 
crushing  of  grains  was  observed  in  the  first  case  at  about  50  MPa.  In  the  second  case,  the 
cemented  grains  stayed  intact,  the  failure  being  localized  within  the  epoxy.  Therefore,  even 
small  amounts  of  cement  at  contacts  prevent  the  failure  of  grains.  Theoretically,  this  effect 
follows  from  our  theory  of  cemented  granular  materials:  stress  concentration  is  high  at  the 
contacts  of  uncemented  grains,  whereas  even  small  amounts  of  relatively  soft  cement  result 
in  a  more  uniform  stress  distribution  over  a  larger  contact  area. 

INTRODUCTION 

The  elastic  moduli  and  the  strength  of  granular  media  in  general,  and  sedimentary  rocks 
in  particular,  may  be  strongly  affected  by  the  properties  and  structure  of  intergranular  bond 
material.  The  rocks  of  interest  span  a  wide  range  from  diagenetic  sediments  and  sand-clay 
mixtures  to  tar  sands. 

In  spite  of  significant  progress  achieved  in  modeling  macroscopic  deformation  of 
granular  media  (e.g.,  Cundall  and  Strack,  1979;  Nemat-Nasser,  1983;  Nelson  et  al.,  1988; 
Cundall  et  al.,  1989),  a  challenge  remains  in  predicting  contact  stiffness  and  strength  from 
the  geometry  of  a  contact  region  and  from  the  mechanical  properties  of  grains  and 
intergranular  cement.  Almost  all  theoretical  investigations  of  deformation  of  granular 
aggregates  are  based  on  the  classical  solutions  to  the  problems  of  normal  (Hertz,  1882)  or 
oblique  (Mindlin,  1949;  Walton,  1978)  interaction  of  elastic  spheres.  Digby  (1981)  solved 
the  problem  of  interaction  of  spherical  particles  that  are  initially  bonded  together  across 
small  areas.  Trent  (1989),  and  Trent  and  Margolin  (1992)  numerically  investigated  the 
behavior  of  cemented  granular  materials  under  both  low-strain  and  high-strain  loads.  The 
materials  are  composed  of  circular  particles  which  are  glued  together  with  elastic  bonds. 
They  show  that  the  effective  properties  of  samples  are  governed  by  the  properties  and 
distribution  of  individual  intergranular  bonds. 

Significant  attention  has  been  paid  to  the  experimental  investigation  of  the  geotechnical 
properties  of  cemented  rocks  (e.g.,  Gallagher  et  al.,  1974,  Clough  et  al.,  1981,  Saxena, 
1982).  Rad  and  Clough  (1982)  experimentally  quantified  the  influence  of  cementation  on 
the  static  and  dynamic  behavior  of  sands.  Their  conclusions  state  that  cementation  can  be  a 
very  important  parameter  relative  to  the  liquefaction  resistance  of  sand.  A  recent 
experimental  study  by  Bernabe  et  al.  (1992),  conducted  on  synthetic  cemented  granular 
materials  (precompacted  Ottawa  sand  with  halite  and  silica  glass  as  the  cement),  included 
triaxial  compression  tests.  It  has  been  found  that  a  small  amount  of  cement  can 
significantly  increase  the  strength  of  a  granular  material  if  it  is  deposited  precisely  at 
previously  formed  grain-to-grain  contacts.  The  main  conclusion  of  this  work  is  that  the 
cement  strengthens  the  rock  by  preventing  sliding  and  rotation  of  grains. 

A  different  (from  grain  sliding  and  rotation)  mode  of  rock  failure  is  related  to  grain 
crushing.  Gallagher  et  al.  (1974)  showed  that  because  grains  in  rock  aggregates  are 
usually  constrained  at  higher  confining  pressures  and  by  stronger,  stiffer  cements  than 
normally  exists  in  soils,  grain  fracturing  and  crushing  in  rocks  cannot  be  neglected.  Bruno 
and  Nelson  (1991)  numerically  investigated  rock  failure  under  tensile,  uniaxial 
compression,  and  biaxial  loading.  Their  conclusion  is  that  internal  grain  fracturing  may 
affect  the  strength  of  the  rock  at  high  stresses.  This  theoretical  result  is  supported  by 
photographs  showing  internally  cleaved  grains. 
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Dvorkin  et  al.  (1991),  theoretically  examined  the  normal  interaction  of  two  spherical 
elastic  grains  and  an  elastic  cementation  layer  between  them.  Their  results  show  that  even 
with  a  small  amount  of  cement,  a  contact  force  is  distributed  over  an  area  that  is  larger  than 
that  in  the  uncemented  case.  Contact  stress  distribution  becomes  more  uniform  as 
compared  to  the  uncemented  case.  These  two  factors  lead  to  a  significant  reduction  in 
stress  concentration  between  cemented  grains.  Even  very  soft  cement  bears  the  load  at 
grain  contacts  thus  reducing  stress  concentration  and  preventing  grain  crushing. 

Both  theoretical  and  experimental  results  indicate  that  grain  crushing  does  not  normally 
occur  at  low  stress  level,  but  becomes  an  increasingly  significant  factor  as  stresses 
increase.  This  mode  of  rock  failure  may  contribute  to  such  instabilities  as  sanding  and 
breakouts  in  deep  wellbores.  This  is  why  we  consider  to  be  important  to  investigate  the 
internal  strength  of  grains  in  a  granular  material. 

In  this  paper,  we  describe  our  compaction  tests  (isotropic  drained  loading)  on  cemented 
and  uncemented  glass  beads.  Uncemented  beads  were  crushed  (brittle  failure)  at  about  50 
MPa.  The  failure  was  associated  with  sharp  porosity  reduction  and  intensive  acoustic 
emission.  The  grains  appeared  crushed  on  the  photo  taken  after  the  loading  cycle. 
Cemented  grains  were  bonded  together  by  epoxy  at  10  percent  epoxy  saturation  of  the  pore 
volume.  These  grains  did  not  fail  during  the  loading  up  to  80  MPa.  The  photo  taken  after 
the  loading  cycle  shows  that  the  grains  stayed  intact  and  only  the  epoxy  cement  failed. 

We  explain  this  effect  from  our  theory  of  cemented  granular  materials:  even  small 
amounts  of  relatively  soft  cement  result  in  a  more  uniform  stress  distribution  over  larger 
contact  areas,  as  compared  to  the  high  stress  concentration  at  small  contact  areas  between 
uncemented  grains. 

Our  conclusion  is  that  even  small  amounts  of  cement  deposited  exactly  at  intergranular 
contacts  significantly  increase  the  strength  of  the  grains.  Cementing  grains  appears  to  be  an 
effective  way  of  creating  strong  high-porosity,  high-permeability  materials.  Such  an 
approach  may  be  used  in  preventing  grain  crushing  for  stabilizing  the  bottom  hole 
formation  zone  of  a  sand-producing  well,  or  in  manufacturing  artificial  aggregates. 

EXPERIMENTS 

In  our  experiments,  we  used  spherical  glass  beads  of  0.42  -  0.5  mm  diameter.  The 
material  contains  no  more  than  two  percent  of  irregularly-shaped  grains,  and  is  reasonably 
free  of  sharp  angular  particles,  or  particles  showing  milkiness,  surface  scoring,  and  foreign 
matter.  The  porosity  of  the  uncemented  glass  bead  sample,  as  measured  before  loading, 
was  about  39  percent  -  close  to  the  porosity  of  a  dense  random  pack  of  identical  spheres. 
The  cemented  sample  was  prepared  by  mixing  glass  beads  with  a  given  volume  of  epoxy 
(in  this  case  it  was  10  percent  of  the  pore  space  volume).  The  mixture  was  placed  into  a 
holder  where  the  epoxy  solidified.  The  liquid  epoxy  is  a  wetting  fluid  and  thus  it 
accumulated  at  grain  contacts. 

The  bulk  and  shear  moduli  of  the  glass  as  calculated  from  the  manufacturer’s 
specifications  are  50.0  GPa  and  26.2  GPa,  respectively.  Those  of  the  solidified  epoxy 
(calculated  from  compressional  and  shear  wave  velocities,  and  density)  are  6.8  GPa  and 
2.0  GPa,  respectively. 

Isotropic  drained  loading  at  confining  pressure  up  to  80  MPa  was  conducted  in  a 
triaxial  device  on  water  saturated  unjacketed  samples.  During  the  loading,  porosity  was 
measured  through  the  volume  of  expelled  pore  water.  We  estimate  the  absolute  error  of 

these  measurements  at  about  ±2  percent. 

The  porosity  of  the  cemented  sample  showed  a  smooth  decrease  with  increasing 
confining  pressure,  whereas  the  porosity  of  the  uncemented  sample  sharply  decreased  from 
35  to  31  percent  at  50  MPa  (Figure  1).  This  sharp  decrease  was  probably  associated  with 
the  crushing  of  glass  beads.  Another  indication  of  the  crushing  was  the  intensive  acoustic 
emission  registered  at  about  50  MPa. 
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Confining  Pressure  (MPa) 

Figure  1.  Porosity  of  glass  beads  versus  pressure:  open  circles  --  cemented  beads,  empty  circles  - 
uncemented  beads.  Porosity  was  measured  by  the  volume  of  pore  fluid  expelled. 

The  photo  of  the  uncemented  grains  taken  after  the  loading  cycle  showed  their  intensive 
brittle  crushing.  On  the  other  hand,  the  cemented  grains  did  not  fail.  This  experiment 
showed  that  cementing  fragile  grains  is  an  effective  way  of  preventing  their  failure  at  high 
loading. 

THEORY 

Our  theory  of  cemented  granular  materials  is  based  on  the  solution  of  the  contact 
problem  of  two  cemented  elastic  grains  (Dvorkin  et  al.,  1991).  The  cement  material  is 
elastic  and  may  be  different  from  the  grain  material.  There  are  two  main  conclusions 
following  from  the  solution. 

The  first  one  is  quite  obvious:  even  a  small  amount  of  cement  at  a  grain  contact  creates 
a  contact  zone  that  is  much  larger  than  that  created  by  pressing  together  two  uncemented 
grains.  For  example,  at  10  percent  epoxy  saturation  of  the  pore  space,  the  ratio  of  the 
radius  of  the  cement  layer  a  (Figure  2a)  to  grain  radius  is  0.43.  At  the  same  time,  the  ratio 
of  the  radius  of  the  contact  area  between  uncemented  grains  (the  Hertz  solution)  to  grain 
radius  is  0. 1 1  at  50  MPa  confining  pressure. 

The  second  conclusion  is  that  the  cement  is  load  bearing  in  the  case  where  there  are 
only  point  glass/glass  contacts.  Normal  contact  stress  (Figure  2b)  is  maximum  at  the  center 
of  the  contact  area  when  cement  material  is  much  softer  than  grain  material,  and  is 
maximum  at  the  periphery  of  the  contact  area  when  the  cement  is  stiff. 

An  example  of  calculated  normal  contact  stresses  between  two  glass  beads  that  are  a) 
epoxy-cemented,  and  b)  uncemented  is  given  in  Figure  3a.  Contact  stresses  at  the 
cemented  interface  are  more  uniformly  distributed  than  those  for  the  uncemented  case.  The 
maximum  stress  at  the  uncemented  interface  is  13.6  times  as  large  as  that  at  the  cemented 
interface.  Figure  3b  shows  a  hypothetical  case  where  the  cement’s  bulk  and  shear  moduli 
are  four  times  smaller  than  those  of  the  epoxy.  In  this  case  we  can  clearly  see  a  maximum 
of  the  contact  stresses  at  the  center  of  the  contact  area.  Even  in  this  case  of  very  soft 
cement,  this  maximum  stress  is  about  0.32  of  that  at  the  cemented  interface. 

These  theoretical  calculations  explain  the  experimentally  observed  effect  of  increasing 
strength  of  grains  with  cement  at  the  contacts. 
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Figure  2.  Two  cemented  grains:  a,  cement  forms  a  circular- shaped  contact  area;  b.  normal  contact 
stresses  in  cement. 
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Figure  5.  Theoretical  prediction  of  stresses  at  grain  contacts  for  cemented  and  uncemented  beads:  a. 
epoxy-cemented  glass  beads,  b.  glass  beads  with  cement  whose  bulk  and  shear  moduli  are  four  times 
smaller  than  those  of  epoxy.  Stresses  are  normalized  by  the  average  stress  at  the  cemented  contact. 
Coordinate  along  cement  layer  is  normalized  by  the  radius  of  the  cement  layer. 

CONCLUSIONS 


We  show  experimentally  and  explain  theoretically  that  even  soft  intergranular  cement,  if 
deposited  at  grain  contacts,  significantly  increases  the  strength  of  individual  grains.  This 
effect  is  due  to  a  more  uniform  distribution  of  contact  stresses  in  the  load-bearing  cement  as 
compared  to  a  high  stress  concentration  at  the  uncemented  interface.  This  effect  holds  even 
for  the  amount  of  cement  as  low  as  10  percent  of  pore  space  volume.  Therefore,  cementing 
grains  provides  an  effective  means  of  creating  strong  high-porosity,  high-permeability 
materials. 
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2.4.  Seismic  detection  of  residual  Contaminants 


ABSTRACT 

The  key  to  effective  characterization  and  treatment  of  contaminated  sites  is  our  ability  to 
delineate  the  spread  of  contaminants  in  the  shallow  subsurface.  A  promising  non-invasive 
technique  that  allows  one  to  identify  contaminants  is  the  geophysical  three-dimensional 
mapping  of  seismic  velocities,  reflections,  and  attenuation.  The  physical  principle  of  this 
technique  is  that  the  same  rock,  if  filled  with  different  fluids,  transmits  sound  waves 
differently.  There  are  two  main  reasons  for  this  effect:  (a)  pore-fluid  compressibility  and 
(b)  pore-fluid  viscosity.  While  the  effect  of  pore-fluid  compressibility  has  been  extensively 
used  for  hydrocarbon  detection  (bright  spots),  the  latter,  pore-fluid  viscosity,  effect  has 
been  neglected,  mainly  because  it  is  noticeable  at  high  frequencies  only.  We  present 
theoretical  models  which  show  that  at  relatively  high,  but  practically  important,  frequencies 
(1-10  kHz)  the  effect  of  pore-fluid  viscosity  becomes  seismically  visible  and  thus  allows 
one  to  locate  residual  viscous  contaminants  left  in  thin  pores  and/or  between  sand  grains. 
This  conclusion  is  supported  by  a  laboratory  ultrasonic  experiment  where  adding  small 
amounts  of  viscous  epoxy  to  a  glass-bead  pack  resulted  in  increasing  compressional-wave 
velocity. 

INTRODUCTION 

The  cost  of  treatment  of  contaminated  sites  may  be  sufficiently  reduced  if  the  location  of 
contaminants  is  accurately  predicted.  This  task  becomes  especially  complicated  if  there  are 
only  small  amounts  of  a  contaminant  left  in  thin  cracks  and  between  grains.  A  promising 
non-invasive  technique  that  allows  one  to  locate  contaminants  in  the  shallow  subsurface  is 
the  geophysical  mapping  of  seismic  velocities,  reflections,  and  attenuation. 

A  critical  missing  ingredient  in  previous  applications  of  seismic  methods  to  contaminant 
and  groundwater  problems  is  the  relations  between  the  geophysical  signatures  (seismic 
velocities  and  attenuation)  and  physical  properties  of  pore  fluid  —  its  viscosity  and 
compressibility.  Finding,  understanding  and  quantifying  these  relations  is  the  key  goal  of 
rock  physics. 

Rock  physics  has  been  successfully  used  in  the  oil  and  gas  industry  in  transforming 
high-resolution  three-  and  two-dimensional  maps  of  seismic  velocities  and  attenuation  into 
important  physical  properties  of  deep  reservoirs  and  pore  fluids.  One  important  effect  that 
allows  for  such  transformations  is  that  seismic  velocities  in  rocks  generally  change  with 
changing  pore  fluid,  and  with  changing  frequency.  There  are  two  principal  reasons  for  this 
velocity  versus  saturation  and  frequency  change:  (1)  Pore  fluid,  when  placed  inside  the 
rock,  acts  as  an  elastic  component  that  adds  stiffness  to  the  solid-fluid  composite,  and  (2) 
At  high  frequencies,  viscous  pore  fluid  is  unable  (does  not  have  time)  to  freely  flow 
between  highly  compliant  thin  pores  and  large  pores  (unrelaxed  state  of  pore  fluid). 
Therefore,  at  high  frequencies,  pore  fluid  is  trapped  inside  thin  pores,  thus  reinforcing  the 
most  compliant  elements  of  the  rock  (Figure  la). 

The  first,  pore-fluid-compressibility,  effect  has  been  widely  used  in  prospecting  for 
hydrocarbons:  often  compressional-wave  velocity  in  gas-saturated  or  oil-saturated 
reservoirs  is  smaller  than  in  surrounding  brine-saturated  rock.  Such  velocity  contrast 
results  in  a  bright  spot  that  directly  indicates  the  location  of  hydrocarbon  (e.g.,  Clark, 
1992).  The  second,  frequency  or  pore-fluid- viscosity,  effect  cannot  be  used  in 
hydrocarbon  reservoirs  because  high-frequency  waves  cannot  penetrate  into  the  deep 
subsurface. 

One  important  advantage  of  applying  seismic  imaging  to  environmental  problems  is  that 
high-frequency  signals  can  be  used  to  characterize  the  shallow  subsurface.  For  example,  a 
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frequency  of  1  kHz  is  routinely  employed  to  obtain  high-resolution  velocity  maps  in  well- 
to-well  tomography. 


Figure  1.  a.  Contaminant  (black)  present  in  thin  cracks,  b.  Contaminant  (black)  coating  sand  grains 
and  creating  cementation  at  the  grain  contacts,  c.  Contaminant  filling  the  gap  at  a  grain-to-grain 
contact. 

In  this  paper,  we  show  that  high-frequency  seismic  imaging  can  be  used  to  detect 
residual  viscous  contaminant  left  in  thin  cracks  and/or  at  grain  contacts.  Our  theoretical 
estimates  are  based  on  the  squirt-flow  mechanism:  Stiff  portions  of  the  pore  space  have 
relatively  small  variations  of  pore  pressure  when  the  rock  is  loaded  by  a  passing  wave. 
Soft  portions  of  the  pore  space  tend  to  transfer  more  of  the  stress  to  the  fluid,  which  results 
in  high  variations  of  induced  pore  pressure.  At  low  frequencies,  one  expects  these 
pressures  to  equilibrate  —  an  assumption  leading  to  Gassmann's  formula  for  the  bulk 
modulus  of  saturated  rock.  At  higher  frequencies,  intensive  cross-flow  between  the  soft 
and  the  stiff  parts  persists,  resulting  in  wave-energy  dissipation  and  velocity/frequency 
dispersion.  At  very  high  frequencies,  the  fluid  is  unrelaxed  and  blocked  in  the  compliant 
pores. 

We  show  that  if  the  contaminant's  viscosity  is  significantly  larger  than  that  of  water, 
noticeable  velocity  dispersion  may  occur  at  frequencies  of  several  kHz.  Therefore,  the 
effect  of  pore-fluid  viscosity  becomes  seismically  visible  allowing  one  to  locate  residual 
viscous  contaminant. 

SANDSTONES:  RESIDUAL  CONTAMINANT  IN  CRACKS 

The  following  theoretical  formulas  have  been  derived  based  on  the  squirt-flow  model 
for  wave-induced  motion  of  the  fluid  present  in  thin  cracks  (Dvorkin  et  ah,  1994).  They 
can  be  used  to  calculate  P-  and  S-wave  velocities  (V^  and  V^,  respectively),  and 

attenuation  (inverse  quality  factors  (2“*  and  ,  respectively)  in  fluid-saturated  rock  with 
residual  contaminant  left  in  the  thin,  compliant  portions  of  the  pore  space  (Figure  la).  In 
this  model,  the  contaminant  may  differ  from  the  saturating  fluid. 

V,=^Rc{K,+ia,)lp.  V.=VRe(G„)/p, 

Q;'  =  |lm  (K,  +  4  G.)|  /  [Re  (K,  +  4  G.)|, 

e;'=|lm(G.)|/|Re(G.)|. 

In  these  formulas 

=  K„,  /  (1  +  a„,  =1  -  K„,  /  K^; 
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l/G-l/G.,=(4/15)(l/X-l/i:^), 


llK.^  =  \lk.„  +  \IK,,-llK,. 

k.,  =  K,„^aK.V-fm. 


IIK=IIK..  +  IIK^-IIK,-, 


^ms  - 


a  =  \-K^J  K^,  M)  =  2  7,(^)  /  ^  Z- 


dP  r  M-1 


l/Fo  =  l/^/  +  l/(^a,),  a^  =  \-KIK^, 

Qo  =  KJ(a,-<t>); 

l/K^,  =  l/K^-l/K,,  +  UK-, 


where  p  is  the  saturated-rock  density;  K  and  G  are  the  dry-rock  bulk  and  shear  moduli, 
respectively;  K^^p  is  the  bulk  modulus  of  dry  rock  at  high  pressure;  is  the  bulk  modulus 
of  the  mineral  phase;  is  the  pore-fluid  bulk  modulus;  p  is  the  dynamic  viscosity  of  the 
contaminant;  ^  is  porosity;  and  Jq  and  J^  are  Bessel  functions  of  zero  and  first  order, 

respectively.  All  above  quantities  are  in  SI  units.  Parameter  Z  describes  the  wave-induced 
squirt  flow  of  the  contaminant  in  thin  cracks  and  is  frequency-  and  viscosity-independent. 

It  can  be  found  by  matching  one  of  the  four  parameters:  V^,  ,  Q~\  or  with  one 

experimental  measurement  (preferably  V^).  The  experiment  can  be  performed  on  a  clean 

water-saturated  sample  at  ultrasonic  frequency.  Typically,  Z  ranges  from  0.03  to  0.09. 

We  used  these  formulas  to  calculate  P-  and  S-wave  velocities  in  water-saturated  Navajo 
sandstone  of  1 1.8  percent  porosity.  The  viscosity  of  the  contaminant  varied  from  1  cPs  to 
10^  cPs.  The  Z  value  calculated  from  ultrasonic-pulse-transmission  experiments  (Coyner, 
1984)  was  0.087. 

Compressional-wave  velocity  versus  frequency  is  shown  in  Figure  2a.  We  can  see 
significant  velocity/viscosity  dispersion  at  1  kHz  -  a  frequency  commonly  used  in  well-to- 
well  tomography.  If  the  contaminant's  viscosity  is  as  low  as  100  cPs  —  the  viscosity  of  30 
deg.  API  oil  (Batzle  and  Wang,  1992),  the  difference  between  the  clean  and  the 
contaminated  rock  is  about  120  m/s.  This  difference  approaches  290  m/s,  and  330  m/s  as 
the  contaminant's  viscosity  becomes  10^  cPs  and  10^  cPs,  respectively.  A  strong 
velocity/viscosity  dispersion  at  1  kHz  can  be  also  observed  for  shear-wave  velocity  (Figure 
2b).  Here  V,  increases  from  2800  m/s  to  2850  m/s,  2905  m/s,  and  2915  m/s  as  the 
contaminant's  viscosity  increases  from  1  cPs  to  10^  cPs,  10^  cPs,  and  10^  cPs, 
respectively. 

Based  on  these  estimates,  we  conclude  that  even  small  amounts  of  viscous  contaminant 
left  in  the  thin,  compliant  portions  of  the  pore  space  can  be  seismically  detectable  at 
frequency  1  kHz  and  larger. 
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Figure  2.  Compressional-and  shear-wave  velocities  in  water-saturated  Navajo  sandstone.  Viscous 
contaminant  is  present  only  in  thin,  compliant  cracks.  The  numbers  indicate  the  contaminant's 

viscosity,  a.  Vp.  b  V^. 

SANDS:  RESIDUAL  CONTAMINANT  AT  GRAIN  CONTACTS 


Consider  unconsolidated  sand  grains  coated  with  viscous  contaminant  (Figure  lb). 
The  rest  of  the  pore  space  may  be  dry  or  water-saturated.  As  a  result  of  such  pore-space 
topology,  the  contaminant  fills  the  gaps  around  grain-to-grain  contacts.  Loading  of  the 
grains  by  a  passing  seismic  wave  results  in  the  squirting  flow  of  the  contaminant  in  these 
contact  gaps  (Figure  Ic).  If  frequency  (and/or  viscosity  of  the  contaminant)  is  high,  the 
fluid  in  a  contact  gap  will  resist  the  squeezing  action  of  the  seismic  excitation  by  developing 
substantial  hydrodynamic  pressure  and  creating  contact  stiffness.  We  calculate  this 
stiffness  by  combining  a  model  of  fluid  flow  in  a  crack  (Dvorkin  et  al.,  1992)  and  a  model 
of  the  normal  deformation  of  two  cemented  grains  (Dvorkin  et  al.,  1991). 

To  obtain  the  stiffness,  one  has  to  numerically  solve  the  following  integro-differential 
equation; 


.  “i, tanh’ l)-P^  = 


X 


S(p+-sl^-r"  sin"  (p 


2/?^  Po  |-1 — -  jd(p  JP(^/r^^i-7^^-2r^cos^)d5-5],  A 

f  TCG  n  n 


2Ri~n 
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where  P  is  hydrodynamic  pressure  in  the  fluid;  (O  is  angular  frequency;  is  the 
contaminant's  density;  rj  is  its  kinematic  viscosity;  Cq  is  its  acoustic  velocity;  R  is  the 

grain's  radius;  G  and  v  are  its  shear  modulus  and  Poisson's  ratio,  respectively;  b  is  the 
radius  of  the  contaminant-filled  gap  between  two  grains,  as  measured  along  the  r 
coordinate  (Figure  Ic);  and  6  is  an  arbitrarily-chosen  non-zero  constant.  The  b  value  can 
be  calculated  from  the  thickness  of  the  contaminant  layer  h  as  b  =  ^2Rh.  The  equation 
has  to  be  solved  with  the  boundary  conditions  P  =  0  at  r  =  b,  and  P'  =  0  at  r  =  0 . 

After  the  equation  is  solved,  the  normal  two-grain  stiffness  can  be  calculated  from 

b 

5„=P/5,  F  =  j2;rrP(r)Jr. 

0 

Finally,  assuming  that  the  shear  stiffness  is  small  as  compared  to  the  normal  stifness, 
the  effective  bulk  and  shear  moduli  of  a  random  pack  of  identical  spherical  grains,  and 
Ggff,  can  be  calculated  from 

12;rP  207tR 

where  C  is  the  coordination  number,  and  0  is  porosity  (Winkler,  1983).  This  dry- 
frame  moduli  can  be  used  in  Gassmann's  equation  when  calculating  P-wave  velocity  in 
water-saturated  sand. 

We  calculated  P-wave  velocity  versus  frequency  and  viscosity  in  unconsolidated  water- 
saturated  glass  beads  coated  with  viscous  contaminant  (Figure  3).  In  this  case 
velocity/viscosity  dispersion  is  not  as  large  as  in  Navajo  sandstone.  For  example,  at 
frequency  1  kHz  the  velocity  increase  between  clean  beads  and  the  beads  coated  with  10^ 
cPs  viscosity  contaminant  is  only  48  m/s.  This  difference  reaches  130  m/s  at  frequency  10 
kHz.  We  conclude  that  in  unconsolidated  sands,  residual  contaminant  is  harder  to  detect 
than  in  Navajo-type  sandstones.  Nevertheless,  using  high-frequency  high-resolution 
seismic  techniques  may  help  locate  this  type  of  residual  contamination  caused  by  high- 
viscosity  substances. 


Figure  3.  Compressional-wave  velocity  in  unconsolidated  water- saturated  glass  beads.  The  beads  are 
coated  with  viscous  contaminant.  The  ratio  h  /  R  is  0.05.  The  numbers  indicate  the  contaminant's 
viscosity. 
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EXPERIMENTS 


We  have  conducted  ultrasonic  pulse-transmission  experiments  (1  MHz)  on  dry  glass 
beads  partially  saturated  with  epoxy.  The  viscosity  of  the  epoxy  was  about  400  cPs.  The 
confining  pressure  varied  from  5  to  20  MPa.  The  microphotographs  of  the  thin  sections 
showed  that  the  epoxy  filled  the  intergranular  contact  gaps.  The  results  (Figure  4)  indicate 
that  there  is  a  small  but  noticeable  velocity  increase  associated  with  the  addition  of  epoxy. 
It  is  interesting  that  this  increase  (about  40  m/s)  only  weakly  depends  on  saturation  (if 
saturation  is  partial).  The  reason  is  that  the  most  intensive  squirt  flow  occurs  in  the  direct 
vicinity  of  the  contact  and  thus  adding  the  epoxy  does  not  significantly  change  wave- 
induced  pressure  in  the  contact  gap.  This  effect  is  also  predicted  by  our  theoretical 
estimates. 


O  Vp  at  5  MPa 

□  Vp  at  10  MPa 

•  Vp  at  15  MPa 

A  Vp  at  20  MPa 


Figure  4.  Ultrasonic  compressional-wave  velocity  in  compacted  partially  epoxy-saturated  glass  beads. 

CONCLUSIONS 

Small  amounts  of  viscous  contaminant,  if  located  in  thin  compliant  pores,  or  at  grain 
contacts,  may  noticeably  affect  high-frequency  (several  kHz)  seismic  velocities  in  rocks. 
This  effect  is  large  in  Navajo-type  sandstones  and  is  small  in  unconsolidated  sands. 
Velocity  increase  in  contaminated  rocks  (as  compared  to  clean  rocks)  strongly  depends  on 
the  viscosity  of  the  contaminant:  the  larger  viscosity,  the  larger  the  increase  (at  a  given 
frequency). 

This  velocity/viscosity  dispersion  may  serve  as  the  theoretical  foundation  for  a 
promising  non-invasive  geophysical  technique  that  allows  one  to  locate  residual 
contaminants  in  the  shallow  subsurface. 
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2.5.  ELASTICITY  OF  PARTIALLY  SATURATED  FROZEN  SAND 


INTRODUCTION 

Seasonal  thawing  and  freezing  of  near-surface  sediments  significantly  affects  the 
interpretation  of  seismic  reflection  surveys  and  vertical  seismic  profiling,  as  well  as  activity 
on  the  surface.  Permafrost  covers  much  of  the  Earth's  colder  regions  and  is  also  subject  to 
periodic  thawing  and  freezing.  The  significance  of  the  regions  with  large  seasonal 
temperature  variations  for  exploration  and  engineering  activities  makes  understanding  the 
elasticity  of  frozen  ground  an  important  practical  and  scientific  goal. 

Previous  theoretical  and  experimental  studies  of  seismic  velocities  in  frozen  rocks  (e.g., 
Timur,  1968)  showed  that  P-wave  velocity  in  fully  water-saturated  rocks  increases  with 
decreasing  temperature  whereas  it  is  almost  independent  of  temperature  in  dry  rocks.  As 
temperature  decreases,  saturated  rocks  become  (given  enough  time)  completely  frozen  and 
velocity  reaches  its  constant  (terminal)  value.  A  theoretical  scheme  for  estimating  both  P- 
and  S-wave  velocities  in  saturated  freezing  rocks  has  been  proposed  by  Zimmerman  and 
Kang  (1986);  the  method  accurately  predicted  the  experimental  data. 

In  this  study  we  are  concerned  with  velocities  in  partially  saturated  unconsolidated 
sands.  Our  experiments  show  that  in  hydrophilic  granular  materials  velocities  may  be  close 
to  their  full-saturation  terminal  values  at  saturations  as  small  as  10  - 15  percent.  The  reason 
for  this  strong  velocity  increase  from  the  dry  state  is  the  cementing  action  of  ice  that  forms 
pendular  rings  at  the  grain  contacts.  This  freezing  effect  is  a  special  case  of  a  more  general 
phenomenon:  even  small  amounts  of  relatively  soft  cement,  if  placed  at  the  grain  contacts, 
act  to  dramatically  increase  the  stiffness  of  the  aggregate. 

By  applying  the  Dvorkin  et  al.  (1994)  intergranular  cementation  theory  we  accurately 
predict  the  measured  experimental  values.  These  accurate  predictions  cannot  be  achieved 
by  effective  medium  theories  that  do  not  take  into  consideration  the  specific  location  of  the 
cement  -  at  the  grain  contacts. 

EXPERIMENTS 

The  laboratory  apparatus  includes  a  76  cm-long,  7.6  cm  in  diameter  Plexiglas  holder 
filled  with  the  dry-matrix  material  (0.7  mm  in  diameter  well-sorted  glass  beads,  or  Ottawa 
sand  with  average  grain  size  0.3  mm).  P-  and  S-wave  velocities  are  measured  within  the 
holder  using  transmitter  and  receiver  transducers  (both  2.5  MHz  central  frequency)  placed 
1.52  cm  apart.  The  transducers  are  held  within  a  rigid  plastic  frame,  which  is  lowered  into 
the  Plexiglas  tube  by  the  cables.  The  holder  is  seated  on  a  metal  base  with  drainage  holes. 

The  matrix  materials  were  saturated  by  suctioning  degassed  water  up  the  tube  through 
the  holes  in  the  base.  This  technique  allowed  air  to  escape  and  prevented  air  bubbles  from 
lodging  in  the  matrix.  Partial  saturation  was  achieved  by  draining  the  fully-saturated  matrix 
materials.  The  amounts  of  water  held  by  capillary  forces  between  the  grains  were  umform 
along  the  upper  part  of  the  holder  and  accounted  for  13.1  and  13.5  percent  saturation  in 
glass  beads  and  Ottawa  sand,  respectively.  Both  matrix  materials  appeared  to  be 
hydrophilic  -  water  formed  pendular  rings  around  grain  contacts. 

P-  and  S-wave  velocities  were  measured  in  the  samples  (dry,  saturated,  and  partially 
saturated)  after  they  spent  a  day  in  a  freezer  at  about  -20°  C  (to  ensure  the  complete  freezing 
of  water).  The  results  of  these  measurements,  as  well  as  the  experimental  values  for  dry- 
matrix  porosity  and  ice  saturation  of  the  pore  space  are  summarized  in  Table  1. 

The  measured  ultrasonic  velocities  are  plotted  versus  ice  saturation  in  Figure  1 .  At  zero 
ice  saturation  both  P-  and  S-wave  velocities  are  zero  since  no  confining  pressure  is  applied 
to  the  samples.  Apparent  is  a  dramatic  increase  in  the  velocities  as  ice  saturation  increases 
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from  zero  to  about  13  percent.  This  increase  is  much  larger  than  that  occurring  between 
low  and  full  saturation. 


Table  1.  Results  of  experimental  measurements. 


Matrix-Cement 

Dry-Matrix  Porosity 

Ice  Saturation 

Vd  (m/s) 

Vs  (m/s) 

Glass-Ice 

0.4066 

I.O 

4583. 

2492. 

Ottawa  Sand-Ice 

0.3875 

1.0 

4480. 

2730. 

Glass-Ice 

0.4066 

0.131 

3418. 

2100. 

Ottawa  Sand-Ice 

0.3875 

0.135 

3040. 

1900. 

Glass-Ice 

0.4066 

0.0 

0. 

0. 

Ottawa  Sand-Ice 

0.3875 

0.0 

0. 

0. 
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Figure  1.  Velocities  in  uncompacted  frozen  glass  beads  and  Ottawa  sand  versus  ice  saturation. 


EFFECT  OF  CEMENTATION 

The  observed  strong  velocity  increase  at  low  ice  saturations  is  apparently  due  to  the 
cementing  effect  of  frozen  pendular  rings.  Indeed,  the  theoretical  and  experimental  results 
of  Dvorkin  et  al.  (1994)  show  that  even  small  amounts  of  relatively  soft  cement,  if  placed  at 
grain  contacts,  dramatically  increase  the  stiffness  of  an  aggregate.  We  applied  the 
cementation  theory  of  Dvorkin  et  al.  (1994)  to  calculate  the  effective  elastic  moduli  of  ice- 
cemented  glass  beads  and  Ottawa  sand  (the  algorithm  is  given  in  Appendix).  Both 
aggregates  were  modeled  as  assemblies  of  randomly  packed  identical  elastic  spheres  (glass 
and  quartz)  with  elastic  cement  (ice)  at  the  contacts.  The  material  properties  used  in  these 
calculations  are  summarized  in  Table  2. 

All  densities,  as  well  as  P-  and  S-velocities  in  ice,  were  measured  by  the  authors.  P- 
and  S-velocities  in  glass  were  taken  from  Berge  et  al.  (1993).  The  velocities  in  Ottawa 
sand  grains  were  obtained  indirectly,  by  using  Hill's  average  (1952)  to  calculate  the 
effective  moduli  of  the  fully-saturated  sand-ice  mixture,  and  adjusting  those  of  the  grains  to 
match  the  experimental  results.  The  equations  for  calculating  an  effective  modulus  ( Mu;,, ) 
are; 


^Hill  ~  (^V  ^r)  ^  ~  ^o)^Sand  ^O^Ice’ 

1/M,=(1-0o)/Ms,„,  +  (/>o/Mj,„ 
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where  and  are  the  moduli  of  the  sand  grain  material  and  ice,  respectively,  and 
(^0  is  the  porosity  of  a  dry  matrix.  These  equations  can  be  used  to  calculate  bulk  (.K”), 
shear  ( G ),  and  P-wave  ( M )  moduli.  The  latter  is  related  to  P-wave  velocity  d&  M  =  pVj-, 
where  p  is  density. 


Table  2.  Material  properties  of  ice,  glass,  and  sand  grains. 


Substance 

Vs  (m/s) 

Ice 

3840. 

1980. 

900. 

Glass 

5860. 

3480. 

2505. 

Sand  Grains 

5372. 

3517. 

2668. 

The  values  thus  obtained  for  and  in  the  sand-grain  material  are  about  86%  to 

89%  those  in  pure  quartz  -  6050  m/s  and  4090  m/s,  respectively  (Carmichael,  1989). 
Nevertheless,  we  decided  to  use  the  calculated  values  because  the  elastic  moduli  of  sand 
grains  may  be  lower  than  those  of  pure  quartz  due  to  microcracks. 

Our  theoretical  calculations  are  compared  with  the  experimental  values  of  the  P-wave 
and  shear  moduli  in  Figures  2  and  3.  The  moduli  are  plotted  versus  the  porosity  (j)  of  an 
aggregate  —  ^  =  ^o(l  -  5) ,  where  S  is  the  volumetric  ice  saturation.  The  theoretical  curves 
begin  at  porosity  0.2  because  the  cementation  theory  can  be  used  only  for  small  amounts  of 
intergranular  cement.  We  observe  a  very  good  match  between  the  experimental  data  and 
the  theoretical  predictions  for  glass  beads.  In  Ottawa  sand,  the  accuracy  of  the  theoretical 
prediction  is  still  satisfactory  (the  error  is  15  percent)  for  the  P-wave  modulus,  whereas  for 
the  shear  modulus  the  error  reaches  34  percent.  Notice,  that  for  and  V^,  these  errors 

become  7  percent  and  16  percent,  respectively.  This  difference  between  the  predicted  and 
the  measured  values  is  most  likely  due  to  the  angularity  of  the  sand  grains:  the  cemented 
contact  area  between  angular  grains  is  smaller  than  that  between  two  spheres  (for  the  same 
amount  of  cement). 


Figure  2.  P-wave  (left)  and  shear  (right)  moduli  in  ice-cemented  glass  beads.  Circles  -  experiments, 
solid  lines  -  the  cementation  theory,  dashed  line  -  modified  Voigt  average.  The  porosity  is  that  of  the 
frozen  sample,  including  ice  as  part  of  the  matrix;  then,  the  0  porosity  point  is  from  the  fully-saturated 
case,  and  the  0.4  porosity  point  is  from  the  dry  case,  with  the  other  point  representing  partial 
saturation. 
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Figure  3.  P-wave  (left)  and  shear  (right)  moduli  in  ice-cemented  Ottawa  sand.  Circles  -  experiments, 
solid  lines  -  the  cementation  theory,  dashed  line  -  modified  Voigt  average.  Porosity  is  defined  in 
Figure  2. 


The  second  reason  for  the  disparity  may  be  a  more  uniform  distribution  of  ice  on  the 
surface  of  the  Ottawa  sand  grains.  In  this  case  only  part  of  the  ice  --  the  intersection  of  ice 
shells  around  the  grains  -  forms  contact  cement.  We  model  the  latter  situation  by 
representing  Ottawa  sand  as  a  pack  of  identical  spheres  with  ice  uniformly  distributed  on 
their  surfaces.  This  model  gives  a  low-bound  estimate  for  the  effective  elastic  moduli. 
This  lower  bound  is  shown  in  Figure  3  for  the  effective  shear  modulus.  The  span  between 
the  upper  and  lower  bounds  is  fairly  large  for  elastic  moduli.  However,  this  span  becomes 
much  smaller  for  P-  and  S-wave  velocities  (Figure  4).  Therefore,  our  theoretical  bounds 
can  be  used  for  practical  estimates  of  velocities  in  frozen  partially  saturated  sands. 
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Figure  4.  P-  (left)  and  S-wave  (right)  velocities  in  ice-cemented  Ottawa  sand.  Circles  -  experiments, 
solid  lines  -  the  cementation  theory,  upper  and  lower  bounds. 
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CRITICAL  POROSITY  AND  EFFECTIVE  MEDIUM  CALCULATIONS 

Advanced  effective  medium  theories,  such  as  self-consistent  approximation  (e.g., 
Berryman,  1980;  Berge  et  ah,  1993)  have  been  successfully  used  for  modeling  the 
effective  properties  of  composites.  However,  these  theories  do  not  make  explicit  use  of  the 
fact  that  at  zero  confining  pressure,  in  dry  unconsolidated,  uncemented  granular  materials 
elastic  waves  do  not  propagate,  and  elastic  moduli  are  zero.  In  this  case,  the  elastic 
modulus  curves  (in  the  modulus-porosity  plane)  have  to  intercept  the  porosity  axis  at  a 
certain  terminal  (critical)  porosity  ( In  sands,  this  critical  porosity  is  about  0.38  -  the 
value  close  to  the  porosity  of  randomly  packed  identical  spheres.  Generally,  in  effective 
medium  theories  this  intercept  occurs  at  values  higher  than  the  realistic  (j)^  value  for  a  given 
aggregate  (e.g.,  Berge  et  al.,  1993).  Nur  et  al.  (1991)  propose  that  this  critical  porosity 
value  is  a  natural  physical  boundary  for  using  effective  medium  relations  in  sandstones  and 
sands.  They  show  that  one  can  arrive  at  very  accurate  estimates  for  velocities  in  sandstones 
by  using  a  heuristic,  modified  Voigt  average,  formula 

•^Effective  ~  ^Solid 


where  may  be  the  solid  phase's  bulk,  shear,  or  P-wave  modulus.  The  solid  phase 
moduli  have  to  be  calculated  from  those  of  mineral  components  by  using  weighted 
averages,  such  as  Hill's. 

We  now  apply  equation  (2)  to  estimating  the  effective  modulus  of  ice-cemented  glass 
beads  and  Ottawa  sand  by  assuming  that  the  grains  and  the  ice  form  a  single  solid  phase 
whose  effective  moduli  can  be  calculated  from  Hill's  average  —  equation  (1).  The  latter 
assumption  is  supported  by  the  fact  that  Hill's  average  gives  fairly  accurate  estimates  to  the 
effective  moduli  of  frozen  glass  beads  at  full  ice  saturation. 

The  results  of  these  calculations  are  given  in  Figures  2  and  3.  In  many  rocks  the 
modified  Voigt  average  formula  can  accurately  predict  the  moduli  at  small  and  medium 
porosities  (Nur  et  al.,  1991).  We  show  here  that  it  may  fail  to  predict  the  dramatic  increase 
of  these  moduli  that  occurs  as  porosity  slightly  reduces  from  its  initial  (dry)  value.  The 
reason  is  that  this  average  does  not  take  into  account  the  special  arrangement  of  the  ice 
cement  among  grains  -  directly  at  their  contacts.  Practically,  this  arrangement  is  relevant  to 
the  case  of  frozen  partially  saturated  sands. 

Zimmerman  and  King  (1985)  give  an  important  example  of  applying  an  effective 
medium  theory  (Kuster  and  Toksoz,  1974)  to  modeling  velocities  in  fully  saturated  freezing 
sands,  shales,  and  silts.  The  theoretical  scheme  is:  First  the  elastic  moduli  of  the  ice-air 
mixture  (the  ice  is  the  matrix)  are  computed  by  using  the  Kuster-Toksoz  formula  for 
spherical  inclusions.  Then  the  effective  moduli  of  the  frozen  aggregate  are  computed  by 
assuming  that  the  ice-air  mixture  is  the  matrix,  and  the  grains  are  spherical  inclusions.  By 
comparing  these  predictions  with  experimental  data,  Zimmerman  and  King  show  that  this 
theoretical  approach  gives  accurate  estimates  to  velocities  in  permafrost.  We  applied  the 
recommended  theoretical  scheme  to  calculating  P-wave  and  shear  moduli  for  partially 
saturated  frozen  glass  beads  (Figure  5). 

The  Kuster-Toksoz  scheme  gives  adequate  estimates  for  the  moduli  at  full  saturation. 
But  again,  it  fails  to  mimic  the  high-porosity  elastic  properties  of  partially  saturated  frozen 
granular  materials.  A  likely  reason  is  that  in  fully  saturated  rocks  ice  forms  first  in  the 
larger  pore  spaces,  and  then  in  progressively  smaller  ones  (Williams  and  Smith,  1989). 
However,  in  partially  saturated  hydrophilic  materials  where  water  forms  pendular  rings,  ice 
forms  at  grain  contacts  and  acts  as  strong  cement.  In  this  case  the  cementation  theory  of 
Dvorkin  et  al.  (1994)  is  apparently  more  appropriate. 


56 


Porosity 

Figure  5.  Kuster-Toksoz-Zimmerman  estimates  for  elastic  P-wave  {M)  and  shear  ( G )  moduli  in 
frozen  glass  beads.  Circles  --  experimental  points  for  M,  triangles  -  experimental  points  for  G. 

CONCLUSIONS 

In  hydrophilic  partially  saturated  frozen  sands,  ice  acts  as  intergranular  cement  by 
forming  pendular  rings  around  grain  contacts.  This  arrangement  leads  to  strong  P-  and  S- 
wave  velocity  increase  between  dry  and  partially  saturated  frozen  samples  at  low 
saturations.  This  velocity  increase  is  much  larger  than  that  occurring  between  low  and  full 
saturation. 

The  experimental  results  can  be  accurately  modeled  by  the  theory  of  cementation  that 
explicitly  takes  into  account  the  special  position  of  intergranular  cement  at  grain  contacts. 
The  effective  medium  theories,  while  giving  accurate  estimates  for  velocities  at  high  ice 
saturation,  fail  to  predict  the  observed  dramatic  velocity  increase  occurring  as  small 
amounts  of  ice  are  added  to  grain  contacts. 

By  assuming  two  geometrical  schemes  of  ice  distribution  in  the  pore  space  (one  in 
which  ice  is  concentrated  at  the  grain  contacts,  and  another  in  which  ice  is  uniformly 
distributed  on  the  grain  surface)  the  theory  of  cementation  provides  fairly  narrow  upper- 
and  lower-bound  estimates  for  P-  and  S-wave  velocities  in  partially  saturated  frozen  sands. 
These  bounds  can  be  used  in  practical  estimates  of  the  elastic  properties  of  permafrost. 
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2.6.  ELASTICITY  OF  HiGH-POROSITY  SANDSTONES:  THEORY 
FOR  Two  NORTH  SEA  DATASETS 


ABSTRACT 

We  have  analyzed  two  laboratory  datasets  obtained  on  high-porosity  rock  samples  from 
the  North  Sea.  The  velocities  observed  are  unusual  in  that  they  seem  to  disagree  with  some 
simple  models  based  on  porosity.  On  the  other  hand,  the  rocks  are  unusually  poorly- 
cemented  (for  laboratory  studies,  at  least),  and  we  investigate  the  likelihood  that  this  is  the 
cause  of  the  disagreement.  One  set  of  rocks,  from  the  Oseberg  field,  is  made  of  slightly 
cemented  quartz  sands.  We  find  that  we  can  model  their  dry-rock  velocities  using  a 
cementation  theory  where  the  grains  mechanically  interact  through  cement  at  the  grain 
boundaries.  This  model  does  not  allow  for  pressure-dependence.  The  other  set  of  rocks, 
from  the  Troll  field,  is  almost  completely  uncemented.  The  grains  are  held  together  by  the 
applied  confining  pressure.  In  this  case,  a  lower  bound  for  the  velocities  can  be  found  by 
using  the  Hertz-Mindlin  contact  theory  (interaction  of  uncemented  spheres)  to  predict 
velocities  at  a  critical  porosity,  combined  with  the  modified  Hashin-Strikman  lower  bound 
for  other  porosities.  This  model,  which  allows  for  pressure-dependence,  also  predicts 
fairly  large  Poisson's  ratios  for  saturated  rocks,  such  as  those  observed  in  the 
measurements.  The  usefulness  of  these  theories  may  be  in  estimating  the  nature  of  cement 
in  rocks  from  measurements  such  as  sonic  logs.  The  theories  could  help  indicate  sand 
strength  in  poorly-consolidated  formations  and  predict  the  likelihood  of  sand  production. 
Both  theoretical  methods  have  analytical  expressions  and  are  ready  for  practical  use. 

INTRODUCTION 

Exploring  seismic-velocity-to-porosity  transformations  in  various  lithologies  has  been 
one  of  the  important  research  areas  in  rock  physics.  The  main  application  of  this 
knowledge  is  to  predict  porosity  from  sonic  logs  and  from  seismic. 

Often  the  elastic  moduli  of  rock  are  used  instead  of  velocities.  Among  them  are  the 
compressional-wave  modulus  (or  the  M -modulus)  and  the  shear  modulus  ( G).  They  can 
be  expressed  through  V^,  and  density  p  as 

M=pv;, 

G  =  pV,\ 

It  has  been  established  (Nur  et  al.,  1991)  that  in  consolidated  sandstones  elastic  moduli 
depend  approximately  linearly  on  porosity.  The  straight  line  connects  two  end  members; 
one  has  zero  porosity  and  a  modulus  close  to  that  of  the  solid  phase,  another  has  "critical" 
porosity  and  (for  dry  rock)  a  modulus  close  to  zero.  The  critical  porosity  value  for 
sandstones  is  between  0.36  and  0.4.  Two  groups  of  high-porosity  sandstones  from  the 
Oseberg  and  the  Troll  fields  defy  this  simple  rule  (Figure  1).  They  represent  two  different 
trends  where  the  modulus  changes  with  changing  porosity  not  as  dramatically  as  it  does  in 
consolidated  sandstones.  Our  goal  is  to  understand  the  physical  principles  that  govern 
velocity-porosity  relations  in  unconsolidated  sandstones  and  to  develop  a  predictive  theory. 

The  two  laboratory  datasets  have  been  obtained  on  high-porosity  rock  samples  from  the 
North  Sea.  The  first  set  is  from  the  Oseberg  field  (Strandenes,  1991).  Some  of  these 
samples  are  shales  and/or  have  strong  intrinsic  anisotropy.  We  select  a  subset  that  includes 
structurally  and  acoustically  isotropic  quartz  sands  with  porosities  between  0.13  and  0.32. 
Quartz  volumetric  content  varies  from  0.6  to  almost  1.0.  The  thin  section  images  reveal 
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slight  quartz  cementation  among  the  grains.  Two  other  dominant  minerals  here  are  mica 
(up  to  0.15  volume  fraction)  and  clay  (up  to  0.25  volume  fraction). 


Figure  1.  The  compressional  wave  modulus  (Af -modulus)  and  the  shear  modulus  versus  porosity  for 
three  groups  of  dry  rocks  at  30  MPa  confining  pressure:  consolidated  sandstones  without  or  with  only 
small  amounts  of  clay  (Han,  1986),  Oseberg  sands  (Strandenes,  1991),  and  Troll  sands  (Blangy,  1992). 

The  second  set  represents  unconsolidated  sands  from  the  Troll  field  with  porosities 
between  0.22  and  0.38  (Blangy,  1992;  and  Blangy  et  al.,  1993).  The  thin  section  iniages 
show  that  intergranular  cementation  here  is  practically  absent.  Quartz  is  the  dominant 
mineral  in  these  sands  with  volume  fraction  varying  between  0.5  and  0.8.  Two  other 
major  mineral  components  are  feldspar  and  mica. 

Based  on  the  apparent  microstructure  of  the  two  rock  groups  we  apply  two  theories  to 
describe  them.  The  first  one  is  the  cementation  theory  (Dvorkin  et  al.,  1991  and  1994). 
This  theory  describes  the  mechanical  interaction  between  grains  that  are  bound  together  by 
intergranular  cement  without  reliance  on  normal  loading  stresses  to  maintain  contact.  An 
important  implication  of  this  theory  is  that  even  very  soft  cement  is  load-bearing  and  thus 
very  small  amounts  of  cementation  can  significantly  increase  the  stiffness  of  a  granular 
composite.  We  use  this  theory  to  describe  the  Oseberg  samples. 

The  second  theory  is  based  on  the  contact  model  of  Hertz  and  Mindlin  (Mindlin,  1949) 
that  describes  the  mechanical  interaction  between  two  smooth  uncemented  elastic  spheres 
based  on  friction  and  loading  by  normal  stresses.  We  use  this  model  to  calculate  the  elastic 
moduli  of  sand  at  critical  porosity.  To  extend  it  for  different  porosities,  we  use  the 
modified  Hashin-Strikman  lower  bound  where  one  end  member  has  zero  porosity  and  the 
modulus  of  the  solid  phase  and  another  end  member  has  critical  porosity  and  a  pressure- 
dependent  modulus  given  by  the  Hertz-Mindlin  theory.  This  Hertz-Mindlin-Hashin- 
Strikman  (HMHS)  model  allows  us  to  account  for  the  noticeable  pressure  dependence 
normally  observed  in  sands.  We  apply  it  to  describe  the  Troll  samples. 

CONTACT  THEORIES 


Cementation 


The  cementation  theory  by  Dvorkin  et  al.  (1991  and  1994)  gives  the  normal  and  shear 
stiffnesses  of  a  two-grain  combination  with  elastic  cement  at  the  contact.  Stiffnesses  thus 
obtained  can  be  used  in  various  effective  medium  approximations.  An  appropriate 
approximation  for  sand  is  a  random  pack  of  identical  spherical  grains.  The  original  theory 
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requires  some  computational  effort  to  obtain  the  stiffnesses.  By  doing  these  computations 
in  a  broad  range  of  the  elastic  constants  of  the  grains  and  those  of  the  cement,  and  for 
varying  amounts  of  cement  we  have  arrived  at  easily  usable  correlations  as  described  below 
and  in  Appendix. 

We  assume  that  the  starting  framework  of  cemented  sand  is  a  dense  random  pack  of 
identical  spherical  grains  with  porosity  (j)^  ~  0.36  and  the  average  number  of  contacts  per 
grain  n~9.  Adding  cement  to  the  grains  acts  to  reduce  porosity  and  to  increase  the 
effective  elastic  moduli  of  the  aggregate.  Then  these  effective  dry-rock  bulk  (K^^)  and 

shear  ( moduli  are  (Dvorkin  et  al.,  1994) 

— —n(l  — (I>q)McS^,  =—K^^ +—n(l- (j>Q)G^S^, 

K=pcVic^  G.=pcyl-, 

where  and  G^  are  the  compressional-wave  and  the  shear  moduli  of  the  cement, 
respectively;  is  its  density;  and  and  are  its  P-  and  5-wave  velocities. 

Parameters  5„  and  5^  are  proportional  to  the  normal  and  shear  stiffness,  respectively,  of  a 
cemented  two-grain  combination.  They  depend  on  the  amount  of  the  contact  cement  and  on 
the  properties  of  the  cement  and  the  grains  (see  Appendix). 

The  amount  of  the  contact  cement  can  be  expressed  through  the  ratio  a  of  the  radius  of 
the  cement  layer  a  to  the  grain  radius  R : 


a  =  a  /  R. 

The  radius  of  the  contact  cement  layer  a  is  not  necessarily  directly  related  to  the  total 
amount  of  eement:  part  of  the  cement  may  be  deposited  away  from  the  intergranular 
contacts  (Figure  2a).  However  by  assuming  that  porosity  reduction  in  sands  is  due  to 
cementation  only,  and  by  adopting  certain  schemes  of  cement  deposition  we  can  relate 
parameter  a  to  the  current  porosity  of  cemented  sand  (p.  For  example,  we  can  use 
Scheme  1  where  all  cement  is  deposited  at  grain  contacts  (Figure  2b): 
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or  we  can  use  Scheme  2  where  cement  is  evenly  deposited  on  the  grain  surface  (Figure  2c): 
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(3) 


In  formulas  (2)  and  (3)  5  is  the  cement  saturation  of  the  pore  space.  It  is  the  fraction  of  the 
pore  space  (of  the  uncemented  sand)  oceupied  by  cement  (in  the  cemented  sand). 

It  follows  from  the  cementation  theory  that  even  soft  contact  cement  is  load-bearing. 
Thus  even  small  amounts  of  cement,  if  deposited  precisely  at  grain  contacts,  act  to 
dramatically  increase  the  stiffness  of  a  particulate  aggregate.  This  theoretical  conclusion  is 
well  supported  by  experimental  facts  (Figure  3). 
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Figure  2.  Cemented  grains,  a.  Contact  and  non-contact  cement,  b.  Scheme  1  of  cement  deposition, 
c.  Scheme  2  of  cement  deposition. 
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Figure  3.  Examples  of  the  theoretical  and  experimental  values  of  velocities  in  cemented  granular 
materials.  Circles  are  from  experiments,  solid  lines  are  theoretical  predictions.  From  left  to  right: 
epoxy-cemented  glass  beads  (data  from  Yin,  1993);  ice-cemented  Ottawa  sand  (data  from  Jacoby  et  al., 
1995);  and  sintered  glass  beads  (data  from  Berge  et  al.,  1993). 

Hprt7-Minfl1in  Theory  with  the  Modified  Hashin-Strikman  Lower  Bound 

The  contact  Hertz-Mindlin  theory  (Mindlin,  1949)  gives  the  following  expressions  for 
the  effective  bulk  ( and  shear  ( )  moduli  of  a  dry  dense  random  pack  of  identical 
spherical  grains  subject  to  a  hydrostatic  pressure  P : 


Ghm  = 


5-Av  3n\\-<^,fG^ 
5(2 -v)^  2;r'(l-v)' 


(4) 


where  v  is  the  grain  Poisson's  ratio  and  G  is  the  grain  shear  modulus. 

Formula  (4)  describes  the  effective  elastic  properties  of  sand  at  critical  porosity  (p^.  In 
order  to  find  the  effective  moduli  and  at  a  different  porosity  (j)  we  propose  a 
heuristic  modified  Hashin-Strikman  lower  bound  based  on  the  original  Hashin-Strikman 
lower  bound  (1963): 
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where  ^  is  the  grain  bulk  modulus. 

This  model  connects  two  end  members:  one  has  zero  porosity  and  the  modulus  of  the 
solid  phase  and  another  has  critical  porosity  and  a  pressure-dependent  modulus  as  given  by 
the  Hertz-Mindlin  theory  (Figure  4).  It  is  apparent  that  this  contact  theory  allows  one  to 
describe  the  noticeable  pressure  dependence  normally  observed  in  sands.  The  Hashin- 
Strikman  low  bound  is  appropriate  to  describe  materials  close  to  suspensions.  Blangy  et 
al.  (1993)  show  that  such  models  give  accurate  estimates  for  velocities  in  unconsolidated 
sands.  The  applicability  of  contact  theories  to  describing  granular  rocks  has  been  explored 
by  many  authors.  A  detailed  investigations  was  done  by  Murphy  (1982). 


Figure  4.  The  modified  Hashin-Strikman  lower  bound  curves  at  varying  confining  pressure  for  pure 
quartz  sand. 

Notice  that  the  dry-rock  Poisson's  ratios,  as  predicted  by  formula  (4),  do  not  vary  with 
pressure.  At  the  same  time,  porosity  correction  given  by  formula  (5)  makes  the  Poisson's 
ratio  slightly  pressure-sensitive. 

APPLYING  THE  THEORY 


Oseberg 

The  Oseberg  samples  have  noticeable  mineralogical  diversity  with  quartz  as  the  main 
component.  Based  on  this  last  fact  we  model  them  as  pure  quartz  sand  bound  by  quartz  or 
clay  cement.  We  assume  that  cement  is  evenly  deposited  on  the  grain  surface,  as  in 
Scheme  2,  and  calculate  the  dry-rock  velocities  from  formulas  (1)  and  (3).  The  bulk 
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modulus,  shear  modulus,  and  density  of  quartz  are  chosen  as  38  GPa,  44  GPa,  and  2.65 
g/cm^,  respectively  (Carmichael,  1990).  Our  estimates  based  on  the  cementation  theory 
capture  dry-rock  velocity  values  reasonably  well  in  the  relatively  clay-free  Oseberg  samples 
at  confining  pressure  ranging  from  20  MPa  to  40  MPa  (Figure  5).  The  theory  explains  the 
small  variation  of  velocities  with  changing  porosity. 


Porosity  Porosity  Porosity 

Figure  5.  P-wave  velocities  (filled  symbols)  and  S-wave  velocities  (open  symbols)  in  the  dry  Oseberg 
samples  at  confining  pressures  40,  30,  and  20  MPa.  Circles  are  for  quartz-cemented  rock,  triangles  are 
for  clay-cemented  rock.  The  solid  lines  are  our  theoretical  estimates  that  are  pressure-independent.  In 
these  estimates  cement  evenly  envelopes  the  grains  (Scheme  2).  Cement  is  either  quartz  or  clay. 

We  also  calculate  dry-rock  velocities  in  the  Oseberg  samples  by  assuming  that  cement  is 
clay.  The  bulk  modulus,  shear  modulus,  and  density  of  clay  are  chosen  as  21  GPa,  7 
GPa,  and  2.58  g/cm^,  respectively.  The  velocity  values  were  derived  by  extrapolating  the 
experimental  linear  velocity-porosity-clay  relations  obtained  by  Tosaya  (1982)  for  clay¬ 
bearing  Gulf  sandstones.  The  density  was  calculated  from  the  measured  density  of 
samples  with  high  clay  content,  and  known  porosity  and  mineralogy.  The  theoretical  line 
is  close  to  the  data  points  for  the  Oseberg  samples  with  10-25  volume  percent  of  clay.  It 
is  interesting  that  reducing  the  cement's  stiffness  does  not  significantly  reduce  the  stiffness 
of  the  cemented  aggregate  (Figure  5). 

It  is  important  to  mention  that  the  cementation  theory  currently  does  not  account  for 
velocity  changes  with  changing  confining  pressure.  The  effective  reservoir  pressure  in  the 
Oseberg  field  is  about  23  MPa.  It  is  crucial  that  the  theory  works  well  at  this  confining 
pressure.  In  fact,  the  theoretical  bounds  slightly  overestimate  the  observed  values  at  20 
MPa.  We  speculate  that  one  reason  for  this  effect  may  be  core  damage  during  pressure 
release  and  thus  a  permanent  velocity  reduction.  Holt  et  al.  (1994)  show  that  if  a  rock 
sample  is  unloaded  and  then  reloaded  to  the  initial  pressure  level,  the  velocities  do  not 
recover  their  initial  values.  Thus  we  conclude  that  the  cementation  theory  has  good 
predictive  power  for  velocities  in  slightly  cemented  sands.  It  allows  one  to  determine  the 
cementing  mineral  (quartz  or  clay  in  this  example). 

The  cementation  theory  lacks  predictive  power  when  applied  to  calculate  Poisson's 
ratios  in  dry,  as  well  as  in  saturated  rock  (Figure  6). 
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Figure  6.  Poisson's  ratios  in  the  Oseberg  samples  at  the  confining  pressures  of  20  MPa  (open  circles) 
and  40  MPa  (filled  circles).  Solid  lines  are  our  theoretical  predictions  (Scheme  2  of  cement  deposition). 
The  water- saturated-rock  values  are  calculated  from  the  dry-rock  values  using  Gassmann's  (1951) 
formula. 


Troll 


Some  of  the  Troll  samples  have  large  amounts  of  feldspar  and  mica.  It  would  be 
logical  to  apply  the  HMHS  scheme  to  grains  made  out  of  a  "composite"  solid  phase  that  is  a 
mixture  of  the  dominant  minerals  present  in  the  samples.  We  find  that  the  theoretical 
velocities  calculated  using  such  mineral  mixtures  are  very  close  to  those  calculated  using 
pure  quartz  grains.  The  theoretical  results  for  dry  rock  generally  serve  as  an  accurate  lower 
bound  for  the  Troll  experimental  values  (Figure  7). 
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Figure  7.  Dry-rock  velocities  in  the  Troll  samples  at  the  confining  pressures  of  5  MPa,  15  MPa,  and 
30  MPa.  Filled  symbols  are  for  P-wave  velocities,  open  symbols  are  for  S-wave  velocities.  Solid 
lines  are  our  theoretical  estimates.  Triangles  are  for  Ottawa  sand  samples  measured  by  Han  (1986),  Yin 
(1993),  and  Domenico  (1984). 

The  data  for  dry  Ottawa  sand  at  30  MPa  confining  pressure  is  plotted  in  these  graphs 
for  reference.  There  is  a  large  disparity  not  only  between  the  theoretical  and  the 
experimental  values  for  Ottawa  sand  but  also  among  the  experimental  values  for  the  three 
samples.  One  reason  may  be  a  variation  in  compaction  and  sorting  among  these  samples. 
As  a  result,  the  average  number  of  contacts  per  grain  may  not  always  be  9  and  may  vary 
among  the  three  samples  of  Ottawa  sand. 
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Low-frequency  water-saturated-rock  velocities  are  computed  from  the  dry-rock  data 
using  Gassmann's  (1951)  formula  (Figure  8).  The  relative  errors  of  the  theoretical 
estimates  are  reduced  as  compared  to  the  dry-rock  results.  There  are  two  apparent  reasons 
for  this  improvement.  First,  pore  fluid  acts  to  reduce  the  dry-rock  velocity  variation  among 
low- velocity  samples.  Second,  it  acts  to  strongly  increase  P-wave  velocities  thus  reducing 
relative  errors. 
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Figure  8.  Low-frequency  saturated-rock  velocities  (from  Gassmann's  formula)  in  the  Troll  samples  at 
the  confining  pressures  of  5  MPa,  15  MPa,  and  30  MPa.  Filled  symbols  are  for  P-wave  velocities, 
open  symbols  are  for  S-wave  velocities.  Solid  lines  are  our  theoretical  estimates. 


Figure  9.  Dry-rock  and  low-frequency  saturated-rock  Poisson's  ratios  in  the  Troll  samples  at  the 
confining  pressures  of  5  MPa,  15  MPa,  and  30  MPa.  Solid  lines  are  our  theoretical  estimates. 

The  last  reason  is  also  responsible  for  the  accurate  predictions  the  HMHS  model 
provides  for  saturated-rock  Poisson's  ratios  (Figure  9).  These  values  are  high  and  match 
those  generally  observed  in  saturated  loose  sediments.  The  physical  meaning  of  these  high 
values  is  clear;  saturated  loose  sediments,  especially  with  no  confining  pressure,  are  close 
to  suspensions  where  Poisson's  ratios  have  their  maximum  value  of  0.5. 

In  the  HMHS  model,  the  saturated-rock  Poisson's  ratio  increases  with  increasing 
porosity.  This  trend  qualitatively  agrees  with  that  in  the  data  and  has  a  clear  physical 
meaning  —  as  porosity  increases,  the  rocks  approach  the  state  of  suspension  and  the 
Poisson's  ratio  approaches  0.5.  At  the  same  time,  it  is  obvious  that  the  Poisson's  ratios  of 
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saturated  rock  obtained  from  formula  (4),  which  does  not  account  for  porosity  variation, 
will  do  as  good  quantitative  predictions  as  those  obtained  from  formula  (5)  where  porosity 
is  accounted  for. 

An  important  feature  of  the  HMHS  model  is  that  while  the  dry-rock  Poisson's  ratios  are 
only  slightly  pressure-dependent,  the  saturated-rock  Poisson's  ratios  have  a  noticeable 
pressure  dependence  (Figure  9). 

We  conclude  that  the  HMHS  model  can  be  used  to  estimate  and  bound  velocities  in 
loose  sands  in  the  Troll  field.  It  can  accurately  predict  saturated-rock  Poisson's  ratios  but 
should  not  be  used  to  estimate  dry-rock  Poisson's  ratios  at  high  pressure. 

DISCUSSION:  PRACTICAL  VALUE  OF  THE  THEORIES 

Oseberg 

We  learn  from  the  above  theoretical  speculations  that  in  the  Oseberg  rocks  a  significant 
part  of  cement  is  deposited  at  grain  contacts.  This  cementation  pattern  ensures  strong 
velocity  increase  with  only  small  porosity  decrease  and  is  likely  to  be  responsible  for  the 
high  velocities  measured  on  these  high-porosity  samples.  We  speculate  that  this  result  can 
be  generalized  to  give  a  rule:  If  high-porosity  rock  has  velocity  values  close  to  those  given 
by  the  cementation  theory  then  it  is  cemented  at  grain  contacts.  The  amount  of  the  contact 
cement  can  be  estimated  by  varying  parameter  a  in  the  cementation  theory  (see  Appendix) 
to  match  the  data  point.  Thus  velocity  measurements  reveal  the  internal  rock  structure. 

The  latter  information  is  important  in  assessing,  for  example,  a  rock's  strength.  Yin 
and  Dvorkin  (1994)  show  that  even  very  small  amounts  of  contact  cement  (even  if  it  is  as 
soft  as,  e.g.,  clay)  prevent  the  grains  from  breaking  at  high  confining  pressure  (Figure  10). 
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Figure  10.  Porosity  versus  hydrostatic  confining  pressure  in  water-saturated  randomly  packed  glass 
beads  that  were  a)  uncemented  and  b)  cemented  by  epoxy  at  their  contacts  (Yin  and  Dvorkin,  1994). 
Porosity  was  measured  by  the  volume  of  expelled  fluid.  In  the  uncemented  case,  a  sharp  porosity 
decrease  is  observed  at  about  50  MPa.  The  decrease  is  associated  with  the  crushing  of  grains.  The 
cemented  grains  (the  volume  of  the  epoxy  accounted  for  only  10  percent  of  the  pore  space)  did  not 
crush.  The  photos  showed  that  in  the  latter  case,  the  grains  stayed  intact  with  the  failure  being 
localized  within  the  epoxy. 

We  conclude  that  if  high-porosity  sandstones  have  velocities  in  the  range  predicted  by 
the  cementation  theory  they  are  mechanically  stable  and  sanding  is  unlikely. 

The  practical  significance  of  the  cementation  theory  is  not  only  in  its  predictive  power  - 
it  is  likely  that  a  simple  linear  velocity-porosity  correlation  obtained  from  the  data  will  serve 
even  better  for  the  purpose  of  estimating  porosity  from  seismic.  The  practical  value  is  also 
in  understanding  the  physical  reasons  for  the  measured  velocities.  This  understanding  may 
help  assess  a  rock's  mechanical  stability  and  lead  even  further  to  estimating,  for  example, 
such  a  crucial  parameter  as  permeability,  once  we  are  aware  of  the  rock's  internal  structure. 
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Troll 


Velocity  values  in  the  Troll  samples  can  be  adequately  approximated  by  the  HMHS 
model.  This  fact  implies  that  the  rock  is  held  together  by  confining  pressure  only.  The 
grains  that  are  uncemented  at  high  porosity  stay  uncemented  as  porosity  decreases.  Thus 
the  additional  material  responsible  for  porosity  reduction  is  deposited  in  the  large  pores, 
away  from  compliant  contacts.  This  in  turn  implies  that  the  rock  is  mechanically  unstable 
and  sanding  is  likely  once  uniform  confining  pressure  is  released.  An  important  practical 
result  is  that  the  predicted  very  high  Poisson's  ratios  in  saturated  rocks  closely  match  the 
observed  values.  Thus  by  combining  seismic  and  depositional  history  one  can  distinguish 
this  unconsohdated  and  mechanically  unstable  rock  from  other  types. 

CONCLUSION 

The  cementation  theory  explains  well  why  the  high  velocity  values  are  encountered  in 
the  high-porosity  Oseberg  samples.  This  fact  means  that  the  grains  mechanically  interact 
through  quartz  cement  that  is  deposited  to  strengthen  the  most  compliant  part  of  the  pore 
space  -  grain  contacts.  The  main  practical  conclusion  is  that  if  at  high  porosities  velocity 
values  are  close  to  those  provided  by  the  cementation  theory  then  the  grain  contacts  are 
cemented.  This  in  turn  means  that  the  rock  is  mechanically  stable  and  sanding  is  not  likely. 

Velocity  values  in  the  Troll  samples  can  be  estimated  from  a  combination  of  the  Hertz- 
Mindlin  contact  theory  and  the  modified  Hashin-Strikman  lower  bound,  the  latter  being 
close  to  an  isostress  model  for  suspensions.  Thus  in  this  rock  contact  cement  is  almost 
absent  and  grains  are  held  together  primarily  by  confining  pressure.  The  model  accurately 
predicts  the  high  Poisson's  ratios  generally  observed  in  saturated  loose  sediments. 

The  main  practical  conclusion  is  that  if  at  high  porosities  velocity  values  are  close  to 
those  provided  by  the  HMHS  model  then  there  is  no  contact  cement  and  the  rock  is  held 
together  by  confining  pressure  only.  Therefore,  the  rock  is  mechanically  unstable  and 
sanding  is  likely  to  occur  once  uniform  confining  pressure  is  released. 

By  analyzing  the  data  and  finding  the  appropriate  theories  to  explain  the  observed 
trends  we  gain  understanding  of  the  internal  rock  structure  and  of  stress  transfer  patterns. 
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APPENDIX:  EXPRESSIONS  FOR  STIFFNESSES 


Parameters  and  in  formulas  (1)  are: 


5„=A„(A„)a"  +  B„(A„)a  +  C„(A„), 

A„(A„)  =  -0.024153  •  A„-'•'^^  fi„(A„)  =  0.20405  • 
C„(A„)  =  0.00024649  • 


5,  =  A,(A,,  v)a^  +  v)a  +  Q(A„  v), 

A,(A„  V)  = -10-' •  (2.26 v' -b  2.07  v-b  2.3)  • 

i5,(A„  V)  =  (0.0573  V'  -b  0.0937  V-b  0.202)  •  a/  '’274vU0.0529v-0.8765^ 

C,(A,,  V)  =  10“"  •(9.654v'  -b4.945v-b3.1)  • 


2G,(l-v)(l-v,)  _  G, 

- - 9  ,  VA/  - 

kG  1-2v  kG  R 


where  G  and  v  are  the  shear  modulus  and  the  Poisson's  ratio  of  the  grains,  respectively; 
G^  and  are  the  shear  modulus  and  the  Poisson's  ratio  of  the  cement,  respectively;  a  is 

the  radius  of  the  contact  cement  layer  (Figure  2b,  c);  and  R  is  the  grain  radius. 

These  formulas  are  statistical  approximations  of  the  rigorous  cementation  theory 
solutions  (Dvorkin  et  al.,  1994).  The  error  does  not  exceed  one  percent. 
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2.7.  LARGE  STRAINS  IN  CEMENTED  GRANULAR  AGGREGATES: 

ELASTIC-PLASTIC  CEMENT 


1.0.  ABSTRACT 

We  describe  large-strain  behavior  of  cemented  geomaterials  by  modeling  the 
deformation  of  a  random  pack  of  identical  cemented  spheres.  In  this  model  we  assume  that 
the  grains  are  elastic  but  that  the  intergranular  cement  becomes  partly  plastic  as  local 
stresses  meet  a  plasticity  condition.  This  plasticity  condition  for  a  thin  elastic-plastic 
cement  layer  is  derived  based  on  the  von  Mises  criterion.  Next  we  solve  the  problems  of 
large-strain  deformation  of  two  cemented  spheres  in  compression  and  in  shear.  This 
solution  allows  for  calculation  of  the  normal  and  shear  strain-dependent  stiffnesses  of  two 
cemented  grains.  Finally,  we  derive  effective  stress-strain  laws  for  an  aggregate  of 
cemented  spheres  where  grain-to-grain  contact  stiffnesses  are  strain-dependent.  These 
theoretical  stress-strain  laws  are  close  to  the  relations  typically  observed  in  clay-cemented 
sands.  An  important  result  is  that  an  initially  isotropic  aggregate  may  become  anisotropic  if 
the  stress  field  is  non-hydrostatic.  This  stress-induced  anisotropy  may  lead  to  an  up  to  10- 
percent  error  in  estimating,  for  example,  an  apparent  isotropic  bulk  modulus  from  a 
uniaxial-strain  experiment. 

2.0.  INTRODUCTION 


2.1.  Background 

The  large-strain  mechanical  properties  of  natural  and  artificial  geomaterials  are 
important  in  many  geomechanical  and  engineering  applications.  Many  granular 
geomaterials  are  cemented  with,  for  example,  clay  (sand-clay  mixtures)  or  asphalt 
(asphaltic  concrete).  Dvorkin  et  al.  (1994),  and  Yin  and  Dvorkin  (1994)  show  that 
intergranular  cementation  may  strongly  affect  the  elastic  properties  of  particulate  materials 
as  well  as  their  strength.  It  is  naturd  to  conclude  that  cementation  affects  the  large-strain 
microscopic  and  macroscopic  behavior  of  such  materials  as  well. 

This  effect  is  apparent  in  the  stress-strain  curves  from  uniaxial-stress  compression  tests 
on  two  tight  gas  sandstone  samples  (Jizba,  1991).  In  one  sample  quartz  grains  are 
cemented  by  a  clay  matrix  that  occupies  63  percent  of  the  sample’s  volume.  In  the  second 
sample  the  volumetric  content  of  the  clay  matrix  is  only  1  percent.  The  stress-strain 
behavior  of  these  two  samples  is  distinctively  different  (Figure  1).  The  low-clay-content 
sample  deforms  almost  linearly  up  to  about  2  percent  strain  and  then  exhibits  brittle  failure. 
The  clay-cemented  sample  shows  ductile  (plastic)  behavior  with  the  stress-strain  curves 
gently  approaching  failure. 

Another  argument  that  supports  the  importance  of  intergranular  cementation  is  the 
difference  between  static  and  dynamic  moduli  in  rocks.  Static  moduli  are  those  obtained 
from  appropriate  stress-strain  curves  (e.g.,  static  Young  modulus  is  the  derivative  of  the 
stress-strain  curve  in  a  uniaxial-stress  experiment).  Dynamic  moduli  are  those  calculated 
from  P-  and  S-wave  velocities  in  rocks.  For  example,  dynamic  bulk  and  shear  moduli  are; 

•^Dynamic  ~  P(^ p  ^Dynamic  ~  P^s  ’ 

where  p  is  density,  and  and  are  P-  and  S-wave  velocities,  respectively. 

Typically,  dynamic  moduli  in  dry  rocks  are  larger  than  the  corresponding  static  moduli. 
One  reason  is  the  effect  of  strain.  A  sound  wave  propagating  through  rock  results  in  minor 
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stress  excursions  around  a  given  point  on  the  stress-strain  curve.  The  slope  of  such  an 
excursion  curve  is  larger  than  that  of  the  major  stress-strain  curve  in  loading  and  is  very 
close  to  the  corresponding  dynamic  modulus  (Hilbert  et  al.,  1994).  This  effect  is  due  to  the 
plastic  behavior  of  rocks  at  high  strain. 


Quartz 
Grain  i 

V 


Clay 


Cement 


Figure  1 .  Stress-strain  curves  from  uniaxial-stress  compression  tests  on  sandstone  samples.  In  each 
test  the  radial  stress  is  constant  (labels  on  the  curves).  Two  upper  curves  correspond  to  a  low-clay- 
content  sample  which  failed  in  a  brittle  manner.  Two  lower  curves  correspond  to  a  sample  with  63 
percent  volumetric  clay  content.  The  sample  exhibited  ductile  deformation  prior  to  failure.  It 
eventually  failed  by  shear  localization. 

If  confining  pressure  is  small,  this  plasticity  is  partly  due  to  grain  compaction, 
redistribution,  and  rotation.  However  at  high  confining  pressure  plasticity  is  mostly  due  to 
the  plastic  nature  of  intergranular  cement.  This  effect  can  be  clearly  seen  in  Figure  2  where 
the  ratios  of  dynamic  to  static  bulk  moduli  (the  latter  from  hydrostatic  compaction 
experiments)  are  plotted  versus  hydrostatic  pressure  (Jizba,  1991). 

In  two  samples  with  almost  no  clay  this  ratio  quickly  approaches  unity  as  pressure 
increases.  We  conclude  that  once  grain  redistribution  and  rotation  is  finished,  rock  behaves 
elastically.  However  in  clay-cemented  samples  this  ratio  stays  above  unity  which  means 
that  the  clay  cement  strongly  contributes  to  the  plasticity  of  such  rocks. 

2.2.  Model  and  Assumptions 

The  goal  of  this  paper  is  to  theoretically  describe  the  apparently  important  effect  of 
cement's  plasticity  on  the  macroscopic  behavior  of  an  aggregate.  Our  principal  approach  is 
to  model  large  strains  in  a  cemented  aggregate  by  assuming  that  cement  is  elastic-plastic, 
that  is  part  of  a  cement  layer  becomes  plastic  as  stress  concentration  increases. 

Contact  stresses  (Figure  3  a)  between  cemented  grains  are  principally  different  from 
those  calculated  using  traditional  Hertzian  models  (e.g..  Hertz,  1882;  Mindlin,  1949; 
Walton,  1978;  Digby,  1981).  Dvorkin  et  al.  (1994)  show  that  if  cement  is  soft  as 
compared  to  grains,  contact  stress  concentrations  occur  at  the  center  of  a  cemented  region 
(Figures  3b  and  3c).  It  reasonable  to  assume  that  as  the  contact  stresses  at  the  center  of  a 
cemented  region  meet  a  plasticity  condition,  the  central  part  of  the  cement  becomes  plastic 
and  is  surrounded  by  elastic  cement  (Figure  3d). 
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Stress  distributions,  such  as  shown  in  Figures  3b  and  3c,  are  obtained  under  the 
assumption  that  the  normal  and  shear  stresses  in  the  elastic  cement  can  be  estimated  using 
the  elastic  foundation  approximation  (Dvorkin  et  ah,  1991).  This  approximation  neglects 
lateral  strain  within  the  cement  layer  and  thus  is  not  adequate  to  describe  stress  distribution 
at  the  periphery  of  the  cement  layer.  As  a  result,  there  may  be  a  possibility  that  plastic  yield 
initiates  from  the  periphery  inward,  instead  of  from  the  center  outwards,  as  assumed  here. 
It  is  clear  that  a  rigorous  analysis  of  contact  stresses  within  cement  is  needed  in  order  to 
refine  the  first-order  approximation  model  presented  here. 

When  deriving  the  plasticity  condition  for  a  thin  cement  layer  we  again  assume  that  the 
normal  and  shear  stresses  in  the  elastic  part  can  be  estimated  from  the  elastic  foundation 
approximation.  In  the  plastic  part  of  the  cement  layer  we  use  the  von  Mises  yield  criterion. 
A  word  of  caution  is  that  this  isotropic  criterion  may  be  not  accurate  for  geomaterials  with 
layered  cement,  such  as  laminated  clay. 

Once  the  plasticity  condition  is  established,  we  find  elastic  and  plastic  contact  stresses 
between  two  grains.  This  problem  requires  solving  an  ordinary  integral  equation.  The 
solution  is  numerical  and  straightforward.  These  contact  stress  distributions  allow  us  to 
find  normal  and  shear  stiffnesses  of  a  two-grain  combination. 

Intrinsic  in  this  analysis  is  the  assumption  of  local  axisymmetry.  Dvorkin  et  al.  (1994) 
show  that  this  assumption  is  aecurate  for  an  elastic  cement  layer.  Additional  future  effort  is 
needed  to  estimate  the  accuracy  of  the  local  axisymmetry  assumption  in  elastic-plastic 
cement. 

Next  we  use  statistical  modeling  to  calculate  the  effective  macroscopic  stresses  in  a 
random  pack  of  identical  cemented  spheres  subject  to  high  strains.  Theoretical  stress-strain 
curve  shapes  are  very  close  to  those  observed  on  real  geomaterials. 

One  important  effect  is  that  an  aggregate  with  elastic-plastic  cement  becomes  anisotropic 
if  macroscopic  stresses  are  non-hydrostatic.  This  effect  follows  from  the  fact  that  the 
contact  stiffnesses  are  strain-dependent  since  larger  local  strains  result  in  larger  plastic 
zones.  This  stress-induced  anisotropy  may  produce  an  up  to  10-percent  error  in 
estimating,  for  example,  an  apparent  isotropic  bulk  modulus  from  a  uniaxial-strain 
experiment. 


Figure  2.  Dynamic-to-static  bulk  modulus  ratios  for  sandstones.  The  labels  on  the  curves  indicate 
volumetric  clay  content. 
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3.0.  ELASTIC-PLASTIC  CEMENT  LAYER 


3.1.  Yield  Criterion 

As  a  yield  criterion  we  use  the  von  Mises  criterion  of  constant  intensity  of  tangential 
stresses  (Kachanov,  1971): 

r=<7,/V3,  (3.1) 

where  T  is  the  intensity  of  tangential  stresses: 

T  =  +  aj  +  6(<  +  +  (3.2) 

are  the  components  of  the  stress  tensor,  and  rr,  is  the  yield  limit  of  the  cement  under 

tension.  .  . 

Equations  (3.1)  and  (3.2)  give  the  following  expression  for  the  von  Mises  critenon: 

aJ  +  aJ  +  +3(<  +  oj  +  £7,/)-  (o-„(T„  +  <7„cr„  +  =  <t/.  (3.3) 

This  last  equation  can  be  also  expressed  in  terms  of  principal  stresses  <7^,  o^,  and  : 


Figure  3.  a.  Contact  stresses  between  two  cemented  grains,  b.  Normal  stresses  between  two 
cemented  grains  with  soft  cement,  c.  Shear  stresses  between  two  cemented  grains  with  soft  cement, 
d.  Cement  layer  with  a  plastic  central  kernel  surrounded  by  elastic  cement.  Stresses  are  normalized  by 
the  average  stress.  Coordinate  is  normalized  by  the  radius  of  the  cement  layer. 
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3.2.  Thin  Cement  Layer  in  Extension  and  Shear 


Consider  a  thin  cement  layer  (Figure  4a)  whose  width  2h  may  gently  change  along  its 
length:  h"(x)«l.  The  z  =  0  plane  of  the  rectangular  (x,y,z)  coordinate  system  is 
parallel  to  the  layer.  The  layer’s  walls  have  displacements  W(x)  that  are  parallel  to  the  z 
axis.  We  assume  that  as  long  as  the  layer  deforms  elastically,  stresses  can  be  estimated 
from  the  elastic  foundation  approximation  (Dvorkin  et  al,  1991  and  1994): 


2Gv  Wix) 


^ XX  ^yy  1  n..  ’  ^zz 


l-2v,  h(x) 


2GSl-vJW(x) 
l-2v^  h(x)' 


(3.4) 


where  G^  and  are  the  shear  modulus  and  the  Poisson's  ratio  of  the  cement,  respectively. 

By  substituting  equation  (3.4)  into  the  criterion  (3.3)  we  arrive  at  the  following 
plasticity  condition  for  a  thin  cement  layer  in  extension  (or  compression): 


W(x,)  _  C7,  ^ 
Kxq)  2G^ 


(3.5) 


where  Xq  is  the  coordinate  at  the  elastic-plastic  boundary.  In  equation  (3.5)  we  neglect 
hardening  and  assume  that  the  yield  limit  in  compression  is  the  same  as  in  tension. 

Consider  now  a  plastic  zone  surrounded  by  two  elastic  zones  (Figure  4c).  Again  we 
assume  that  the  elastic  foundation  formulas  (3.4)  are  valid  in  the  elastic  zone.  By  doing  so 
we  neglect  the  condition  of  incompressibility  that  will  influence  stresses  at  the  elastic-plastic 
boundary.  The  equations  of  balance 


dx  By  Bz 


Bo^ 

=  0,-^ 

ox 


<9(7., Be. 


zy 


By  Bz 


=  0, 


Ba^,  .  Ba^^  Bo., 


Bx 


By  Bz 


=  0 


will  be  satisfied  if  we  assume  that  all  stresses  are  constant  within  the  plastic  zone.  If  the 
plastic  zone  is  symmetrical  with  respect  to  its  center,  then  on  its  boundaries  AA  and  BB 
(Figure  4c)  where  condition  (3.5)  is  met,  stresses  will  be  the  same: 


2G,v,  Wix)  (T,v, 


i_2v,  Kx)  l-2v/ 

_  2G,(1-vJW(x)_ct,,(1-vJ_ 

<7„  -  ^  I  rT~  -  ,  ~ 


l-2v^  h(x)  l-2v^ 


Figure  4.  a.  Thin  cement  layer  in  extension,  b.  Thin  cement  layer  in  shear,  c.  Thin  cement  layer  - 
plastic  zone  surrounded  by  elastic  cement. 
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Therefore,  within  the  plastic  zone  stresses  are  constant: 


^xx  -  ^yy  - 


Q  —  _ ^•<7  =0 


(3.6) 


Consider  next  a  thin  cement  layer  deforming  in  shear  with  tangential  displacement  IJ (x) 
(Figure  4b).  From  the  elastic  foundation  approximation  we  have  the  following  stresses  in 
the  elastic  part  of  the  layer: 


_  GTiJ{x) 
h{x) 


(3.7) 


with  all  other  stresses  being  zero.  Now  the  von  Mises  criterion  (3.3)  yields  the  following 
plasticity  condition: 


<7(Xo)  ^  0-, 

/z(Xo)  V3G^ 


(3.8) 


Again,  putting  stresses  in  the  plastic  part  of  the  layer  constant,  we  find  that  within  the 
plastic  zone 


«^x^  =  o-./V3,  (3.9) 

with  all  other  stresses  being  zero. 

3.3.  Thin  Cement  Laver:  Combined  Extension  and  Shear 


If  the  layer  deforms  simultaneously  in  extension  and  shear,  the  elastic  foundation 
approximation  gives  the  following  stresses  within  the  elastic  part: 


^  W(x)  2G,(1-v,)W(x)  G,U{x) 

^  \-2v,h{x)'  "  l-2v,  Mx)’  h{x)  ’ 


(3.10) 


with  all  other  stresses  being  zero.  From  equations  (3.10)  and  (3.3)  the  plasticity  condition 
follows  as 


4 


Wjx,) 

KXq) 


G(Xo) 

/j(Xo) 


(3.11) 


and  the  stresses  within  the  plastic  zone  are 


— 

«  »  i_2v 


2G.V.  W(Xo)_ 
/l(Xo)  ’ 


2G.(l-v,)W(Xo)  ^ 

l-2v,  Mxo)’  ^ 


G,G(Xo) 

KXq) 


(3.12) 


It  is  clear  that  if  stresses  are  axisymmetrical  then  a  plasticity  kernel  surrounded  by  the 
elastic  material  will  be  axisymmetrical  as  well.  In  this  case  we  can  use  the  radial  coordinate 
r  instead  of  x  in  equations  (3.5),  (3.6),  (3.8),  (3.9),  (3.11),  and  (3.12). 
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4.0.  TWO  CEMENTED  SPHERES 


4.1.  Assumptions 

Consider  two  identical  elastic  spheres  (grains)  cemented  at  their  contact.  The  cement 
layer  is  much  thinner  than  the  sphere  radius.  The  configuration  is  axisymmetrical  with 
respect  to  the  line  connecting  the  centers  of  the  spheres.  The  spheres  are  subjected  to 
normal  and  tangential  displacements  with  respect  to  the  contact  point  (Figure  5a).  The 
cement  may  deform  in  both  elastic  and  plastic  modes  (Figure  5b). 

We  assume  that  shear  stresses  that  are  imposed  on  the  grain  surfaces  due  to  the 
tangential  displacements  of  the  spheres,  and  also  produced  by  the  normal  displacements  do 
not  affect  normal  stresses  on  the  grain  surfaces.  Also,  normal  stresses  that  are  imposed  on 
the  grain  surfaces  due  to  the  normal  displacements  of  the  spheres,  and  produced  by  the 
shear  displacements  do  not  affect  shear  stresses  on  the  grain  surfaces.  A  similar 
assumption  is  inherently  contained  in  the  elastic  foundation  model  that  we  use  for  the 
cement  layer.  As  far  as  the  grains  are  concerned,  the  assuniption  can  be  supported  by  some 
contact  mechanics  solutions,  e.g.,  a  rigid  punch  on  an  elastic  half-space  (Johnson,  1985). 

In  addition  we  assume  that  all  stresses  are  axisymmetrical  with  respect  to  the  line 
connecting  the  centers  of  the  spheres.  Dvorkin  et  al.  (1994)  show  that  this  assumption  is 
accurate. 


Figure  5.  a.  Two  cemented  elastic  spheres  subject  to  normal  and  tangential  deformation,  b.  Thin 
cement  layer  between  the  spheres. 


4.2.  Stresses  in  Cement 

The  cement  layer  may  have  a  plastic  kernel  and  an  elastic  ring  surrounding  this  kernel 
(Figure  5b).  The  normal  W(r)  and  the  tangential  U{r)  displacements  of  the  surface  of  the 
cement  layer  are  axisymmetrical  (with  respect  to  the  line  connecting  the  centers  of  the 
spheres,  and  depend  on  the  radial  coordinate  r . 

To  relate  W{r)  and  U{r)\.o  the  normal  p{r)  and  tangential  q{r)  stresses  acting  on  the 
grain  surface  in  the  cemented  area,  we  use  equations  (3.10)  and  (3.12).  In  the  elastic  part 
of  the  cement  we  have 


p{r)  =  - 


2G,(1- V,)W(r) 
l-2v,  h{r) 


q(r)  =  - 


GMr). 

h(r) 


hir)  =  R[£  +  ^(-^fle  = 
L  K 


h(0) 

R  ’ 


(4.1) 


where  e  is  the  minimum  half-thickness  of  the  cement  layer.  In  the  plastic  region  these 
stresses  do  not  depend  on  the  radial  coordinate: 
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(4.2) 


The  normal  displacement  of  the  grain  surface  w{r)  and  the  tangential  displacement  of 
the  grain  surface  u{r)  can  be  approximately  related  to  the  normal  p{r)  and  tangential  q{r) 
stresses  as  (Dvorkin  et  al.,  1994) 


where  G  and  v  are  the  grain  shear  modulus  and  Poisson's  ratio,  respectively.  The 
integration  is  conducted  inside  the  circle  r<b  (Figure  6),  where  b  is  the  radius  of  the 
cement  layer.  Note  that  w(r)  and  u{r)  are  displacements  relative  to  the  sphere's  center. 


Figure  6.  The  region  of  cemented  contact  on  the  grain  surface. 

These  displacements  can  be  related  to  the  displacements  W(r)  and  U{r)  of  the  surface 
of  the  cement  layer  as 

w{r)  =  W(r)  +  5,  u{r)  =  U{r)  +  T,  (4.5) 

where  5  and  T  are  the  normal  displacement  and  the  tangential  displacement,  respectively, 
of  the  sphere's  center  relative  to  the  contact  point  (Figure  5a). 

4.4.  Contact  Stresses  and  Stiffnesses 

By  combining  equations  (4.1)  -  (4.5)  we  arrive  at  the  following  system  of  equations  for 
the  contact  stresses  p{r)  and  q{r) : 
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.  -  TT  rcos(p+^jh^~r'^s^n^^ _ 

— ~  - h(r)p(r)+ \d(p  f  p(-\jr^+s^-2rscos(p)ds  - 

2G,(l-v,)(l-v)  ^  i  ^  i  ^ 

JtG  ^  , 

- 0  at  c<r<b, 

1-  V 


ji  rcos(p+^jb^ -r^ sin^  (p _ 

—h(r)q{r)  +  jd(p  jq{^jr'^+s^-2rscos(p)(l-  vsin^ (p)ds  = 

0  0 

kGt  at  c<r<b; 


(4.6) 


p(r)  =  p(c),  q{r)  =  q{c)  at  0  <  r  <  c; 


Pi.c) 


1-2V, 

1-V, 


+  3^^(c) 


2 


System  (4.6)  can  be  easily  solved  numerically  using,  for  example,  the  quadrature  method 
(Delves  and  Mohamed,  1985). 

We  define  the  normal  (S„)  and  the  tangential  (5^)  contact  stiffnesses  of  a  two-grain 
combination  as 


S,=FJS,S,  =  F,/r,  (4.7) 

where  and  are  the  resulting  normal  and  tangential  forces,  respectively: 

h  b 

F„  =  27rj p(r)rdr,  F^  =  2;r| q{r)rdr. 

0  0 


4.5.  Examples 

We  explore  contact  stress  distribution  between  two  quartz  grains  cemented  by  clay. 
The  ratio  of  the  cement  radius  to  the  grain  radius,  &  /  /?,  is  0.7.  The  separation  distance 
between  the  grains  at  their  contact,  e,  is  zero.  The  elastic  constants  are  chosen  as  follows: 
G  =  45  GPa,  V  =  0.064,  G^  =  1.98  GPa,  =  0.29  (Yin,  1993).  The  yield  limit  of  the 
cement,  (T^,was  chosen  at  450  MPa,  following  the  failure  experiments  by  Jizba  (1991)  on 
high-clay-content  sandstones. 

In  Figure  7  we  present  normal  stresses  along  the  radius  of  the  cement  layer  at  varying 
normal  force  As  this  force  increases,  the  plastic  kernel  starts  developing  at  the  center 
of  the  contact.  In  this  kernel  normal  stresses  are  constant  as  given  by  formula  (3.6). 
Outside  the  plastic  kernel  normal  stresses  gradually  decrease  towards  the  periphery  of  the 
cement  layer. 

Similar  shapes  are  obtained  for  tangential  contact  stresses  at  varying  tangential  force  F^ 
(Figure  8).  In  this  case  tangential  stresses  in  the  plastic  kernel  are  constant  as  given  by 
formula  (3.9).  The  decrease  of  the  tangential  stresses  towards  the  periphery  of  the  cement 
layer  is  steeper  than  that  of  the  normal  stresses. 

Next  we  explore  the  effect  of  grain  separation,  e  =  /j(0)  !  R,  on  the  contact  stresses. 
We  compare  the  case  with  no  separation  (g  =  0)  to  the  case  where  £  =  0.1  (Figure  9). 
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Normal  stress  distribution  appears  to  be  more  uniform  in  the  case  where  the  grains  are 
separated  as  compared  to  the  no-separation  case. 

By  changing  the  b!  R  ratio,  and  the  yield  stress  value,  one  can  arrive  at  the  results 
qualitatively  identical  to  those  presented  in  Figures  7  -  9. 


0  Center -I - ^ - 1 - Periphery 

0.2  0.6  1 
Coordinate  along  cement  layer 


0  Center _j - 1 - ^ - Periphery 

0.2  0.6  1 
Coordinate  along  cement  layer 


Coordinate  along  cement  layer 


Coordinate  along  cement  layer 


Figure  7.  Normal  stress  distribution  along  a  cement  layer.  The  coordinate  is  normalized  by  the 
layer’s  radius.  The  stresses  are  normalized  by  the  grain’s  shear  modulus.  The  plastic  kernel 
progresses  with  increasing  load  (left  to  right,  top  to  bottom). 
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Normal  Stress  Shear  Stress  Shear  Stress 


0.2  0.6  1 


Coordinate  along  cement  layer 


0.2  0.6  1 


Coordinate  along  cement  layer 


0,2  0.6  1 


Coordinate  along  cement  layer 


0.2  0.6  1 


Coordinate  along  cement  layer 


Figure  8.  Tangential  stress  distribution  along  a  cement  layer.  The  coordinate  is  normalized  by  the 
layer’s  radius.  The  stresses  are  normalized  by  the  grain’s  shear  modulus.  The  plastic  kernel 
progresses  with  increasing  load  (left  to  right,  top  to  bottom). 


0.2  0.6  1  0.2  0.6  1 


Coordinate  along  cement  layer  Coordinate  along  cement  layer 

Figure  9.  Normal  stress  distribution  along  a  cement  layer.  Grains  are  separated  ( £  =  1)  or  have  a 
point  contact  (£  =  0).  The  coordinate  is  normalized  by  the  layer's  radius.  The  stresses  are 
normalized  by  the  grain's  shear  modulus.  The  plastic  kernel  progresses  with  increasing  load  (left 
to  right). 

5,0.  EFFECTIVE  MODULI  OF  AN  AGGREGATE 


5.1.  Assumptions  and  Definitions 


Consider  a  randona  pack  of  identical  elastic  spheres  of  radius  R  cemented  with  elastic- 
plastic  cement.  The  pack  is  statistically  homogeneous  and  isotropic.  The  average  number 
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n  of  contacts  per  sphere  (coordination  number)  is  about  9.  The  spheres  occupy 
approximately  64  percent  of  the  pack's  volume  which  means  that  porosity  of  the  pack  (not 
counting  cement)  is  about  0.36  (Bourbie  et  al.,  1987).  In  our  further  derivations  we  follow 
the  statistical  scheme  used  by  Walton  (1987). 

The  i-th  sphere's  center  has  position  vector  in  the  rectangular  (1,2,3)  coordinate 

system  (Figure  10a).  The  displacement  of  the  i-th  center  is  u^''\  The  unit  vector  7^"'"’  that 
connects  the  center  of  the  m-th  sphere  and  the  center  of  the  n-th  sphere  (Figure  10b)  is 

Y(«)  _  y(w) 

/(»»■)  _ d —  (5.1) 

2R 

We  also  assume  that  the  displacements  of  the  sphere  centers  are  consistent  with  the 
applied  uniform  strain  field  e-^ : 

=  (5.2) 


Figure  10.  Displacements  and  forces  in  a  sphere  pack. 

5.2.  Relative  Deformation  of  Grains 

The  displacement  vector  of  the  m  -th  sphere's  center  relative  to  the  contact  point  is 

2 

By  decomposing  this  vector  into  the  normal  part  Wq,  along  and  into  the  tangential 
part  Uq,  normal  to  we  find: 


Wo-/ 


Ur) 


-I 


(nm) 


(5.3) 


Formulas  (5.1)  -  (5.3)  yield  the  following  expressions  for  the  absolute  values  and  for 
the  components  of  vectors  Wq  and  : 


w,=  5  =  Uo  =  t^R 


■'pq  p 


^  T  («^)  _  ^  T  («^)r  (nm)j{nm)\2  . 
Z^^^ip^p  ^pq^p  ^q  i  ^ 


/=l 


(5.4) 


w, 


0/ 
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5.3.  Contact  Forces 


We  define  the  normal  contact  stiffness  between  two  spheres  as  the  ratio  of  the 
normal  component  of  the  contact  force  to  the  normal  displacement  Wq  (Figure  10c). 
The  tangential  contact  stiffness  is  defined  similarly: 

S.  =  FJw„,S,  =  F,lu„.  (5.5) 

Generally,  each  of  these  stiffnesses  depends  on  both  normal  and  tangential  displacement. 

It  follows  from  formulas  (5.4)  and  (5.5)  that  the  components  of  the  contact  force  vector 

in  the  (1,2,3)  coordinate  system  are 


’^r)^pq^p 


(nm)  T  {nm)j(nm) 


(nm) 


]. 


(5.6) 


5.4.  Macro.scopic  Stresses 

We  define  average  macroscopic  stresses  in  volume  V  that  is  occupied  by  a  sphere  pack 
as 


Any  average  value  is  meant  to  be  a  volumetric  average  as  in  the  latter  equation.  By  means 
of  the  Gauss  theorem  this  volume  integral  can  be  transformed  into  a  surface  integral  as: 


V  s 

where  S  is  the  surface  of  volume  V,  x,.  are  coordinates  on  the  surface,  and  t.  are  surface 
traction  components. 

Consider  now  volume  that  contains  the  «-th  sphere.  By  integrating  the  stress- 
tensor  components  inside  this  volume  we  find: 


=i  {(a:,!,.  +x,t,)dS,=^j(.x:ti  +x’/,)dS., 


(5.7) 


where  x'  =  3c-X'"\  The  last  equity  in  equation  (5.7)  is  due  to  the  fact  that  surface 
tractions  are  self-equilibrating  and  thus  their  surface  integral  is  zero.  Tractions  on  the 
sphere's  surface  are  contact  forces  that  act  at  the  contacts  only.  Therefore, 

x'  =  (X/"'^  -  X/"^)  /  2  =  (5-8) 

as  follows  from  equation  (5.1).  By  substituting  equation  (5.8)  into  equation  (5.7)  we  find 
that  for  one  sphere 


J  =  -R  /  2. 

y  Contacts 


Therefore,  average  stress  in  the  aggregate  is 


(5.9) 
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(5.10) 


(nm)  p  (.nm)  j  (.nm)  p  (nm) 

^  All  Contacts 

where  the  factor  1/2  from  equation  (5.9)  disappears  since  the  summation  in  equation  (5.10) 
is  through  all  spheres  and  all  contacts  on  each  sphere,  and  thus  each  contact  is  counted 
twice.  Next  we  substimte  equation  (5.6)  into  equation  (5.10)  and  find  that  average  stress 
in  the  aggregate  is 


K)  =  -V  2  + 

^  All  Contacts 

If  the  average  number  of  contacts  per  sphere  is  n  and  the  number  of  spheres  is  N,  then 
the  total  number  of  contacts  \s  nN  12.  The  volume  occupied  by  N  spheres  is 


4  NkR^ 

3  !-(/> 

where  (p  is  the  porosity  of  the  aggregate.  Now,  by  introducing  statistically  average 
quantities,  we  can  replace  equation  (5.11)  with  the  following  equation  for  the  macroscopic 
stress  tensor  components: 

oKK 

(5,(Wo,Mo)(ey//  +  (5-12) 

_ 


5.5.  Averaging  Technique 

The  expression  that  has  to  be  averaged  in  equation  (5.12)  is  a  scalar  function  of  a  unit 
vector  I  Consider  such  a  function  0(7)  =  <&(/i,/2,/3)  in  a  spherical 

coordinate  system  (r,  6,  (p)  (Figure  1  la).  In  this  coordinate  system 

/j  =  sin  6  cos  (p,  I2  =  sin  6  sin  (p,  /j  =  cos  0. 


We  assume  that  the  position  of  the  head  of  the  vector  I  is  uniformly  distributed  on  the 
surface  of  a  unit-radius  sphere  (Figure  1  lb).  Then  the  average  of  function  O  is 


In  n 


(o)  =  1  f  <D Ju  =  —  j  J  0(/„  /2 ,  /j )  sin  Odcpde  = 


In  n 

cos  <p,  sin  0  sin  cos  6)  sin  Bd<pdd, 

0  0 


(5.13) 
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where  Z  is  the  surface  of  the  unit  sphere,  and  t?  is  an  element  of  the  surface  (Figure  1  lb). 
By  using  formula  (5.13)  we  can  find,  for  example,  that 


(5.14) 


Figure  1 1 .  Spherical  coordinate  system  for  averaging  a  function  of  a  unit  vector. 

5.6.  Constant  Contact  Stiffnesses 

In  this  section  we  derive  the  effective  elastic  constants  of  an  aggregate  where  eement  is 
purely  elastic  and  thus  contact  stiffnesses  S„  and  are  stress-independent  (Dvorkin  et  al., 
1994).  In  this  case  formula  (5.12)  can  be  written  as 


(5.15) 


Consider  a  hydrostatic  stress  field  where  strain-tensor  components  are  given  by 
e  =  e5  .  Then  formulas  (5.14)  and  (5.15)  yield 


3n(l-0) 


sA>")= 


n(l  —  0)c 


4;ri?  "  '  ■  '  47tR 

Thus  the  effective  bulk  modulus  is 


.  «(l-«^)5„ 

^  \1kR 


(5.16) 


Consider  now  a  pure  shear  field  where  e,2  =^21“^’  other  strain-tensor 

components  are  zero.  Formula  (5.15)  now  yields 


'12 


=  _  s,)t{iK)  +  S,r{/,=  +  /,“>]  =  ^^(.5.  + 1^.)- 


?>7tR 


Thus  the  effective  shear  modulus  is 


^(1-^) 

20kR 


(5.17) 


Formulas  (5.16)  and  (5.17)  are  consistent  with  those  given  by  Winkler  (1983). 
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6.0.  EXAMPLES  OF  AGGREGATE  DEFORMATION 

6.1.  Hydrostatic  Compression 

Consider  a  random  pack  of  identical  elastic  spheres  cemented  at  their  contacts.  The 
cement  is  initially  elastic  but  becomes  plastic  as  contact  loads  increase.  Consider 
hydrostatic  compression  where  strain-tensor  components  are  given  as  =  -eS-j.  Equation 
(5.12)  gives  in  this  case: 

WQ  =  6  =  Re,Uo  =  t  =  0,(Jij  =  S„iRe,0)5y. 

Therefore,  it  follows  from  equations  (5.12)  and  (4.6)  that 


(7,1  _  n(l-^)g^ 

G  2  " 


(6.1) 


and  function  f(t)  can  be  found  from  equation 


2  jr  fcos(p+vp  -rsm  (p _ 

+  j d(p  J f (y/^  cos (p )d^  —  (52) 

2  0  0 

m  =  f{x\  0<t<;t. 

where  x  is  found  from  equation 


fix) 


l-2v, 

1-v, 


G  ■ 


The  notations  in  equations  (6.1)  -  (6.3)  are: 


(6.3) 


;r(l-2v,)G 

2(l-v,)(l-v)G, 


1  ^  o  ^ 


f(t\=PiM  t  =  -  t=-  S  = 

G  ’  R’^  R’^"  2nRG' 


I 


We  use  equation  (6.1)  to  calculate  the  hydrostatic  deformation  of  a  random  pack  of 
identical  cemented  spheres  with  elastic  and  plastic  parameters  as  in  the  example  in  Section 
4.5.  Again,  we  explore  the  effect  of  grain  separation,  e  =  h(0)  /  i?,  on  the  stress-strain 
laws  (Figure  12a). 

Separation  of  grains  with  e  =  0.1  results,  as  could  be  expected,  in  smaller  macroscopic 
stresses  at  given  strain,  as  compared  to  the  no-separation  case  e  =  0 .  The  elasticity  region 
in  the  stress-strain  plane  is  larger  for  £  =  0.1  than  for  e  =  0.  This  is  a  result  of  reduced 
contact  stress  concentration  at  the  center  of  contact  in  the  former  case  (Figure  9). 

The  effective  bulk  modulus  is  constant  in  the  elastic  region  and  gradually  decreases  as 
the  plastic  kernel  progresses  from  the  center  of  the  cement  layer  (Figure  12b).  This 
modulus  is  small  where  the  grains  are  separated  as  compared  to  the  case  of  contacting 
grains. 
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6.2.  Uniaxial-Strain  Compression 


Consider  a  uniaxial-strain  compression  of  a  random  sphere  pack  with  gjj  =  -e  and  all 
other  strain-tensor  components  vanishing.  Equations  (5.12)  and  (4.6)  give  in  this  case 
(using  the  averaging  technique  described  in  Section  5.6): 

CTi,  _(722  _  2>n{\-^)e 

G  G  A 

nil  _ 

j [S„(ecos^  0,ecos0sin  6)  -  S^(ecos^  d,ecos6sin0)]cos^  Osin^  Odd, 

0 


g-33  ^  3«(l-<^)g^^ 

G  2 

nil  __ 

J[5„(ecos^  0, c cos 0 sin 0)cos^  0-1-5^ (ccos^  0,ccos0sin0)sin^  0]cos^  6s,m0dd. 
0 

where 


5„  = 


1 

ccos^  6 


jf{t)tdt,S^  = 
0 


1 

ccos0sin0 


Functions  f(t)  and  g(t)  are  found  from  equations 

2  )r  <  cos  <p+^lpi^-t^  sin^  <p _ 

A^(e  +  —)f{t)  +  jd(p  j  -2t^ cos (p)d^  =  X„e cos^  0, 

^0  0 

2  n  t  cos  <P+-\]p^ sin^  (p _ 

A,(e  +  +  jd(p  j g(^|t^  +  -2t^  cos  (p)(\  -  vsin^  (p)d^  = 

^0  0 

A^ecos0sin0, 

m=fixi  g(t)=g{xx  o<t<x. 

The  value  of  x  is  found  from  equation 

f(x)(^^^y + ^8"(x) = (%'• 

We  use  the  above  equations  to  calculate  macroscopic  stresses  in  a  random  pack  of 
identical  cemented  spheres  from  the  previous  example.  The  grains  are  in  direct  contact. 
The  axial  and  radial  stresses  are  plotted  versus  the  axial  strain  in  Figure  13.  The  plastic 
deformation  region  is  clearly  visible  in  the  <733  versus  strain  curve.  However  the  cr,j 
versus  strain  curve  is  practically  straight,  apparently  due  to  smaller  stresses  acting  in  the 
radial  direction. 
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Figure  12.  Hydrostatic  compression  of  a  random  pack  of  cemented  elastic  grains.  Grains  are  in 
directcontact  (e  =  0)  or  are  separated  ( e  =  0.1).  a.  Stress-strain  curves,  b.  Effective  bulk 
modulus.  Dots  indicate  the  boundary  of  the  elastic  deformation  region. 


Figure  13.  Uniaxial-strain  compression  of  a  random  pack  of  cemented  elastic  grains.  Grains  are  in 
direct  contact.  The  axial  and  radial  stresses  are  plotted  versus  the  axial  strain.  The  dot  indicates  the 
boundary  of  the  elastic  deformation  region. 

6.3.  Complex  State  of  Stress.  Real  Geomaterials 

Formula  (5.12)  gives  a  straightforward  expression  for  macroscopic  stresses  if  the 
boundary  conditions  are  formulated  in  strains,  such  as  in  the  above  examples  of  hydrostatic 
compression  and  uniaxial-strain  compression.  The  situation  becomes  more  complicated  if 
the  boundary  condition  is  formulated  partly  in  stresses  as,  for  example,  in  a  uniaxial-stress 
compression  experiment.  In  this  case  the  lateral  stress,  tr,  ,  is  constant.  A  scheme  of 
solving  such  a  problem  is  as  follows.  Use  equation  (5. 12)  with  =  -^3 ,  =  ^22  =  ^1 » 

and  all  other  strain-tensor  components  vanishing.  By  computing  tTi,  select  Cj  (at  fixed  e^) 
such  that  CT„  =  The  axial  stress  computed  at  the  determined  e,,  and  ^33  is  the 

desired  (733(^33). 
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Expressions  (4.7)  for  the  normal  and  tangential  eontact  stiffnesses  of  a  two-grain 
combination  can  be  directly  plugged  into  a  numerical  discrete  element  scheme  to  calculate 
the  effective  macroscopic  laws  in  realistic  geomaterials  under  a  complex  state  of  stress 
(e.g.,  Trent  and  Margolin,  1992). 

Certain  confidence  in  the  appropriateness  of  our  micromechanical  model  for  an  elastic- 
plastic  deformation  of  cemented  grains  follows  from  the  fact  that  the  solution  produces 
stress-strain  laws  similar  to  those  experimentally  obtained  in  ductile  rocks  (Figure  1). 

6.4.  Stress-Induced  Anisotropy 

Due  to  cement’s  plastic  deformation,  normal  and  tangential  contact  stiffnesses  are 
stress-dependent.  This  effect  may  lead  to  stress-induced  anisotropy  in  an  otherwise 
isotropic  aggregate.  Consider  two  experiments  performed  on  the  same  sample:  hydrostatic 
compression,  and  uniaxial-strain  compression.  In  the  first  experiment  the  effective  bulk 
modulus  can  be  calculated  from  the  stress-strain  curve  as 

_  ^^11  _  ^^22  _  ^^33  (5  4) 

Hydrostadc  -  3^^^^  3^^^^ ' 

If  we  assume  that  the  aggregate  is  isotropic  in  the  second  experiment,  then  the  effective 
bulk  modulus  is: 

d(<7^^  -f  2(T)])  (6.5) 

^Uniaxial  o 

By  comparing  the  "hydrostatic"  stress  from  the  uniaxial-compression  experiment, 

cr33-i-2cr„^  with  the  true  hydrostatic  stress,  (7„,  from  the  hydrostatic-compression 
experiment  (Figure  14a)  we  see  that  the  uniaxial-compression  "hydrostatic"  stress  is  larger. 
The  results  of  calculating  the  effective  moduli  as  given  by  formulas  (6.4)  and  (6.5)  show 
that  the  uniaxial  bulk  modulus  is  larger  than  the  "true"  hydrostatic-compression  bulk 
modulus  (Figure  14b).  At  high  strain  this  difference  may  become  as  large  as  10  percent. 
This  observed  disparity  is  an  indication  of  stress-induced  anisotropy  in  the  uniaxial- 
compression  experiment. 

7.0.  CONCLUSION 

A  theoretical  solution  is  given  to  the  problems  of  normal  and  tangential  elastic-plastic 
deformation  of  two  cemented  elastic  grains.  Cement  serves  as  an  elastic-plastic  element. 
As  resulting  forces  between  two  grains  increase,  stress  concentration  develops  in  the  center 
of  the  cement  layer.  Then  a  plastic  kernel  develops  as  local  stresses  meet  a  plasticity 
condition.  This  plasticity  condition  for  a  thin  elastic-plastic  cement  layer  is  derived  from 
the  von  Mises  criterion.  The  plastic  kernel  in  the  cement  layer  is  surrounded  by  an  elastic 
ring  where  stresses  and  strains  are  linked  through  the  elastic  foundation  model.  The 
solution  is  reduced  to  ordinary  integral  equations.  It  gives  normal  and  tangential  stiffnesses 
of  a  two-grain  combination  for  large  microscopic  strains. 

This  microscopic  model  is  used  to  obtain  macroscopic  stress-strain  laws  in  elastic- 
plastic  packs  of  identical  cemented  spheres.  To  do  so  we  statistieally  derive  effective  stress- 
strain  laws  for  an  aggregate  of  cemented  spheres  where  grain-to-grain  contact  stiffnesses 
are  strain-dependent. 

Theoretical  stress-strain  laws  thus  obtained  are  close  to  the  relations  typically  observed 
in  clay-eemented  sands.  Therefore,  our  microscopic  model  ean  be  used  to  realistically 
model  large-strain  deformation  of  cemented  geomaterials. 
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An  important  result  is  stress-induced  anisotropy  in  initially  isotropic  aggregates.  Such 
an  aggregate  may  become  anisotropic  if  the  stress  field  is  non-hydrostatic.  This  stress- 
induced  anisotropy  may  lead  to  an  up  to  10-percent  error  in  estimating,  for  example,  an 
apparent  isotropic  bulk  modulus  from  a  uniaxial-strain  experiment. 


Figure  14.  Stress-induced  anisotropy  in  a  random  pack  of  identical  cemented  spheres,  a.  "Hydrostatic" 
stress  versus  strain  from  the  uniaxial  compaction  and  the  hydrostatic  compaction  experiments,  b. 
Bulk  moduli  from  the  uniaxial  compaction  and  the  hydrostatic  compaction  experiments.  The  dots 
indicate  the  boundary  of  the  elastic  deformation  region. 
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2.8.  PRESSURE  SENSITIVITY  OF  CEMENTED  GRANULAR 

MATERIALS 


ABSTRACT 

In  this  paper,  we  explore  the  mechanical  contact  interaction  of  two  identical  elastic 
spheres  uniformly  coated  with  thin  layers  of  a  different  elastic  material.  These  two  coating 
layers  intersect  over  a  finite  contact  area  thus  bonding  the  spheres.  The  normal  contact 
stiffness  and  the  shear  contact  stiffness  increase  when  the  spheres  are  axially  pressed 
together,  due  to  the  increasing  contact  area.  We  calculate  the  dependence  of  these 
stiffnesses  on  the  axial  load  by  using  a  new  approximate  analytical  solution.  The  solution 
also  gives  the  distributions  of  the  normal  and  shear  stress  components  on  the  cemented 
contact.  We  use  this  solution  to  calculate  the  pressure  dependence  of  the  effective  elastic 
moduli  of  a  random  pack  of  identical  cemented  spheres.  This  pressure  dependence  may  be 
large  if  the  initial  contact  radius  is  small.  It  is  insignificant  for  large  contact  radii.  If  the 
spheres  are  in  direct  contact  and  the  initial  contact  radius  is  small,  the  elastic  properties  of 
the  cement  have  little  effect  on  the  pack's  elastic  moduli.  However,  if  the  spheres  are 
separated  by  even  a  small  cemented  gap,  the  elastic  properties  of  the  cement  may  have  a 
considerable  effect  on  the  pack's  elastic  moduli. 

1.  INTRODUCTION 

The  mechanical  properties  of  brittle  granular  materials  are  strongly  affected  by  local 
intergranular  effects,  such  as  the  geometrical  nonlinearity  due  to  pressure  dependence  of  the 
contact  area.  These  nonlinear  effects  prevail  even  if  the  material  consists  of  linearly  elastic 
grains.  Granular  materials  are  often  modeled  as  a  pack  of  linearly  elastic  spheres,  and  the 
fundamental  problem  of  the  contact  between  two  elastic  spheres  has  been  the  focus  of  many 
studies  for  over  a  century.  Notable  are  the  solutions  of  stress  distribution  and  contact 
force-displacement  relations  for  the  problems  of  a  relative  axial  displacement  (Hertz,  1882); 
a  relative  tangential  displacement  (Mindlin,  1949);  and  a  relative  oblique  displacement 
between  two  elastic  spheres  (Walton,  1978).  Digby  (1981)  gives  a  solution  of  the  contact 
problem  for  two  elastic  spheres  initially  bonded  over  a  finite  area. 

A  problem  that  is  relevant  to  geologic  granular  materials  is  the  contact  between  two 
elastic  spheres  cemented  by  a  different  elastic  material.  Dvorkin  et  al.  (1991)  showed  that 
the  cement  layer  may  be  approximately  treated  as  an  elastic  foundation,  and  used  this  result 
to  solve  for  the  stress  distribution  and  contact  force-displacement  relations  for  cemented 
grains  (Dvorkin  et  al.,  1994).  These  solutions  did  not  account  for  any  geometrical 
nonlinearity  beeause  the  area  of  a  cemented  contact  was  constant  and  independent  of 
pressure. 

Once  a  two-sphere  contact  problem  is  solved,  the  effective  elastic  constants  of  a  dense 
random  pack  of  identical  spheres  can  be  related  to  the  average  number  of  contacts  per 
sphere,  the  pack's  porosity,  and  the  normal  and  tangential  contact  stiffnesses  (e.g.,  Digby, 
1981;  Walton,  1987).  These  relations  are  based  on  a  statistical  approach  assuming  that  no 
relative  roll  or  torsion  occurs  between  any  pair  of  spheres  in  contact,  and  that  the 
displacements  of  the  sphere  centers  are  consistent  with  the  uniform  displacement  field 
applied  to  the  pack. 

In  this  paper  we  give  an  approximate  analytical  pressure-dependent  solution  for  the 
contact  problem  of  two  cemented  spheres.  Specifically,  we  assume  that  the  spheres  are 
uniformly  coated  with  elastic  cement  layers  (different  from  the  sphere  material),  which 
intersect  over  a  finite  contact  area  thus  bonding  the  spheres.  As  the  spheres  are  axially 
pressed  together,  the  contact  area  increases  resulting  in  increasing  normal  and  shear 
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stiffnesses.  We  use  this  solution  to  calculate  the  pressure  dependence  of  the  effective 
elastic  moduli  of  a  cemented  pack. 

The  solution  shows  that  this  pressure  dependence  may  be  large  if  the  contact  radius  is 
small  and  becomes  insignificant  for  large  contact  radii.  If  the  spheres  are  in  direct  contact 
and  the  contact  radius  is  small,  the  elastic  constants  of  the  pack  are  hardly  affected  by  the 
elastic  properties  of  the  cement.  In  this  case  Digby's  (1981)  solution  can  be  used  as  a  close 
approximation.  However,  if  two  contacting  spheres  are  separated  by  soft  cement,  its 
properties  strongly  affect  the  elastic  constants  of  the  pack. 

2.  TWO  COATED  SPHERES 


2.1.  Axial  Loading 


Consider  two  identical  elastic  spheres  of  radius  R  separated  by  a  gap  2e  and  uniformly 
coated  by  an  elastic  cement  (Figure  1).  The  contact  area  is  the  circular  cross-section  of  the 
cement  coatings  at  the  mid  plane  A-A.  The  radial  coordinate  originating  at  the  center  of  the 
contact  is  r,  and  the  initial  radius  of  the  contact  area  is  r  =  a.  When  the  spheres  are  axially 
pressed  together  the  radius  of  the  contact  area  increases  to  r  =  b.  Due  to  this  axial 
compaction,  the  relative  displacement  of  the  mid  plain  with  respect  to  the  sphere  center  is 
5^.  The  resulting  displacement  u^{r)  of  the  surface  of  a  sphere  relative  to  its  center  is 
given  by  (Johnson,  1985): 


M^(r)  = 


.  t=b  6=271 

ML 

2m, to 


tdd. 


(1) 


where  v  and  ji  are  Poisson's  ratio  and  the  shear  modulus  of  the  elastic  spheres, 
respectively;  q^(r)  is  the  normal  axisymmetric  contact  traction  on  the  surface  of  the 
spheres;  and  j  is  the  distance  between  t  and  r  (Figure  2): 


s(t,r)  =  ^|t^  +  r^  —  Itrcosd.  (2) 

The  angular  integration  in  equation  (1)  can  be  solved  analytically  (Gradshteyn  and 
Ryzhik,  1980): 


M^(r)  = 


1-  V 
Kfl 


J 

^=f 


where 


/*(x)  = 


\2K{x\ 

X  X 


X  <  1; 
X  >  1; 


and  K{k)  is  the  complete  elliptic  integral  of  the  first  kind: 


nil 

Kik)^  J 

0 


dd 


(3) 


(4) 


(5) 
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with  the  series  representation 


A  different  form  of  integral  equation  (3)  appears  in  Johnson  (1985).  In  Appendix  we 
show  that  the  form  used  here  is  computationdly  more  efficient. 


1  1 

r  ^ 

L 

n 

is  1 

n 

Figure  1  (left).  The  contact  between  two  elastic  spheres  coated  by  elastic  cement.  Figure  2  (right). 
The  coordinate  system  of  the  contact. 

Dvorkin  et  al.  (1991)  show  that  the  contact  stress  in  the  cement  can  be  approximated  by 
the  elastic  foundation  model: 


2jtf^(l-  v^)u^{r)  +  g{r)-8^ 
l-2v,  h{r)-g{r) 


where  and  are  Poisson's  ratio  and  the  shear  modulus  of  the  elastic  cement, 
respectively;  h(r)  is  half  the  distance  between  the  surfaces  of  the  spheres  (without  the 
coating);  and  g(r)  is  half  the  distance  between  the  external  surfaces  of  the  coated  layers: 


/i(r)  =  e  +  l-Vl-r 


1  +  e  -  J(l  +  ef  +a^  -r^ ,  r  >  a; 


Notice  that  at  r  =  we  have 
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S,=u^(b)  +  g(b). 

By  substituting  equations  (7)  and  (10)  into  equation  (3)  we  arrive  at 


(10) 


.  s  ^  _ A'f  ^At)  +  8(0-u,{b)-g(b) 
^  ^,lo  Kt)-g(t) 


(11) 


where 

2/i^(l~  ^c)(l~  (12) 

M1-2V,) 

We  solve  equation  (11)  for  a  given  b  to  find  M,(r).  Then  we  use  equations  (7)  and 
(10)  to  calculate  the  resultant  contact  force  and  the  contact  normal  stiffness  defined 
by: 

r=b 

F^=  ^q^{r)2nrdr,  (13) 

r=0 

In  this  way  we  relate  the  contact  stiffness  and  the  new  radius  of  contact  to  the  applied  axial 
force. 


2.2.  Tangential  Loading 

Let  us  now  consider  a  small  relative  tangential  displacement  15^  (in  the  x  direction) 
between  the  centers  of  the  axially  pressed  spheres,  such  that  there  is  no  relative  torsion  or 
rolling  between  the  spheres.  In  this  case  the  normal  and  shear  contact  tractions  can  be 
treated  independently,  and  therefore  the  tangential  loading  does  not  change  the  shape  or  size 
of  the  contact  area.  We  also  assume  that  the  free  surface  of  the  cement  is  infinitely  rough 
and  that  no  slip  occurs.  Finally,  we  use  approximations  suggested  by  Dvorkin  et  al. 
(1994):  (1)  both  the  tangential  displacement  of  the  surface  of  the  sphere  relative  to  the 
sphere's  center,  and  the  resulting  tangential  traction  are  axisymmetric;  and  (2)  the 
induced  surface  displacement  and  traction  in  the  y  direction  are  relatively  small  and  may 
therefore  be  ignored.  Note  that  for  the  special  case  of  v  =  =  0  these  two  assumptions 

are  strictly  correct. 

Accordingly,  the  displacement  due  to  the  stress  q^  is  given  by  (Johnson,  1985) 


u(r)  = 


t=b  6=271 

jdt  j 


27CM,io  eio 


{1  -  V  +  v[ 


r  —  t  cos  6 


f}tde. 


(15) 


where  s  is  defined  by  equation  (2). 

The  angular  integration  in  equation  (15)  can  be  solved  analytically  (Gradshteyn  and 
Ryzhik,  1980): 
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where  /  *  is  given  by  equation  (4), 


\[Kix)-E(x)l 

X 

2  1  1 


lx  X  X 


a:  <  1; 

X  >  1; 


and  E(k)  is  the  complete  elliptic  integral  of  the  second  kind, 


(17) 


tc/2 


E(k)-  sin^  ddO, 

0 

(18) 

with  the  following  series  representation: 

2  2^  2^-4^  2"n\  2n-l 

(19) 

From  the  elastic  foundation  approximation  (Dvorkin  et  al., 
contact  stress  q^(r)  is: 

1991)  we  find  that  the 

9x(0=  l^c  .V 

h{r)  -  g(r) 

(20) 

where  h(r)  and  g(r)  are  defined  by  equations  (8)  and  (9). 

By  substituting  formula  (20)  into  equation  (16)  we  arrive  at 

M,(r)  =  -^T  *  (-) -  V/  *  *{-)]dt. 

n/l  h(r)  -  g(r)  t  t 

(21) 

Once  integral  equation  (21)  is  solved,  equation  (20)  is  used  to  calculate  the  resultant 
contact  force  and  the  contact  tangential  stiffness  defined  by 

r^b 

Fx=  j  q^(r)27trdr, 

r=0 

(22) 

II 

(23) 

2.3.  Numerical  Procedure 


In  order  to  solve  integral  equations  (11)  and  (21)  we  uniformly  divide  domain 
0<r<^?  into  small  intervals.  Within  each  of  these  intervals  u^(r)  and  u^(r)  are 


approximated  by  quadratic  functions.  The  integration  in  every  interval  is  carried  out  by  a 
four-point  Gauss-Legendre  quadrature.  The  error  in  and  is  evaluated  as  a 
function  of  the  interval  size  £  in  the  form 


\\el=c£\  (24) 

where  c  and  a  are  constants.  For  the  example  problems  considered  in  this  work,  it  was 
found  that  the  order  of  convergence  is  a  ~  1.1  regardless  of  high  gradients  in  the  solution. 

3.  EFFECTIVE  MEDIUM  MODEL 

We  use  the  solutions  obtained  in  the  previous  section  to  calculate  the  pressure 
dependence  of  the  wave  velocities  of  a  random  pack  of  identical  cemented  spheres.  We 
assume  that  the  wavelength  is  sufficiently  larger  than  the  grain  size,  so  that  the  material  can 
be  considered  continuum. 

The  normal  contact  forces  between  the  particles  can  be  related  to  the  macroscopic 
pressure  applied  to  the  aggregate  (Digby,  1981): 


P  = 


C(l-(l>) 
47rR\l  +  e/  Rf 


(25) 


The  effective  bulk  modulus  and  the  effective  shear  modulus  can  be 

presented  as  functions  of  the  pack's  porosity  (p,  the  coordination  number  C  (the  average 
number  of  contacts  per  sphere),  the  radius  R  of  a  sphere,  and  the  normal  and  tangential 
stiffnesses  of  a  two-sphere  combination  (Digby,  1981).  Then  the  compressional-wave 
velocity  (V^)  and  the  shear- wave  velocity  ( )  are 


C(l-«  2 

^  4kR(1  +  £/R)P  ^  3 

(26) 

C(l-(/>)  (/)+3 

^  207tR(l  + e  /  R)p  "  2 

(27) 

where  p  is  the  density  of  the  aggregate,  and  and  are  given  by  (14)  and  (23). 

It  is  important  to  mention  that  equations  (25),  (26)  and  (27)  are  based  on  the 
displacement  affinity  assumption,  i.e.,  that  the  displacements  of  the  sphere  centers  are 
consistent  with  the  uniform  displacement  field  applied  to  the  pack. 

4.  EXAMPLES 

Some  limitations  of  our  theoretical  solution  are  related  to  the  size  of  the  initial  contact 
radius.  Specifically,  if  this  radius  is  large  the  spheres  may  not  be  treated  as  elastic  half 
spaces  and  the  cement  may  not  be  treated  as  an  elastic  foundation.  In  addition  we  find  from 
our  solutions  that  for  the  case  of  large  initial  radii,  yielding  of  the  cement  may  be  expected 
to  occur  before  any  noticeable  pressure  dependence  can  be  observed.  This  is  why  in  all 
presented  examples  the  initial  cemented  contact  radii  are  small  with  respect  to  the  grain 
radius. 

We  use  the  above  solutions  to  model  the  elastic  wave  velocities  in  cemented  sand.  To 
do  so  we  assume  that  the  grains  are  identical  spheres  arranged  in  a  random  dense  pack  with 
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porosity  (j)  =  0.36  and  coordination  number  C  =  9.  There  is  experimental  evidence  that 
these  assumptions  are  valid  when  estimating  the  effective  elastic  properties  of  well-sorted 
cemented  sand  (Jacoby  et  ah,  1995;  and  Dvorkin  and  Nur,  1995). 

We  fix  Poisson's  ratio  of  the  spheres  and  the  cement  at  v  =  =  0.2  and  compute  the 

wave  velocities  for  the  ratios  [i/  =0.1,  1.0,  and  10.0.  To  explore  the  effect  of  initial 

grain  separation  we  consider  two  cases  of  the  initial  gap  size:  e!  R  =  0  and  e  /  /?  =  0.01 . 

4.1.  Contact  Stresses 


The  distributions  of  normal  contact  stresses  are  given  in  Figure  3  for  increasing  (due  to 
pressure)  contact  radii.  In  all  cases  the  contact  radius  is  initially  rlR  =  0.05,  and  it 
increases  up  to  three  times  its  initial  value  as  the  external  pressure  increases.  When  the 
cement  is  as  stiff  as  the  grain  material,  or  stiffer,  the  stress  distributions  are  nearly 
independent  of  the  cement-to-grain  stiffness  ratio,  or  on  the  gap  size  (Figure  3b,  c,  e,  and 
f).  On  the  other  hand,  when  the  cement  is  significantly  softer  than  the  grain  material 
(Figure  3a  and  d),  stress  distribution  curves  are  different  from  those  with  the  stiff  cement. 
Also,  a  non-zero  gap  yields  smaller  normal  contact  stresses  for  the  same  contact  radii. 

The  distributions  of  shear  contact  stresses  are  given  in  Figure  4  for  increasing  (due  to 
pressure)  contact  radii.  When  the  cement  is  as  stiff  as  the  grain  material  (Figure  4b  and  e), 
or  stiffer  (Figure  4c  and  f),  shear  stress  concentrations  are  well  pronounced  at  the 
circumference  of  the  contact  region.  These  stress  concentrations  resemble  singularities  at 
the  edge  of  a  rigid  punch  on  an  elastic  half-space  (Johnson,  1985).  The  gap  acts  to  reduce 
these  stress  concentrations.  In  the  soft  cement  case  (Figure  4a  and  d)  we  observe  stress 
concentrations  at  the  center  of  the  contaet,  similar  to  those  observed  by  Dvorkin  et  al. 
(1994). 


4.2.  Pressure-Sensitivity  of  Velocities 

If  the  cement  and  the  grain  materials  are  identical,  our  results  (Figure  5a)  are  in  close 
agreement  with  those  of  Digby  (1981),  as  expected.  It  is  apparent  that  if  no  gap  exists 
between  the  grains,  and  the  initial  contact  radius  is  small,  the  properties  of  the  cement 
material  have  little  effect  on  the  elastic  wave  velocities  in  a  cemented  pack.  In  this  case 
Digby's  solution  is  a  good  approximation.  In  contrast,  even  a  small  gap  between  the  grains 
may  affect  the  velocities  if  the  cement  is  relatively  soft  (Figure  5b).  At  the  same  time  if  the 
cement  is  stiff,  its  properties,  as  well  as  the  gap  size,  hardly  affect  the  velocities. 

It  is  important  to  note  that  when  the  cement  stiffness  is  equal  to  that  of  the  grain,  the 
results  should  be  identical  to  those  of  Digby  -  the  gap  size  should  not  affect  the  velocities. 
However,  there  is  a  disparity  between  the  results  shown  in  Figure  5  a  and  Figure  5b  (solid 
lines).  The  error  in  the  velocities  (about  10%)  may  be  attributed  to  the  elastic  foundation 
approximation. 

APPENDIX 

Integral  equation  (3)  appears  in  Johnson  (1985)  in  the  alternative  form 

t~~b 

u^(r)  =  [  -^q  (t)K(k)dt,  (Al) 


where 
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(t  +  rf 


(A2) 


The  difference  between  the  two  forms  is  that  for  any  given  r  and  t,  the  argument  of  K  in 
(3)  is  always  smaller  than  the  argument  of  /sT  in  (Al).  Consequently,  for  any  given  r  and 
t,  the  evaluation  of  the  integrand  of  equation  (3)  will  converge  faster  and  require  less 
computation  effort  than  the  evaluation  of  the  integrand  of  equation  (Al).  Since  the 
integration  in  t  is  carried  out  numerically  there  is  a  clear  advantage  in  using  form  (3)  of  the 
integral  equation. 
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Figure  5.  Normalized  wave  velocities  in  a  pack  of  cemented  spheres  as  a  function  of  normalized 
pressure. 
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2.9.  PLASTIC  Compaction  of  Cemented  granular 

Materials 


ABSTRACT 

We  analytically  relate  hydrostatic  stress  to  strain  in  a  random  dense  pack  of  identical 
spheres  cemented  at  their  contacts.  The  spheres  are  elastic  and  the  cement  is  perfectly 
plastic.  This  solution  for  the  sphere  pack  is  based  on  a  solution  for  the  normal  interaction 
of  two  cemented  spheres.  Initially,  the  two  spheres  touch  each  other  at  a  point.  We  show 
that  as  loading  increases  and  cement  becomes  plastic,  a  finite  (Hertzian)  direct-contact  area 
between  the  spheres  necessarily  has  to  develop  and  progress.  The  stress-strain  behavior  of 
the  pack  depends  on  the  cement's  yield  limit  and  on  the  amount  of  cement.  At  the  same 
hydrostatic  stress,  the  deformation  of  the  cemented  aggregate  is  smaller  than  that  of  the 
uncemented  one.  This  difference  becomes  large  as  the  yield  limit  increases.  We  calculate 
the  bulk  modulus  of  an  aggregate  from  the  stress-strain  curve.  In  the  plasticity  domain,  the 
bulk  modulus  of  the  cemented  aggregate  is  smaller  than  that  of  the  uncemented  one.  The 
difference  between  the  two  may  easily  reach  50  percent.  Of  course,  as  the  cement's  yield 
limit  decreases,  the  aggregate's  stress-strain  curve  and  the  bulk  modulus  approach  those  of 
the  uncemented  sphere  pack.  This  theoretical  conclusion  is  qualitatively  supported  by 
experiments  on  epoxy-cemented  glass  beads.  The  maximum  contact  stress  in  the  cemented 
aggregate  may  be  less  than  a  half  of  that  in  the  uncemented  one.  This  result  explains  an 
experiment  where  an  uncemented  glass  bead  sample  failed  at  the  hydrostatic  stress  of  50 
MPa,  whereas  an  epoxy-cemented  sample  stayed  intact. 

INTRODUCTION 

Calcite-,  quartz-,  and  clay-cemented  sands  carry  static  structural  loads  and  may  be 
subjected  to  dynamic  earthquake  or  explosive  stresses.  Asphalt  concrete  is  often  used  for 
air  strips.  Understanding  the  mechanical  behavior  of  such  natural  and  artificial  cemented 
particulate  materials  is  important  in  geomechanics  and  engineering. 

A  theory  for  calculating  the  effective  elastic  properties  of  cemented  aggregates  where 
grains  and  cement  are  elastic  was  developed  by  Dvorkin  et  al.  [1,  2].  The  results  indicate 
that  intergranular  cement  may  strongly  affect  stress  distribution  at  the  contacts  (Figure  la). 
If  cement  is  soft,  contact  stress  concentration  occurs,  as  could  be  expected,  at  the  center  of 
the  contact  (Figure  Ic).  For  stiff  cement,  however,  contact  stress  concentration  occurs  at 
the  periphery  of  a  cemented  contact  (Figure  Id).  This  conclusion  is  supported  by 
photoelastic  images  of  Sienkiewicz  et  al.  [3]. 

The  effective  elastic  moduli,  as  calculated  from  the  cementation  theory,  are  close  to 
those  measured  in  artificial  cemented  aggregates  (Figure  2).  Of  course,  particles  in  natural 
granular  materials  are  rarely  spheres.  Nevertheless,  the  sphere  pack  approximation  proves 
to  be  reliable  in  calculating  the  elastic  moduli  of  natural  rocks  made  out  of  well-sorted 
grains  (Figure  3). 

Dvorkin  [4]  describes  large-strain  behavior  of  cemented  particulate  materials  by 
assuming  that  the  grains  are  elastic  and  the  cement  becomes  partly  plastic  as  local  stresses 
meet  a  plasticity  condition.  This  plasticity  condition  for  a  thin  elastic-plastic  cement  layer  is 
derived  based  on  the  von  Mises  criterion.  The  stresses  in  the  elastic  part  of  the  cement  are 
modeled  using  the  elastic  foundation  approximation.  It  is  natural  to  assume  that  in  this  case 
the  cement  is  softer  that  the  grains.  Then  the  plasticity  zone  in  the  cement  initiates  at  the 
center  of  the  contact  and  propagates  outwards  as  the  load  increases  (Figure  lb). 

In  this  paper  we  explore  the  process  of  hydrostatic  compaction  of  a  cemented  particulate 
aggregate  where  all  cement  becomes  plastic.  The  theoretical  stress-strain  law  in  such  an 
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aggregate  is  based  on  a  solution  for  the  normal  compaction  of  two  identical  elastic  spheres 
with  plastic  cement  at  their  contact. 


Coordinate  along  cement  layer  Coordinate  along  cement  layer 


Figure  1.  a.  Contact  stresses  between  two  cemented  grains,  b.  Cement  layer  with  a  plastic  central 
kernel  surrounded  by  elastic  cement,  c.  Normal  stresses  between  two  cemented  grains  with  soft 
cement,  d.  Normal  stresses  between  two  cemented  grains  with  stiff  cement.  Stresses  are  normalized 
by  the  average  stress.  Coordinate  is  normalized  by  the  radius  of  the  cement  layer. 


Porosity  Porosity 


0.3  0.4 


Porosity 


Figure  2.  Examples  of  theoretical  and  experimental  values  of  P-wave  velocity  in  cemented  granular 
materials.  Circles  are  from  the  experiments,  solid  lines  are  theoretical  predictions.  From  left  to  right: 
epoxy-cemented  glass  beads,  data  from  Yin  [5];  ice-cemented  Ottawa  sand,  data  from  Jacoby  et  al.  [6]; 
and  sintered  glass  beads,  data  from  Berge  et  al.  [7]. 


We  show  that  if  the  yield  limit  of  the  cement  is  high,  the  stress-strain  curve  noticeably 
differs  from  that  for  the  non-cemented  aggregate  with  Hertzian  grain-to-grain  interaction. 
Consequently,  the  bulk  moduli  calculated  from  the  stress-strain  curves  differ  as  well:  in  the 
plasticity  domain,  the  bulk  modulus  of  the  cemented  aggregate  is  smaller  than  that  of  the 
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uncemented  one.  The  difference  may  become  as  large  as  50  percent.  If  the  yield  limit  of 
the  cement  is  small,  the  stress-strain  curve  and  the  bulk  modulus  are  close  to  those  in  the 
uncemented  aggregate.  This  theoretical  conclusion  is  qualitatively  supported  by 
experiments  on  epoxy-cemented  glass  beads. 

We  find  that  plastic  cement  acts  to  reduce  the  contact  stresses  between  two  grains. 
This  result  explains  an  experiment  where  an  uncemented  glass  bead  sample  failed  at  the 
hydrostatic  stress  of  50  MPa,  whereas  an  epoxy-cemented  sample  stayed  intact  [8]. 


Figure  3.  P-wave  velocity  (filled  symbols)  and  S-wave  velocity  (open  symbols)  in  dry  high-porosity 
rock  samples.  Solid  lines  are  our  theoretical  predictions.  Examples  are  from  [9].  a.  Samples 
cemented  by  quartz  and  clay  (as  indicated  on  the  lines).  The  lines  are  obtained  from  the  cementation 
theory,  b.  Uncemented  samples.  Theoretical  lines  are  based  on  the  Hertz-Mindlin  contact  theory. 

FORMULATION  OF  THE  PROBLEM 

Plastic  Cement  Laver:  Geometry  and  Stresses 

Consider  two  identical  undeformed  elastic  spheres  in  contact.  The  space  around  the 
point  of  contact  is  fdled  with  cement  in  a  way  that  the  cement  layer  is  axisymmetrical  with 
respect  to  the  line  connecting  the  centers  of  the  spheres.  Let  z  be  the  coordinate  along  this 
line  with  the  origin  at  the  point  of  contact.  Then  r  is  the  radial  coordinate  in  the  plane 
perpendicular  to  the  z  axis.  The  radius  of  the  cement  layer  is  (Figure  4a). 

A  confining  force  F  is  applied  to  the  spheres  along  the  z  axis.  As  this  force  increases, 
the  cement  becomes  plastic  and  is  squeezed  sideways.  The  periphery  of  the  cement  layer 
assumes  a  new  position  at  r  =  a . 

Let  us  assume  that  the  cement  is  soft  as  compared  to  the  spheres,  and  is  linear  elastic, 
perfectly  plastic.  Its  yield  limit  is  As  force  F  increases,  a  plastic  zone  develops 
around  the  center  of  the  contact  (Figure  lb).  This  zone  propagates  outwards  [4]  and  finally 
reaches  the  free  surface  of  the  cement.  At  this  moment  of  equilibrium  which  precedes  the 
cement's  plastic  flow,  the  equations  of  balance  in  the  cement  can  be  satisfied  if  we  assume 
that  all  stresses  are  constant.  Then  from  the  von  Mises  criterion  [10]  we  conclude  that  the 
normal  stress  in  the  z  direction  is  -<7,,  and  all  other  stresses  disappear.  As  the  cement 
flows  sideways,  a  finite  direct  contact  area  may  develop  between  the  spheres.  The  radius 
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of  this  area  is  b  (Figure  4b).  If  the  loading  now  is  kept  constant,  the  plastic  flow  stops  and 
the  cement  comes  to  equilibrium  with 


Figure  4.  Cemented  contact  between  two  spheres. 


Deformation  of  the  Sphere 

We  adopt  a  common  contact  mechanics  approximation  that  if  the  cemented  contact  is 
small  (as  compared  to  the  sphere),  the  deformation  of  the  sphere  can  be  approximately 
calculated  as  that  of  an  elastic  half-space  subjected  to  the  contact  stresses  on  its  surface 
(Figure  4c).  Then  the  displacement  of  the  sphere's  surface  w(r)  due  to  axisymmetrical 
pressure  p{r)  is  [11] 


w(r)  = 


2(1 -V) 
nG 


4tr 

(t  +  rf 


]dt. 


(2) 


where  G  and  v  are  the  shear  modulus  and  Poisson's  ratio  of  the  sphere,  respectively;  and 
K(r)  is  the  complete  elliptic  integral  of  the  first  kind. 

In  the  direct  contact  area  r<b  the  displacement  of  the  sphere's  surface  satisfies  the 
Hertzian  condition 


w(r)  =  w(0)-^, 

2/v 


(3) 


where  R  is  the  sphere's  radius. 

Governing  Equations 

Equations  (1),  (2),  and  (3)  can  be  combined  to  obtain  the  following  governing 
equations  for  the  pressure  acting  on  the  sphere's  surface: 


2(1- V) 
kG 
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t  +  r 
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Atr 


{t  +  rf 


p{r)  =  <7^,  b<r<a;  p{r)  =  0,  r>a. 


r<b\ 


(4) 
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Contact  Force  and  Hydrostatic  Pressure 


If  hydrostatic  pressure  P  is  applied  to  a  dry  random  pack  of  identical  spheres  of 
porosity  (j) ,  the  stress  exerted  on  the  solid  phase  is  P  /  (1  -  0) .  Then  the  total  force  acting 
on  the  surface  of  a  sphere  is  47tR^P / (I- (j)).  This  force  is  evenly  divided  among  C 
contacts  (the  coordination  number  C  is  assumed  to  be  the  same  for  each  sphere). 
Therefore,  contact  force  F  is  related  to  hydrostatic  pressure  P  as 


F  =  (5) 

C(1 -(/>)■ 

The  porosity  of  a  random  pack  of  identical  spheres  is  about  0.36,  and  the  coordination 
number  is  between  8  and  9  (e.g.,  [12]). 

SOLUTION  OF  THE  PROBLEM 

Direct  Contact  Area  Development 

Consider  the  case  where  the  cement  is  in  plastic  equilibrium  and  b  =  0.  The  contact 
pressure  acting  on  the  grain  surface  is  cr^.  For  such  constant  pressure,  the  displacement  of 
the  grain  surface  is  [1 1] 


w(r)  = 


2(l-v)(T,a  £ 

TtG  a 


r<a, 


(6) 


where  E(r)  is  the  complete  elliptic  integral  of  the  second  kind.  Then  the  half-thickness  of 
the  cement  layer  (Figure  4a)  is 


^  +  vv(r)  -  w(0)  =  ^  +  2a^^[E(I)  - 1]. 


(7) 


The  half- volume  of  the  cement  at  a  contact  is 


V  =  j27trhir)dr  =  —  +  {--—)— — - .  (8) 

Provided  that  cement  is  incompressible,  this  volume  is  the  same  as  the  initial  volume: 
^Tca^  2  4(l-v)gy 

4R  4R  ^3  4  G 

From  Equation  (9)  and  for  a,,  /  a  «  1  we  find  that 


=  H-- 


Clr. 


4-3A 
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Gan 


(10) 


The  corresponding  contact  force  F  is  and  the  hydrostatic  pressure,  as  given  by 

formula  (5)  is 
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(11) 


^,  =  0-,(l+ 
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If  hydrostatic  pressure  exceeds  ,  a  finite  zone  of  direct  contact  between  the  grains 
necessarily  has  to  develop. 


Complete  Solution 

After  a  finite  direct  contact  zone  between  the  spheres  has  developed,  pressure 
distribution  on  the  sphere's  surface  can  be  found  from  Equations  (4).  By  using  Equation 
(6)  we  can  transform  the  integral  in  Equations  (4)  as 


2(1 -V)} 
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The  complete  elliptic  integral  of  the  second  kind  has  the  following  series  representation 
[13]: 


2^  64  a"  256  a"  "" 


By  calculating  a  polynomial  approximation  of  E(r  /  a) 


XT  ^  n  1  ^ 

^2(-)  =  :r(l-T^). 


a 


4  a 


(13) 


(14) 


we  find  that  the  relative  error  (E-E2)/  E  does  not  exceed  5  percent  if  r  /  a  <  0.4  (Figure 
5).  Therefore,  approximation  (14)  is  accurate  as  long  as  b/  a<0A. 


Figure  5.  Relative  error  of  using  approximation  (14)  for  the  complete  elliptic  integral  of  the  second 
kind. 


Now  Equations  (4)  can  be  presented  as 
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TtG  r  2  t  +  r  \t  +  r) 


j]}dt  =  —  -—,  r<b- 
^  2R  2R, 


f(r)  =  0,  r>b; 


where 


f{r)  =  p{r)-o^,  r<a;  f{r)  =  0,  r>a; 


/Vi  -  .  N  ^ 

(1  -  V)(T, 

Equations  (15)  can  be  satisfied  by  the  Hertzian  pressure  distribution  [1 1].  Therefore, 
we  arrive  at  the  following  solution  of  the  problem: 


AG-^b^  -r^ 


■(l-e)  +  cT,,  r<b-. 


^  kR{X-v) 

p{r)  =  a.,  b<r<a\  p{r)  =  Q,  r>a; 


where 

^_a-v)a,R  (,5) 

2G  Cl 

The  sphere's  surface  displacement  due  to  pressure  given  by  Equation  (18)  is  for 
b<r<a  [11]: 


w(r)  = - [{2b^  -r^)arcsin(— )  +  Zjr-.|l-(— )^]  +  e — —E{—),  (20) 

kR  r  \  r  kR  a 

and  the  half-thickness  of  the  cement  layer  is 
A(r)  =  — 

^  (21) 

+-i — -[(2b'^  - ) arcsin(-)  +  hrJl-(-f]+  s^^E(—). 

TtR  r  \  r  kR  a 

Stress-Strain  Law 

The  total  force  F  at  a  contact  can  be  calculated  from  Equations  (18)  and  then  related  to 
the  hydrostatic  pressure  P  using  Equation  (5).  The  resulting  force-balance  equation  is: 

3^_3^[  AR^P - (22) 

4(1 -£)  Ca^{\- 
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where  fi-b/ a. 

The  cement  volume  can  be  calculated  by  integrating  the  half-thickness  of  the  cement 
layer,  as  given  by  Equation  (21).  By  assuming  that  the  cement  is  incompressible  and  its 
current  volume  is  the  same  as  its  initid  volume,  we  arrive  at  the  following  equation: 


a 

+^-^[2P(l  -t- -  37tP‘^  -  2(1  -  4j8^)arcsin()3)]  (23) 

7t 

+^[2  -  (i+p^)E(p) + 

3k 


Equations  (22)  and  (23)  can  be  used  to  find  a  and  ^  at  a  given  hydrostatic  pressure  P 
exerted  on  the  aggregate.  Once  a  and  h  are  known,  we  can  calculate  the  sphere's 
deformation  e  as 


w(0) 

R 


(24) 


Finally,  the  aggregate's  bulk  modulus  is 


3  d£ 


EXAMPLES 


(25) 


Stresses  and  Strains 


Let  us  assume  that  the  spheres  are  glass  with  G  =  26.2  GPa  and  v  =  0.277.  First 
consider  a  case  where  the  amount  of  cement  is  small  {a!  R  =  0.1).  The  results  of 
calculating  strain  in  the  aggregate  versus  hydrostatic  stress  are  given  in  Figure  6a.  All 
curves  initiate  at  pressure  P^,  as  given  by  Formula  (11),  and  terminate  when  the  ratio  b!  a 

becomes  large  and  the  approximation  for  E,  as  given  by  Formula  (14),  is  not  valid.  The 
corresponding  stress-strain  curves  are  given  in  Figure  6b.  The  results  for  a  case  where  the 
amount  of  cement  is  relatively  high  (al  R=  0.3)  are  given  in  Figure  7. 

At  the  same  stress,  the  deformation  of  a  cemented  aggregate  is  smaller  than  that  of  the 
uncemented  one.  As  the  yield  limit  of  the  cement  becomes  smaller,  the  strain-stress  curves 
approach  the  curve  for  the  uncemented  aggregate.  At  the  same  yield  limit  and  at  a  given 
pressure,  the  smaller  the  amount  of  cement,  the  larger  the  deformation  (Figure  8). 

Bulk  Modulus 


The  bulk  modulus  of  the  glass  bead  pack,  as  given  by  Formula  (25),  is  calculated  for 
the  case  al  R  =  0.3  for  varying  yield  limit  (Figure  9).  The  bulk  modulus  strongly  depends 
on  the  cement  yield  limit  and  may  significantly  deviate  from  the  modulus  of  the  uncemented 
aggregate. 

Notice  that  the  bulk  modulus  decreases  with  the  increasing  cement  yield  limit.  The 
reason  is  that  at  a  given  pressure  the  radius  b  of  the  direct  contact  area  between  the  spheres 
decreases  with  the  increasing  yield  limit  (Figure  10).  The  normal  compliance  of  a  two- 
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Deformation  c,  _  Deformation 


Sphere  combination  depends  mostly  on  this  radius  and,  therefore,  increases  with  the 
increasing  yield  limit. 


a 


0.01 


0.006 


0.002 

0 


Figure  6.  a.  Strain-stress  curves  for  random  packs  of  identical  glass  spheres  with  plastic  cement. 
a!  R  =  0.1.  The  curves  are  marked  with  the  values  of  the  cement  yield  limit  in  GPa.  The 
enveloping  curve  marked  "H"  is  from  the  Hertzian  solution  for  the  uncemented  aggregate,  b.  The 
corresponding  stress-strain  curves. 


Figure  7.  Same  as  Figure  6,  but  tor  a!  R  -  0.3. 

The  moduli  calculated  from  stress-strain  curves  are  called  "static"  moduli,  as  opposed 
to  "dynamic"  moduli  which  are  calculated  from  the  P-  and  S-wave  velocities  and  density  of 


a  material.  We  measured  the  static  bulk  modulus  for  two  epoxy-cemented  glass  bead 
samples.  The  epoxy  saturation  of  the  pore  space  in  the  first  sample  was  0.15,  which 
translates  into  the  a  /  i?  value  being  between  0.24  and  0.48.  The  epoxy  saturation  of  the 
second  sample  was  0.25.  Accordingly,  the  a!  R  value  is  between  0.3  and  0.55.  Of 
course,  epoxy  is  not  a  linear  elastic,  perfectly  plastic  material  -  it  exhibits  strong 
viscoelasticity.  The  apparent  plasticity  of  epoxy  strongly  depends  on  the  strain  rate. 
Nevertheless,  it  is  interesting  to  observe  that  the  static  bulk  moduli  of  the  epoxy-cemented 
samples  have  the  same  order  of  magnitude  as  the  theoretical  bulk  moduli  for  the  pack  of 
uncemented  glass  beads,  and  for  a  cemented  pack  (Figure  11).  This  result  qualitatively 
supports  our  theoretical  solution. 


Figure  8.  Strain  versus  stress  for  random  packs  of  identical  glass  spheres  with  plastic  cement.  The 
cement  yield  limit  is  constant  at  0.1  GPa.  The  amount  of  cement  varies  as  shown  on  the  graph. 
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Figure  9.  Bulk  modulus  versus  stress  for  random  packs  of  identical  glass  spheres  with  plastic  cement 
for  a  /  R  =  0.3  and  varying  yield  limit.  The  curves  are  marked  with  the  values  of  the  cement  yield 
limit  in  GPa.  The  enveloping  curve  marked  "H”  is  from  the  Hertzian  solution  for  the  uncemented 
aggregate. 
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Figure  10.  The  b!  R  ratio  versus  hydrostatic  stress  for  uncemented  grains  (curve  "H")  and  for 
cemented  grains  with  a/  R  =  03  and  the  yield  limit  0.5  GPa  and  1  GPa  (as  marked  on  the  curves). 


Figure  11.  Bulk  modulus  versus  stress  for  two  epoxy-cemented  glass  head  samples.  The  epoxy 
saturation  of  the  pore  space  is  0.15  (filled  symbols)  and  0.25  (open  symbols).  The  theoretical  curve 
marked  "H"  is  from  the  Hertzian  solution  for  the  uncemented  aggregate.  The  theoretical  curve  marked 
"0.1"  is  for  a  cemented  aggregate  with  al  R  =  03  and  the  cement  yield  limit  of  0.1  GPa. 

INTEGRITY  OF  GRAINS 

The  maximum  pressure  on  the  sphere  surface  is  at  the  center  of  the  contact  area.  It  is 


_  4G(l-e)  b 

for  a  cemented  aggregate  and 


_  4G  bj,  _  I  3;r(l-v)~: 
7r(l-v)R’  "  ^2C(l-</>)G 


(26) 


(27) 


no 


for  the  uncemented  aggregate. 

We  calculate  for  the  uncemented  aggregate  as  well  as  for  the  case  where  a/  b  = 
0.4,  and  cr^  =0.1  GPa.  The  maximum  normal  stress  between  two  uncemented  spheres 
may  be  twice  that  in  the  cemented  case  (Figure  12). 
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Figure  12.  Maximum  contact  pressure  on  the  sphere  surface  versus  hydrostatic  stress  for  uncemented 
glass  spheres  (curve  ”H")  and  for  cemented  spheres  (curve  "0.1")- 

This  conclusion  may  explain  the  results  of  an  experiment  where  two  glass  bead  samples 
were  hydrostatically  compressed  [8].  One  sample  was  uncemented,  whereas  the  other  was 
cemented  with  epoxy.  The  volume  of  the  epoxy  accounted  for  10  percent  of  the  pore  space 
volume  which  translates  into  the  a/  R  value  being  between  0.2  and  0.43.  The  samples 
were  filled  with  water  kept  at  atmospheric  pressure.  Porosity  versus  hydrostatic  stress  was 
measured  by  the  volume  of  expelled  fluid. 

In  the  uncemented  case,  a  sharp  porosity  decrease  was  observed  at  about  50  MPa 
(Figure  13).  This  decrease  is  associated  with  the  crushing  of  grains.  The  cemented  grains 
did  not  crush.  The  photos  showed  that  in  the  latter  case,  the  grains  stayed  intact  with  the 
failure  being  loealized  within  the  epoxy. 
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Figure  13.  Porosity  versus  hydrostatic  confining  stress  in  two  water-saturated  glass  bead  samples. 
Open  symbols  are  for  the  uncemented  sample,  closed  symbols  are  for  the  cemented  one. 


Ill 


CONCLUSION 


•  Plastic  cement  has  to  be  squeezed  sideways  from  the  contact  area  between  two  elastic 
spheres  in  order  to  balance  a  confining  force.  A  finite  direct-contact  area  between  the 
spheres  necessarily  has  to  develop. 

•  The  static  bulk  modulus  of  a  cemented  aggregate  with  plastic  cement  is  smaller  than  that 
of  the  uncemented  one  (at  a  given  confining  pressure).  This  modulus  decreases  with 
the  increasing  cement  yield  limit.  The  reason  is  that  the  higher  the  yield  limit,  the 
smaller  the  direct  contact  area  between  the  spheres. 

•  The  maximum  contact  stress  in  the  cemented  aggregate  may  be  much  smaller  than  that 
in  the  uncemented  one.  This  effect  may  contribute  to  the  increased  strength  of  the 
grains  observed  in  a  cemented  material. 
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